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19.  ABSTRACT  (Gont'd) 


to  shear  banding  existed.  Preliminary  armor  penetration  calculations  using  the 
SHEARS  model  were  performed  to  simulate  this  BEIL  experiment,  and  preliminary  results 
indicate  that  a  phenomenologically  based  model,  describing  fragmentation  as  the 
nucleation,  growth,  and  coalescence  of  shear  bands,  can  be  used  to  predict  a 
fragmentation  envlroment. 

The  phenomenological  experiments  indicated  the  need  for  a  set  of  specific 
computational  features  for  armor  penetration  simulations,  such  as  finite-element 
numbering,  and  intersecting  and  wandering  slide  lines.  Because  these  features  were 
not  all  contained  in  one  existing  program,  a  program  called  Composite  HEMP  was  written 
and  ass^bled,  incorporating  many  of  those  features.  Testing  of  this  code  continues. 


SUMMARY 


Goal  and  Approach 

The  goal  of  this  three-year  program  is  to  develop,  in  cooperation 
with  the  IJ.S.  Army  Ballistic  Research  Laboratory  (BRL)  and  the  U.S. 

Army  Materials  and  Mechanics  Research  Center  (AMMRC) ,  a  computational 
method  for  predicting  armor  penetration  by  projectile  Impact,  and  in 
particular,  the  back-of-the-armor  fragment  size  and  velocity  distributions. 
An  example  of  the  extensive  fragmentation  produced  at  high  Impact  velocities 
is  illustrated  in  Figure  SI,  which  shows  three  consecutive  flash  x-ray 
exposures  of  a  steel  rod  penetrating  a  steel  armor  place. 

The  research  has  focused  on  long  rod  steel  and  depleted  uranium  (DU) 
penecracors  attacking  rolled  homogeneous  steel  armor  (RHA)  places.  The 
ultimate  goal  of  the  program  is  to  be  able  to  computationally  simulate 
normal  and  oblique  rod  iaqjacts  of  both  single  and  multiple  target  places; 
however,  the  work  to  date  has  concentrated  on  single  target  place 
configurations. 

The  computations  must  be  based  on  chose  material  properties  chat 
govern  the  fracture  and  fragmentation;  only  by  basing  the  computational 
material  models  on  the  actual  microscructural  failure  modes  and  associated 
material  properties  can  the  experinencal  results  be  extrapolated  to 
arbitrary  geometries  and  designs.  Therefore,  the  approach  comprises 
the  following  steps; 

(1)  Perform  penetration  phenomenology  experiments  to  establish 
the  key  operative  microscructural  failure  modes. 

(2)  Perform  additional  experiments  and  analyses  to  measure  the 
kinetics  of  the  pertinent  failure  modes. 

(3)  Develop  computational  models  of  those  kinetics. 

(4)  Calibrate  the  models  with  laboratory  experimental  data. 

(5)  Modify  and  adapt  a  computer  code  to  handle  the  failure 
kinetics  models. 
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Both  bdckscatter  and  downrange  fragments  are  seen.  The  penetrator  travels  from  right  to  left. 


(6)  Perform  computational  simulations  of  key  penetration 
experiments  to  test  Che  resultant  capability. 

Ac  the  end  of  the  second  year  some  progress  has  been  made  on  each 
of  these  steps  and  two  key  results  have  been  obtained.  First,  the 
dominant  fragmentation  mechanism  In  both  target  and  penecracor  has 
been  found  to  be  the  nucleatlon,  growth,  and  coalescence  of  adiabatic 
shear  bands.  Therefore,  our  efforts  have  focused  on  gaining  a  quanti¬ 
tative  understanding  of  the  shear  banding  behavior  and  on  constructing 
a  computational  model  of  that  behavior. 

Second,  the  resultant  computational  model  gives  strong  promise  of 
successfully  predicting  the  back-of-che-armor  fragmentation,  as  shown 
in  Figure  S2.  In  this  figure  we  compare  the  observed  cumulative  back- 
of-che-armor  fragment  size  distribution  from  the  experiment  of  Figure  SI 
with  Che  calculated  fragment  size  distributions  for  both  penecracor 
and  target  just  before  break-out  of  the  penecracor.  Not  only  is  the 
predicted  distribution  of  the  same  order  as  chat  observed,  but  the 
predicted  total  fragment  masses  from  both  penecracor  and  target  are 
in  fair  to  good  agreement  with  observation. 

As  discussed  below  in  Che  section  on  computational  simulations, 
the  comparison  can  only  be  semiquantltative  at  this  state  for  several 
reasons.  However,  this  preliminary  comparison  is  extrenely  encouraging 
because  It  indicates  that  the  program  approach  is  working  and  that  it 
is  possible  to  predict  back-of-the-armor  fragmentation  based  on  basic, 
measurable  properties  of  the  target  and  penetrator  materials. 

The  program  status  is  summarized  in  more  detail  below. 

Phenomenology  Experiments 

A  penetration  phenomenology  data  base  has  been  generated  by  BRL 
and  SRI  for  4340  steel  and  DU  long  rod  penetrators  striking  single  RHA 
target  places  at  normal  and  oblique  impact  angles.  The  data  include 
back-of-che-armor  fragment  size  and  velocity  distributions  for  comparison 
with  computational  simulations.  The  key  ml c rose rue Cura 1  failure  mode 
Identified  for  long  rods  striking  relatively  thick  places  is  adiabatic 
shear  banding. 
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FIGURE  S2  COMPARISON  OF  MEASURED  OOWNRANGE  FRAGMENT  SIZE  DISTRIBUTION 
FROM  EXPERIMENT  SHOWN  IN  FIGURE  SI  WITH  DISTRIBUTIONS 
CALCULATED  FROM  COMPUTER  SIMULATION 
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Tha  role  played  by  Che  adiabatic  shear  bands  changes  significantly 
as  the  Itnpacc  velocity  increases.  At  velocities  below  the  ballistic 
limit,  little  shear  banding  occurs  In  the  RHA  target.  Instead,  the 
target  deforms  plastically  to  form  a  cavity  as  the  penetracor  material 
flows  around  and  back  out  the  cavity.  However,  shear  banding  does 
dominate  in  the  penetracor  deformation  and  backflow.  This  is  shown  in 
Figure  S3,  which  is  a  section  through  a  thick  steel  target  struck  by  a 
long  steel  rod.  Shear  banding  caused  the  penetracor  nose  to  divide  into 
fragments,  which  then  flowed  back  out  of  the  cavity  to  form  back-scatter 
fragmentation.  Figure  S4  shows  similar  pictures  for  thinner  targets 
at  normal  impact  velocities  Just  below  the  ballistic  limit.  No  shear 
banding  in  the  targets  is  seen,  and  the  beginning  of  the  shear-band- 
induced  penetracor  fragmentation  and  backflow  is  evident. 

For  oblique  impacts  Just  below  the  ballistic  limit,  the  situation 
is  similar  to  chat  for  normal  impact,  as  shown  in  Figure  SS.  Backflow 
of  shear-banded  fragments  from  the  penetracor  causes  a  scalloped  cavity 
boundary,  similar  to  those  seen  in  the  normal  impact  case. 

Ac  impact  velocities  Just  above  the  ballistic  limit,  the  picture 
changes  dramatically,  as  shown  in  Figure  S6.  The  plastic  flow  in  the 
target  has  now  become  localized  by  adiabatic  shear  banding  to  form  a 
target  plug.  The  shear  banding  in  the  target  is  confined  to  the  narrow 
region  chat  defines  the  plug;  the  plug  diameter  is  close  to  that  of  the 
expanded  nose  of  the  penetracor. 

As  the  iaq>act  velocities  Increase  further  to  values  well  above  the 
ballistic  limit  for  both  oblique  and  normal  In^acts,  the  picture  changes 
again.  Mow  the  penetration  velocity  is  cong)arable  to,  or  greater  than, 

Che  plastic  wave  velocities  in  the  target,  and  high  plastic  strains  and 
associated  shear  banding  are  concentrated  in  a  bow  wave  around  the 
penetracor  nose.  Thus,  extensive  shear  banding  occurs  throughout  the 
plug  and  adjacent  to  the  plug  region,  as  shown  in  Figure  S7  for  the 
experiment  depicted  in  Figure  SI.  The  plug  and  collateral  material  are 
shattered  by  shear  banding,  and  a  spray  of  target  and  penetrator  fragments 
emerges  behind  the  armor. 
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MP-7893-ai 


FIGURE  S3  CROSS  SECTION  THROUGH  A  TARGET  IMPACTED  PERPENDICULARLY 

BY  A  LONG-ROD  PENETRATOR,  SHOWING  BACKWARD  FLOW  OF  SHEAR 
BAND-ERODED  PENETRATOR  MATERIAL  (Photograph  courtesy  of  BRL) 

Note;  Horizor>tal  hole  near  bottom  was  pre-drilled  for  instrumentation  access. 
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MP-7893-121 

FIGURE  S4  CROSS  SECTIONS  THROUGH  CENTERS  OF  RHA  TARGETS 
IMPACTED  BY  BLUNT-NOSED  4340  STEEL  PROJECTILES 


IX 


I _ I 

1  cm 


MA-7a93-89 

FIGURE  So  CROSS  SECTION  OF  CRATER  FORMED  BY  IMPACT  OF  FULL-SCALE  STEEL 
LONG-ROD  PENETRATOR  INTO  4-inch-THICK  RHA  TARGET  AT  VELOCITY 
OF  ^  1 .5  km/s  AND  OBLIQUITY  OF  =s:  70° 
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(b) 


MP-7893-42A 

FIGURE  S6  BLUNT-NOSED  4340  STEEL  PROJECTILE  EMBEDDED  IN  RHA  TARGET 
PLUG  (a)  AND  CROSS  SECTION  THROUGH  IMPACT  AXIS  (h) 

FOLLOWING  IMPACT  AT  792  m/s 
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(b)  (c) 

MP-7893-t14 

FIGURE  S7  CROSS  SECTION  THROUGH  CENTER  OF  RHA  TARGET  (a)  IMPACTED 
BY  BLUNT-NOSED  4340  STEEL  PROJECTILE  AT  1.55  km/sec  AND  TWO 
DETAILED  VIEWS  (b)  AND  (c)  SHOWING  DENSE  SHEAR  BANDING  IN 
REGION  ADJACENT  TO  PENETRATION  HOLE 

Target  is  from  experiment  depicted  in  Figure  SI. 
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Experiments  co  Measure  Shear  Band  Kinetics 


After  adiabatic  shear  banding  was  identified  as  the  dominant  failure 
mode,  the  next  step  in  the  program  was  to  perform  a  series  of  experiments 
to  measure  the  shear  band  kinetics.  In  addition  to  the  usual  mechanical 
tests,  these  experiments  consisted  of  applying  internal  explosive  loads 
to  thick-walled  cylinders  of  the  test  material  so  chat  they  expanded  co 
final  dimensions  controlled  by  massive  containment  cylinders.  The  setup 
for  this  contained  fragmenting  cylinder  (CFC)  experiment  is  shown  in 
Figure  S8.  Shear  bands  originating  at  the  inner  surface  of  the  cylinder 
were  counted,  measured,  and  correlated  with  computed  stress  and  strain 
histories,  as  shown  in  Figure  S9,  co  deduce  the  shear  band  kinetics  formulae. 

This  procedure  has  now  been  completed  for  a  previously  available 
data  base  for  4340  steel  and  is  under  way  for  a  data  base  generated 
during  the  past  year  for  KHA  material.  Similar  data  for  DU  will  be 
taken  during  the  upcoming  year. 

Computational  Model  Development 

A  computational  model  of  the  observed  shear  band  kinetics  has  also 
been  completed.  The  model,  called  SHEARS,  describes  the  nucleation, 
growth,  and  coalescence  of  the  shear  bands  to  form  fragments.  SHEARS 
produces  good  co  fair  agreement  with  a  variety  of  CFC  data  for  4340  steel 
and  wHl  soon  be  calibrated  against  RHA  steel  CFC  data. 

To  fulfill  the  program  goals,  SHEARS  must  be  based  on  the  fundamental 
material  properties  that  govern  the  shear  band  kinetics.  The 
status  of  relating  measured  SHEARS  parameters  co  more  fundamental 
properties  is  as  follows.  First,  Che  shear  band  nucleation  process 
appears  at  least  partially  understood.  Localization  into  shear  bands 
occurs  in  4340  steel  at  the  value  of  equivalent  plastic  strain  predicted 
by  plastic  instability  theory,  namely,  the  strain  at  which  the  race  of 
strength  reduction  from  thermal  softening  exceeds  the  race  of  strength 
Increase  from  the  work  hardening.  Thus,  the  nucleation  threshold  strain 
depends  on  basic  material  properties,  as  desired.  However,  ocher  observed 
features  of  the  nucleation  behavior  have  not  yet  been  related  to  more 
fundamental  material  properties.  For  example,  the  shear  band  propagation 
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FIGURE  S8  SCHEMATIC  OF  CONTAINED  FRAGMENTING  CYLINDER  (CFC) 
EXPERIMENTS  FOR  STUDYING  SHEAR  BAND  KINETICS 
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threshold  is  governed  by  the  resolved  shear  strain  in  the  plane  of  the 
band  rather  than  by  the  equivalent  plastic  strain,  and  the  nucleatlon 
sites  have  not  yet  been  related  to  mlcrostructural  heterogeneities  or 
to  instability  wavelengths.  Work  Is  continuing  to  better  relate  these 
aspects  of  nucleatlon  to  ba.slc  properties. 

Second,  the  observed  shear  band  propagation  is  consistent  with  a 
viscous  growth  relation  similar  to  that  measured  previously  for  tensile 
voids  and  cracks.  However,  the  dependence  of  this  propagation  relation 
on  more  fundamental  material  properties  has  not  yet  been  established, 
and  work  is  under  way  to  gain  a  better  understanding  of  this  process. 

Finally,  as  the  shear  bands  propagate  and  coalesce  to  form  fragments, 
the  material  softens  and  the  amount  of  shear  stress  that  can  be  supported 
decreases.  SHEAR3  contains  a  description  of  this  work  softening  process, 
which  is  assumed  to  depend  on  the  frictional  forces  in  the  sliding  bands 
and  on  the  material  anisotropy.  SHEAR3  currently  describes  the  frictional 
effect,  but  not  the  anisotropy.  Work  is  therefore  continuing  to  improve 
the  model  to  handle  material  anisotropy. 

Calibration  of  the  Model 

As  mentioned  above,  the  SHEAR3  model  has  been  calibrated  with  labo¬ 
ratory  mechanical  and  CFC  data  for  43A0  steel,  and  calibration  with 
mechanical  and  CFC  data  for  RHA  is  in  progress.  During  the  upcoming 
year,  we  plan  to  perform  similar  calibration  experiments  for  DU. 

Computer  Code  Development 

A  wave  propagation  computer  code  has  been  developed  during  the  past 
year  to  be  used  in  conjunction  with  SHEAR3  for  computationally  simulating 
penetration  events.  The  code  is  a  modification  of  the  HEMP  code  and  is 
called  Composite  HEMP.  This  code,  which  is  now  in  the  debugging  stage, 
has  the  following  features.  First,  the  Input/output  is  similar  to  that 
for  the  versions  of  HEMP  now  in  use  at  BRL  and  AMMRC.  Second,  the  cell 
numbering  algorithm  has  been  changed  to  a  finite  element  scheme  used. 
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for  example,  by  the  EPIC  code.  This  scheme  has  several  advantages, 

Including  easier  inclusion  of  double  slide  lints  near  the  penetrator/ 
target  Interface.  Third,  an  automatic  rezonlng  algorithm  has  been 
included.  Fourth,  a  self-locating  and  self-orienting  internal  slide 
line  subroutine  has  been  planned,  but  has  not  yet  been  written.  Finally, 
the  NAG/FRAG  family  of  failure  models,  including  SHEAR3,  has  been  inserted 
into  the  code. 

During  the  upcoming  year.  Composite  HEMP  will  be  used  in  computational 
simulations  of  selected  penetration  experiments. 

Computational  Simulations 

Preliminary  computational  simulations  of  several  normal  impacts  of 
4340  steel  rods  on  RHA  targets  have  been  performed  during  the  past  year 
with  the  TROTT  code,  an  in-house  SRI  code  that  is  similar  to  HEMP. 

TROTT  contains  the  SHEAR3  failure  model,  but  does  not  contain  the  finite 
element  numbering  scheme  or  a  double  slide  line  algorithm.  The  purpose 
of  these  calculations  was,  first,  to  test  the  potential  of  SHEAR3,  using 
material  parameters  based  on  CFC  experiments,  for  simulating  long  rod 
impact  events,  and,  second,  to  help  design  the  experiment  shown  in 
Figures  SI  and  S7.  The  results  of  the  calculations  are  shown  in  Figures  SIO 
through  S12.  The  conditions  of  the  calculations  were  similar  to  those  of 
the  experiment  shown  in  Figures  SI  and  S7,  where  a  1. 02-cm-diameter  by 
10. 16-cm-long  4340  steel  (R^40)  rod  normally  Impacted  a  2 . 54-cm-thick  by 
19 . 24-co^diameter  RHA  target  plate  at  1550  m/s. 

Figure  SIO  shows  a  plot  of  computed  equivalent  plastic  strain  contours 
in  the  penetrator  and  target  17  ys  after  impact.  The  contours  with 
equivalent  plastic  strains  of  0.2  represent  approximately  the  onset  of 
shear  banding,  and  values  of  0.6  correspond  roughly  to  complete  fragmen¬ 
tation.  (SHEAR3  computes  different  levels  of  damage  in  four  separate 
orientations,  but  detailed  plots  of  the  damage  contours  in  these 
orientations  would  be  difficult  to  present  in  a  single  illustration.) 
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INCLUDING  CONTOURS  OF  EQUIVALENT  PLASTIC  STRAIN 


unzippiuq  slide 


Ic  can  be  seen  that  a  large  strain  gradient  develops  at  the  boundary 
of  what  can  be  considered  a  developing  plug,  but  the  plug  is  being 
fragnented  completely  by  a  fragmentation  front  that  looks  like  a  "bow 
wave"  ahead  of  the  mushrooming  pcnetrator  nose.  In  addition,  the 
penctrator  nose  itself  is  heavily  fragmented,  as  is  a  collateral  region 
outside  the  plug,  although  the  f ragnencatlon  wave  in  this  region  lags 
behind  that  inside  the  plug. 

Figures  Sll  and  S12  show  similar  plots  at  17  and  23  us  after  impact, 
respectively,  for  a  case  when  the  impact  velocity  Ls  2000  m/s.  The 
results  are  similar  to  those  of  Figure  SIO.  In  particular,  in  Figure  S12 
it  la  seen  chat  at  23  us  after  impact,  the  fragmentation  front  has  nearly 
reached  the  rear  surface.  We  would  expect  that  when  the  front  reached 
the  rear  surface,  a  spray  of  fragments  would  be  released  from  the  target 
and  penctrator.  The  total  number  of  fragments  and  their  mass  distribution 
for  both  the  penctrator  and  the  target  at  the  time  shown  in  Figure  S12 
were  shown  previously  in  Figure  S2,  along  with  the  observed  mass  distribu¬ 
tion  from  Che  experiment  of  Figure  SI. 

The  computed  mass  fraction  of  the  penctrator  that  was  completely 
fragmented  at  the  time  shown  in  Figure  S12  is  31X,  in  fair  agreement 
with  the  value  of  40%  obtained  by  measuring  the  mass  of  the  recovered 
unfragmented  distal  portion.  Furthermore,  the  computed  mass  lost  from 
the  target  in  the  form  of  the  plug  and  collaterally  fragmented  material 
is  35  grams,  in  good  agreement  with  the  34  grams  obtained  by  weighing 
the  armor  place  before  and  after  the  experiment. 

Comparison  of  computations  with  observations  in  Figure  S2  shows 
that  we  compute  coo  many  small  fragments  and  too  few  large  ones.  This 
discrepancy  can  have  several  causes.  First,  an  unknown  number  of  small 
fragments  were  not  recovered  in  the  experiment.  Second,  since  the  SKEAR3 
parameters  for  RHA  have  not  yet  been  derived  from  the  CFC  data,  we  used 
in  the  computations  the  measured  nucleatlon  threshold  strain  for  RHA 
but  the  nucleatlon  rate  function  appropriate  for  the  4340  penetracor 
steel.  On  the  basis  of  our  experience  with  another  ductile  steel  (Armco 
iron),  we  expect  the  nucleatlon  rate  to  be  less  for  the  RHA  than  for  the 
4340  steel.  Accounting  for  both  of  these  factors  would  reduce  the 
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dlffer#nc«  b«tw«en  observation  and  compucacion.  Another  cause  for 
conpuca clonal  error  Ir  che  as-yec  unmodeled  effect  of  the  significant 
anisotropy  of  che  RHA.  Furthermore,  che  calculations  of  Figures  Sll  and 
S12  used  a  somewhat  higher  Impact  velocity  than  was  achieved  In  the 
experiment. 

These  compuca cions  were  cermlnaced  because  of  excessive  cell  dis- 
corcion— a  problem  chat  Illustrates  che  need  for  che  rezoning  and  slide 
line  features  contained  in  the  Composite  HEMP  code. 

The  above  uncertainties  mean  chat  these  comparisons  of  computed 
and  observed  fragment  size  dlscributlons  can  only  be  semlquancicaclve 
Nevertheless,  the  computations  provided  the  Information  desired. 

The  overall  approach  Is  showT'  to  be  promising;  SHEAR3  parameters  obtained 
from  laboratory  CFC  experiments  appear  to  produce  encouraging  results  in 
che  penetration  corpucacions.  Ac  an  impact  velocity  well  above  che 
balllsclc  limlc,  SHEAR3  will  produce  a  spray  of  target  and  penecrator 
fragments  with  a  reasonable  mass  or  size  distribution.  Work  is  continuing 
to  perform  a  similar  computational  simulation  of  the  case  In  which  a  plug 
formed  at  a  velocity  Just  above  the  ballistic  limit  (Figure  56).  Further 
calculations  must  then  await  che  decerroinaclon  of  che  SHEAR3  parameters 
for  RHa  material. 

Future  Work 

During  che  third  year  of  che  program,  the  following  work  is  planned: 

•  Conclude  analyses  of  the  RHA  CFC  experiments  to  obtain  the 
SHEAR3  parameters  for  chat  material. 

•  Generate  CFC  and  tensile  test  data  for  DU,  and  analyze  che 
data  to  obtain  che  SHEAR3  parameters. 

•  Continue  che  effort  to  relate  che  measured  SHEAR3  parameters 
to  more  fundamental  material  properties. 

•  Continue  development  of  che  Composite  HEMP  code. 

•  Perform  Composite  HEMP  simulations  of  normal  and  oblique 
penetration  experiments. 
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1  INTRODUCTION 


The  tradiclonal  empirical  approach  for  designing  armor  and  weapons 
has  often  been  effective  in  gaining  successive  incremental  improvements 
over  the  established  state  of  the  art.  However,  as  more  complex  weapon 
and  armor  systems  are  introduced,  the  usually  very  lengthy  and  expensive 
empirical  approach  becomes  less  and  less  tractable.  In  recent  years 
numerical  finite  difference  and  finite  element  methods  have  been 
developed  to  the  point  where  they  are  capable  of  handling  complex  shapes, 
deformations,  and  material  behavior  and  therefore  promise  to  be  useful 
in  armor  and  penetrator  design. 

The  computational  approach  to  ordnance  design  is  an  iterative 
procedure  in  which  approximate  computational  simulations  of  an  ordnance 
event  are  first  performed  using  estimatod  values  of  material  properties. 
Test  firings  are  then  carried  out  to  evaluate  the  performance  of  the 
design,  to  provide  data  for  comparison  with  calculated  results,  and  to 
indicate  the  adequacy  of  the  material  property  values.  The  test  results 
then  suggest  changes  in  the  initial  ordnance  design  and  guide  the 
selection  of  material  descriptions  and  dynamic  material  property  tests 
for  the  next  series  of  computational  simulations.  Test  firings  of  a 
second  generation  design  are  then  performed,  and  the  results  are 
compared  with  computational  results,  as  before.  The  iterative  process 
is  repeated  until  the  ordnance  refinements  produce  the  desired  performance. 
This  approach,  currently  being  developed  by  the  DoD,  reduces  the  overall 
design  effort  because  the  interaction  of  experiments  and  calculation.^ 
leads  more  rapidly  to  an  acceptable  design. 

Material  response  chat  must  be  modeled  In  contpucaclonal  simulations 
of  ordnance  events  Include  a  thermodynamic  equation  of  state,  plastic 
flow,  and  material  disintegration.  Existing  descriptions  of  thermo¬ 
dynamic  equations  of  state  are  adeguate  for  most  calculations  in  ordnance 
design,  and  concepts  of  perfect  plasticity  are  adequate  for  many  ordnance 
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applications.  The  most  serious  limitation  in  efficient  ordnance  design 
is  in  the  modeling  of  material  disintegration.  Disintegration  can  occur 
by  several  mechanisms,  depending  on  the  material  and  the  conditions  of 
stress,  temperature,  and  strain  race.  Detailed  examinations  of  material 
recovered  from  many  and  varied  ordnance  events  have  identified  three 
main  material  failure  modes:  brittle  fracture,  ductile  fracture,  and 
localized  shear  failure.  Of  these  three  microscopic  failure  modes, 
the  last  two  are  more  common  In  metallic  armor. 

Although  simple  failure  criteria  may  suffice  to  obtain  reliable 
computational  simulations  of  certain  ordnance  events,  ocher  situations 
require  a  detailed  model  and  hence  a  detailed  physical  understanding 
of  nuclcacion,  growth,  and  coalescence  of  cracks,  voids,  and  shear 
bands.  Whereas  a  significant  amount  of  relevant  research  has  already 
been  done  in  this  area,  much  more  is  required  to  allow  adequate  compu¬ 
tations  to  be  performed. 

A.  Approach 

The  objective  of  this  three-year  program  is  to  develop  for  the 
U.S.  Army  Materials  and  Mechanics  Research  Center  (AMMRC)  and  Ballistic 
Research  Laboratory  (BRL)  a  computational  macerial-dlsincegraclon  model 
for  simulating  the  long-rod  penetration  of  spaced  armor.  This  model 
will  provide  the  capability  for  predicting  the  behlud-the-armor  fragment 
environment  for  given  penetrator  attacks  and  will  aid  in  armor  and 
penetrator  design. 

Our  approach  to  meet  this  objective  has  the  following  six  main 
parts : 

(1)  Materials  are  selected  for  study,  and  penetration  phenomenology 
experiments  are  performed  In  which  rods  are  fired  at  various 
velocities  and  angles  of  incidence  against  target  places. 

The  rod  and  target  specimens  arc  recovered,  sectioned,  and 
metallographlcally  examined  to  reveal  the  key  mlcrostructural 
damage  mechanisms  that  govern  failure  and  disintegration  of 
both  penetrator  and  target.  The  baseline  materials  chosen  for 
study  are  depleted  uranium  (DU}-3/4  vcZ  titanium  alloy  rods 
fired  against  rolled  homogeneous  armor  (RHA)  targets.  Addi¬ 
tional  materials  that  may  be  studied  in  the  future  are  tungsten 
alloy  rods  and  electroslag-remelc  (ESR)  steel  armor. 


(2)  Laboratory  material  property  characterization  experiments  are 
perforaed,  where  necessary,  to  measure  the  properties  that 
govern  the  damage  mechanisms.  £xperlment.s  are  chosen  to 
exercise  a  particular  failure  mechanism  In  a  geometry  that 

Is  relatively  simple  to  model.  The  contained  fragmenting 
cylinder  (CFC)  test,  for  example,  is  used  to  study  shear  band 
nucleaclon  and  growth. 

(3)  Previously  developed  computational  models  of  the  various 
damage  processes  are  modified  and  adapted  to  describe  the 
observed  failure  mechanisms.  The  models  are  all  based  on 
Che  nucleaclon  and  growth  (NAG)  approach,  which  atcempts  to 
simulate  qualitatively  and  quantitatively  the  actual  mlcro'- 
mechanlcal  failure  processes  involved. 

(4)  These  models  are  calibrated  against  the  experimental  data 
by  computationally  simulating  the  material  characterization 
experiments.  Parameters  are  obtained  that  quancltatlvely 
describe  each  material's  behavior,  for  use  in  subsequent 
armor  penetration  computations. 

(5)  A  computer  code  la  assembled,  containing  the  features  that 
have  been  seen  to  be  necessary  to  simulate  the  long-rod 
penetration  of  armor.  Although  this  code  includes  many 

of  the  same  features  found  in  the  various  HEMP  codes,  it 
iBusc  also  contain  some  new  logic  to  handle  the  complex 
phenomenology  of  armor  penetration,  and  it  must  incorporate 
the  material  models  developed  above. 

(6)  Finally,  the  ability  of  the  models  to  correctly  predict 
the  main  features  of  long-rod  penetration  of  spaced 
armor  Is  assessed,  and  if  Che  results  are  promising,  the 
models  are  used  to  predict  behlnd-che-armor  fragment 
environments  and  to  aid  in  armor  and  penetracor  design. 

It  should  be  emphasized  that  these  six  parts  to  our  approach  are  not 

performed  sequentially,  but  are  developed  simultaneously,  with  a  large 

amount  of  interaction  necessary  between  them. 


B.  Summary  of  First  Year's  Progress 

The  experimental  progress  made  during  the  first  year  of  the  program, 
as  detailed  In  the  first  annual  report,  is  shown  graphically  in  the 

Program  Technical  Milestones  chart.  RHA  material  was  obtained  from  BRL 
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early  in  the  program,  and  a  metallurgical  characterization  was  performed. 
Quasi-static  tensile  tests  were  performed  with  samples  cut  from  three 
orthogonal  directions,  and  stress-strain  curves  were  obtained  chat  showed 
the  RHA  to  be  very  similar  to  that  described  by  Benck*  and  Hauver  and  to 
exhibit  a  large  degree  of  anisotropy  in  its  ductility. 


Phenomenology  experiments  were  performed,  involving  normal-incident 
rod  impaccs  onto  1/4-inch- chick  RHA  places  at  velocicies  near  the 
ballistic  limit  (655  to  792  m/s).  Ac  these  moderate  velocicies,  the 
target  failed  by  catastrophic  shear  band-induced  plugging  following 
substantial  homogeneous  deformation,  but  without  any  significant  internal 
damage  or  fragmentation  of  Che  plug  or  the  adjacent  target  regions  (see 
Figures  S4-and  S6  of  the  Summary) . 

A  aeries  of  CFG  experiments  was  begun  using  RHA  specimens  of 
different  orientations.  The  series  was  completed  during  the  second 
year,  and  the  results  are  reported  in  Section  III. 

Finally,  an  experimental  technique  was  developed  to  enable  us  to 
determine  the  constitutive  relations,  or  flow  curve,  of  a  materia]  in 
compression  at  the  same  high  strain  races  existing  in  a  long-rod  armor 
penetration.  This  technique,  called  the  symmetric  rod  impact  test 
(or  symmetric  Taylor  test),  involves  photographing  (at  s  10^  frames/s) 
the  silhouettes  of  the  cylindrical  specimens  during  the  actual  deforma¬ 
tion  phase  following  a  high-speed  inqjact.  The  event  is  then  simulated 
by  a  two-dimensional  computer  code  (such  as  HEMP),  varying  the  flow 
parameters  until  the  computed  and  experimental  deformations  agree.  The 
feasibility  of  this  approach  was  tested  successfully  using  aluminum 
specimens . ^ 

Progress  made  in  the  computational  part  of  the  program  during  the 
first  year  included  development  of  fracture  models,  acquisition  of  a 
HEMP  program,  and  derivation  of  a  criterion  for  shear  band  nucleation. 

A  greatly  expanded  version  of  the  shear  band  model  was  developed 
and  preliminary  efforts  were  made  to  correlate  the  model  with  the  4340 
data  from  contained  fragmenting  cylinders.  Major  changes  were  in  the 
nucleation  and  growth  processes,  and  in  the  algorithm  for  stress  relax - 
tion  caused  by  the  developing  damage. 

A  ductile  fracture  model  for  intermediate  strain  rates  (durations 
of  10  to  100  usee)  was  derived  by  combining  features  from  fracture  models 
applicable  to  either  shock  or  quasi-static  loading.  The  model  was 
not  implemented  in  a  computer  subroutine. 
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Several  compucer  slmulaclons  were  made  wich  che  TROT!  code  of 
blunc-nosed  pcnetrators  impacting  steel  targets.  These  simulations 
showed  many  general  features  of  penetration  calculations,  including  the 
needs  for  rezoning  and  slide-line  capabilities.  In  attempt  to  satisfy 
those  needs,  we  acquired  a  version  of  the  two-dimensional  wave  propaga¬ 
tion  code  HEMP  and  made  arrangements  to  obtain  a  more  recent  version  of 
Che  code  with  intersecting  (or  double)  slide  lines. 

A  single  criterion  for  the  initiation  of  shear  bands  was  derived 
from  stability  considerations:  when  the  rate  of  stress  reduction  caused 
by  thermal  softening  equals  or  exceeds  the  rate  of  work-hardening,  the 
material  becomes  constitutively  unstable.  This  criterion  was  able  to 
correctly  rank  materials  according  to  their  propensity  to  form  shear 
bands . 

C.  Report  Organization 

This  report  describes  our  progress  during  che  second  year  of  the 
program  and  our  plans  for  che  third  year.  Section  II  describes  che 
phenomenological  experiments,  and  Section  III  describes  che  materials 
characterization  experiments  (parts  1  and  2  of  our  approach  as  outlined 
above).  Section  IV  describes  an  additional  experiment  performed  to 
provide  a  downrange  fragmentation  data  base  for  early  comparison  with 
our  developing  computer  predictions  capability.  This  allows  us  to  determine 
at  an  early  state  if  our  phenomenologically  based  computational  models 
can  indeed  predict  fragmentation  and  velocity  distributions. 

Development  of  the  shear  band  model  (part  3  of  the  program)  and 
calibration  of  this  model  against  che  experimental  data  (part  ^)  are 
described  in  Section  V.  Synthesis  of  a  composite  code  for  penetracor 
simulations  (part  S)  is  discussed  in  Section  VI,  and  preliminary  simula¬ 
tions  of  projectile  penetrations  (part  6)  is  described  in  Section  VII. 
Section  VIII  gives  our  plans  for  completing  this  work  during  che  third 
year  of  che  program. 


II  PHENOMENOLOGICAL  STUDIES 


During  this  past  year,  we  studied  the  phenomenology  of  armor 
penetration  by  examining  recovered  specimens  from  three  different  types 
of  ballistics  experiments:  (1)  a  full-scale,  high-velocity,  high-obliquity 
impact  of  hardened  steel  into  RHA,  stopping  Just  short  of  full  penetration, 
(2)  a  series  of  quarter-scale  impacts  of  DU  into  RHA  triple  armor  at 
various  velocities,  and  (3)  a  quarter-scale,  zero-obliquicy  impact  of 
43A0  steel  into  RHA,  well  above  the  ballistic  limit.  In  all  cases  studied, 
shear  banding  was  seen  to  be  the  dominant  failure  mechanism  of  both  the 
target  and  penetrator. 

A.  Full-Scale  Impact  of  Steel  into  RHA  at  High  Obliquity 

We  performed  a  metallographic  examination  of  sections  of  the  target 
and  penetrator  recovered  from  the  ballistics  test  witnessed  by  attendees 
of  the  February  1980  project  meeting  at  BRL.  In  this  full-scale  test, 
a  l-inch-diamecer  long-rod  steel  penetrator,  impacting  at  a  high  velocity 
(=1500  m/s)  and  a  high  angle  of  obliquity  (=70“),  stopped  just  short  of 
penetrating  a  A-inch-thlck  slab  of  RHA.  This  incomplete  penetration 
allowed  us  to  observe  some  Important  phenomenology  that  would  have  been 
difficult  or  impossible  to  see  if  complete  penetration  had  occurred.  The 
sections  we  examined  included  a  distal  2-lnch  section  of  the  penetrator 
(the  only  large  penetrator  piece  recovered)  and  a  slab  of  the  target  that 
Included  half  of  the  damaged  region.  The  following  is  a  detailed  report 
on  the  results  of  this  examination  and  the  phenomena  we  observed. 

1.  Analysis  of  the  Penetrator 

The  recovered  penetrator  piece  is  shown  in  Figure  1  along  with 
a  polished  and  etched  cross  section  through  the  center  of  the  proximal 
end.  A  hardness  traverse  across  the  specimen  yielded  values  of  8^51  +  1 
for  the  undamaged  end.  Our  past  experience  with  steels  of  this  hardness 
led  us  to  expect  shear  banding  at  relatively  low  shear  strains  (<20X), 
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MP-7893-86 

FIGURE  1  RECOVERED  SEGMENT  OP  STEEL  PENETRATOR  (a)  AND  AXIAL 
CROSS  SECTION  (b)  OF  PROXIMAL  END 
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and  this  is  precisely  what  we  sec  when  we  look  at  the  proximal  region 
under  higher  magnification  (see  Figure  2).  A  consistent  pattern  of 
shear  bands  extending  into  the  penetrator  at  angles  of  45°  +  10°  from 
the  axis  is  seen,  with  the  bands  cutting  through  material  having  under¬ 
gone  relatively  small  shear  strain  (as  shown  by  distortion  of  the 
extrusion- induced  texture).  The  maximum  shear  displacement  along  a  band 
appears  to  be  about  50  um.  Cracking  has  occurred  along  or  very  near 
a  few  of  the  bands,  and  it  is  by  separation  along  these  cracks  that  the 
penetrator  retains  its  pointed  nose  while  eroding  during  a  penetration. 
Arrows  in  Figure  2  point  to  material  surrounded  by  shear  bands;  we 
believe  this  material  would  have  been  eroded  next  if  the  penetration 
process  had  continued. 

The  photographed  cross  section  (Figure  lb)  shows  two  cracks 
that  are  apparently  not  connected  with  a  shear  band.  The  largest  of 
these,  a  crack  perpendicular  to  the  axis,  transverses  material  that  has 
undergone  negligible  shear  deformation.  The  role  of  these  tensile  cracks 
in  the  penetrator  erosion  is  not  clearly  understood.  They  may  not  even 
have  occurred  during  the  erosion  regime.  In  any  case,  the  proximal 
surface  of  the  penetrator  consists  of  material  transformed  during  shear 
banding  and/or  subsequent  sliding  of  the  penetrator  across  the  target; 
the  surface  contains  no  tensile  cracks.  Therefore,  although  we  cannot 
entirely  eliminate  tensile  cracking  in  the  penetrator  from  the 
scenario,  we  believe  that  shear  banding  and  subsequent  separation  along 
the  banded  surface  are  the  important  mechanisms  in  penetrator  erosion. 

2.  Analysis  of  the  Target 

The  piece  of  the  recovered  4-inch-thick  RHA  target  plate  sent 
to  us  by  BRL  is  seen  in  Figure  3(a)  and  (b).  The  penetrator  entered 
at  the  upper  left  at  an  angle  of  approximately  70°,  deforming  and 
removing  significant  amounts  of  the  RHA  as  it  travelled  into  the  target, 
until  it  finally  stopped  fairly  close  to  full  perforation  with  its 
axis  rotated  to  a  more  perpendicular  angle  with  respect  to  the  target 
face.  Since  the  piece  we  received  was  cut  on  a  plane  going  through 
the  axis  of  impact  (referred  to  here  as  the  axial  plane),  it  contained 
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MB-7893-87 

FIGURE  2  DETAILED  VIEW  OF  PROXIMAL 
REGION  SHOWING  SHEAR  BANDS 
AND  ASSOCIATED  CRACKS 

Arrowi  point  to  material  completely 
surrounded  by  shear  bands. 
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FIGURE  3  TWO  VIEWS  (a)  AND  (b)  OF  THE  PIECE  OF  THE  RECOVERED 
RHA  TARGET  PLATE.  AND  THREE  SLABS  (cl  OF  THE  PLATE 
MADE  BY  SLICES  AT  LOCATIONS  SHOWN  IN  (bl 
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roughly  half  of  the  damaged  target.  Assuming  that  the  remaining  half 
of  the  damaged  target  Is  similar,  rou^ly  5  kg  of  RHA  was  eroded  from 
the  Impact  side  by  the  tl  kg  penetrator,  leaving  an  irregular  scalloped 
and  scoured  surface. 

Three  slices  perpendicular  to  the  axial  plane  were  made  in  the 
SHA  target  piece  so  as  to  view  damage  farther  away  from  the  impact  axis. 
Three  of  the  slabs  created  are  shown  in  Figure  3(c) .  The  planes  made  by 
the  slices  are  referred  to  here  as  the  transverse  planes  1,  2,  and  3. 

A  detailed  view  of  part  of  the  axial  plane  Is  shown  in  Figure 
A,  which  reveals  most  of  the  important  aspects  of  the  damage  phenomenology. 
Because  of  the  extreme  difficulty  In  polishing,  etching,  and  metallograph- 
Ing  a  32-kg  piece  of  steel,  the  surface  shown  in  Figure  A  is  polished 
but  unetched.  Hence  shear  bands  do  not  show  up  white  in  this  figure. 
However,  subsequent  photographs  of  etched  surfaces  of  the  easier-to- 
handle  sliced  slabs  (shown  in  Figure  3c)  confirmed  the  presence  of  the 
shear  banding  failures  reported  below. 

The  two  principal  features  of  target  damage,  as  seen  in  Figure 
A,  are  shear  failure  and  separation  along  planes  of  weakness  and  the 
shearing  off  and  removal  of  the  more  damaged  material.  It  had  been 
determined  (see  Section  III)  that  a  A-inch-thick  RHA  plate  has  a 
significant  anisotropy  in  its  resistance  to  shear  banding.  Planes 
parallel  to  the  rolling  plane  undergo  shear  banding  at  strains  well 
below  those  required  in  any  other  planes.  Thus  we  see,  inside  the 
target,  material  that  has  sheared  and  separated  along  the  planes  of 
weakness.  These  shear  failures  often  extend  10  cm  or  more,  and  by  studying 
Figure  3(c),  one  can  see  the  three-dimensional  nature  of  the  distorted 
sheared  planes. 

The  damage  process  is  somewhat  analogous  to  a  deck  of  cards 
that  is  held  firmly  at  the  left  end  and  bent  downward  at  the  right. 
Individual  cards  will  slide  with  respect  to  their  nearest  neighbors 
and  severe  bending  may  cause  a  gap  to  grow  between  the  cards.  The 
upper  cards  will  slide  to  the  left  with  respect  to  the  lower  cards, 
just  as  the  upper  layers  of  the  RHA  slide  toward  the  left  with  respect 
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to  the  lower  layers  (as  Is  shown  by  the  arrow  next  to  the  longest 
sheared  planes  In  Figure  A) . 

This  relative  motion  Is  opposite  to  what  one  might  first 
expect.  A  long  rod  projectile  at  a  high  angle  of  obliquity  impacting  a 
target  with  planes  of  weakness  parallel  to  the  target  face  might  be 
expected  to  cause  the  motion  shown  schematically  in  Figure  5(a),  with 
the  projectile  dragging  the  upper  layer  of  the  target  to  the  right. 
Actually,  the  upper  layers  are  pushed  in  the  direction  opposite  to 
that  of  the  projectile.  We  believe  this  to  be  caused  (as  shown  in 
Figure  S(b)  by  eroded  material  from  the  front  of  the  projectile  traveling 
backward  along  the  body  of  the  projectile,  in  a  manner  similar  to  that 
shown  for  a  rlghc^angle  impact  In  Figure  6.  The  heavily  distorted 
material  from  the  front  of  penetrator  is  forced  to  flow  backward  along 
the  still  advancing  penetrator  body,  damaging  and  then  shearing  off 
adjacent  target  material.  The  areas  circled  in  Figure  A  are  highly 
damaged  target  regions  almost  sheared  off  in  the  direction  opposite  to 
the  motion  of  the  projectile. 

As  the  projectile  penetrates  farther  into  the  target  and  as 
the  remaining,  unpenetrated  layer  of  the  target  becomes  thinner,  there 
is  less  resistance  to  deformation  in  the  direction  perpendicular  to  the 
target  face  chan  in  the  original  direction  of  the  impact.  Hence,  the 
uneroded  portion  of  the  projectile  turns  in  the  direction  of  least 
resistance,  and  further  penetration  begins  to  look  like  a  more  normal 
Impact. 

The  phenomenology  that  occurs  in  the  final  stages  leading  to 

full  penetration  can  be  best  seen  by  examining  transverse  plane  2,  which 

* 

is  shown  in  detail  in  Figure  7.  This  transverse  plane,  which  was  cut 
perpendicular  to  the  point  of  deepest  penetration  by  the  projectile  into 
the  target,  exhibits  several  sheared  and  separated  planes  of  weakness 
extending  as  much  as  5  cm  away  from  the  axis  of  penetration.  Vertical 
cracks,  perpendicular  to  the  planes  of  weakness,  have  opened  very  near 
Che  farthest  extremity  of  these  failed  planes.  The  lowermost  of  these 
cracks  has  reached  the  back  surface  and  forms  the  lateral  boundary  of  a 

*The  image  is  reversed  with  respect  to  that  in  Figure  3(c). 
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S  TWO  CONCEPTUAL  SCENARIOS  FOR  RELATIVE 
MOTION  ACROSS  SHEARED  PLANES  OF  WEAKNESS 
IN  A  TARGET  IMPACTED  BY  A  LONG-ROD 
PENETRATOR  AT  A  HIGH  ANGLE  OF  OBLIQUITY 
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MP-7893-91 

FIGURE  6  CROSS  SECTION  THROUGH  A  TARGET  IMPACTED  PERPENDICULARLY 

BY  A  LONG-ROD  PENETRATOR,  SHOWING  BACKWARD  FLOW  OF  ERODED 
PENETRATOR  MATERIAL  (Photograph  courtesy  of  BRL) 

Note:  Horizontal  hole  near  bottom  was  predrilled  for  instrumentation  access. 
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MP-7893-92A 

FIGURE  7  TRANSVERSE  PLANE  2,  SHOWING  VERTICAL  CRACK  DEVELOPMENT 
IN  FINAL  STAGES  PRIOR  TO  COMPLETE  PENETRATION 
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large  piece  chat  vaa  very  close  co  being  ejected  from  the  back  of  the 
target.  This  would  have  been  the  large  piece,  often  called  the  "plug," 
that  has  been  seen  for  oblique  Impacts  just  above  the  ballistic  limit 
(as  reported  by  Rlcchlazzl  of  BRL  at  the  February  1980  project  meeting). 
Another  of  these  vertical  cracks,  which  already  extends  halfway  across 
rhe  material  between  two  failed  planes  of  weakness,  would  have  formed 
the  next  section  of  target  material  to  be  ejected.  Thus  It  appears 
that  the  final  stages  of  penetration,  at  least  In  cases  close  co  the 
ballistic  limit  (but  possibly  also  In  cases  well  above  the  ballistic 
limit),  are  controlled  by  the  growth  of  vertical  cracks  across  already 
failed  planes  of  weakness  and  their  Intersection  with  the  back  surface 
of  the  target. 

Figures  8  and  9  contain  higher  magnification  photographs 
of  the  two  vertical  cracks  discussed  above.  The  relative  smoothness 
of  the  crack  edge  and  the  presence  of  a  white  transformed  region  In 
Figure  8  seem  tc  Indicate  that  this  lowermost  vertical  crack  is  a 
shear  band.  However,  the  rou^ness  of  the  crack  edge  and  the  absence 
of  any  transformed  region  In  Figure  9  Indicate  the  possibility  of 
tensile  fracture  In  the  second  vertical  crack. 

Both  of  these  mechanisms  are  distinct  possibilities:  shear 
banding  has  clearly  occurred  Jn  the  plane  of  weakness,  and  continued 
shear  distortions  will  eventually  trigger  banding  across  the  planes 
of  weakness.  Tensile  failure  Is  also  a  possibility.  After  the 
material  Is  separated  Into  parallel  slabs  by  shear  failure  along 
planes  of  weakness,  further  bending  of  the  slabs  may  Induce  local 
tensile  stresses  parallel  co  the  slab  surface  at  the  upper,  or  convex, 
edge  of  the  slab  (Just  as  local  compressive  stresses  will  be  induced 
at  the  lower,  or  concave,  edge  of  the  slab).  If  the  region  Is  far 
enough  away  from  the  axis  of  penetration,  the  local  tensile  stress 
may  exceed  the  more  global  compressive  stresses  by  enough  to  cause  a 
Mode  I  tensile  crack  to  open  and  propagate  downward.  Which  mechanism 
occurs  at  any  particular  vertical  crack  site  depends  on  which  failure 
criteria — shear  banding  or  tensile — are  satisfied  first. 
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MP-7893-94 


FIGURE  9 


ofr?«ef»  ^  extending  along  an  entire  length 

OF  A  SHEAR  PLANE  OF  WEAKNESS 
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3. 


Sutmnary 


The  Information  obtained  about  the  damage  phenomenology  for 
oblique  penetration  near  the  ballistic  limit  of  a  hardened  steel 
projectile  into  an  RHA  target  can  be  summarized  as  follows.  The 
projectile  bores  into  the  target,  deforming  the  RHA  downward,  activating 
shear  failure  along  planes  of  weakness  in  the  target  and  in  the  front 
of  the  projectile.  Further  penetration  causes  shear  failure  across  the 
planes  of  weakness,  and  the  backflowing  material  from  the  front  of  the 
projectile  shears  out  the  damaged  target  material. 

Further  penetration  causes  the  projectile  to  turn  more  normal 
to  the  target  face,  which  increases  the  distortion  of  the  remaining 
ui^eroded  layers  of  the  target.  Finally,  vertical  tensile  or  shear 
cracks  at  positions  away  from  the  penetration  axis  propagate  across 
the  planes  of  weakness  and  intersect  the  back  surface,  removing  a 
large  "plug"  of  material  and  completing  the  penetration  process. 

B.  Quarter-Scale  Impact  of  DU  into  RHA  Triple  Armor 

In  October  1979,  we  received  from  BRL  a  shipment  of  recovered 
fragments  from  quarter-scale  experiments  involving  the  impact  of 
DU-3/4  wtZTi  penetrators  into  RHA  triple  armor.  Although  complete 
Information  was  not  available  on  the  test  conditions  of  the  experiments, 
we  examined  two  of  the  penetrators  metallographically  to  study  the 
failure  phenomenology. 

The  first  DU  penetrator  fragment  we  examined  was  embedded  in  what 
we  believe  is  the  last  of  the  three  RHA  target  plates,  as  seen  in 
Figure  10(a).  An  electrolytic  polishing  technique  was  used  that  did 
not  affect  the  RHA  place.  The  crater  shape  suggests  a  highly  oblique 
impact. 

Figure  10(b)  is  a  more  highly  magnified  view  of  the  embedded 
penetrator  fragment.  The  wavy  lines  running  through  the  fragment  are 
prior  Y  grain  boundaries  resulting  from  the  SSCC  solutionizlng 
treatment.  Ac  this  temperature,  the  equilibrium  phase  for  a  U-3/4  wtX 

4 

Ti  alloy  is  a  solid  solution  of  Ti  in  y  (body  centered  cubic)  uranium. 
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FIGURE  10  RHA  TARGET  PLATE  IMPACTED  OBLIQUELY  BY  A  URANIUM  3/4  wt%  TITANIUM  ALLOY  PENETRATOR 


Grain  growth  in  the  solid  solution  at  this  temperature  is  fairly  rapid, 
permitting  the  formation  of  y  grains.  Quenching  from  the  solutionizlng 
treatment  transforms  the  y  (b.c.c.)  uranium  to  a  (orthorhombic)  uranium 
by  a  diffusionless  (martensitic)  transformation.^  The  location  of  the 
prior  y  grain  boundaries  can  still  be  detected,  however,  because  the 
resultant  lens-shaped  a  grains  Initiate  and  propagate  from  the  y  grain 
boundaries . 

The  flat  surface  oriented  perpendicular  to  the  penetrator  axis  is 
most  probably  Che  call  end  of  the  penetrator.  The  fragment  to  the  right 
of  this  surface  was  probably  one  of  the  many  eroded  from  the  penetrator 
during  the  perforation  of  the  first  two  plates.  This  particular  piece, 
traveling  in  the  cloud  of  particles  trailing  the  main  body  of  the 
penetrator,  struck  and  welded  Itself  to  the  call  end  of  the  penetrator 
already  embedded  in  the  third  target  plate.  The  structure  of  the  inter¬ 
face  between  Che  residual  penetrator  and  this  fragment  is  detailed  in 
Figure  11.  Figure  11(a)  shows  that  the  area  of  contact  between  the  two 
fragments  consists  of  a  wide,  featureless  zone,  which  suggests  that 
excessive  localized  heating  occurred  when  the  fragment  impacted  the 
tail  end  of  the  penetrator,  possibly  enough  heat  to  induce  melting  to 
resolldlflcaclon.  Figure  11(b)  shows  how  the  a-uranlum  grains  adjacent 
to  the  featureless  zones  are  severely  sheared,  suggesting  the  the  localized 
heating  was  caused  by  frictional  drag  between  the  fragment  and  the  tail 
end  of  the  penetrator. 

Examination  of  the  residual  penetrator  itself  revealed  little 
information  on  the  mechanisms  by  which  it  broke  up.  However,  examination 
of  Che  fragment  attached  to  the  call  end  of  the  penetrator  revealed  a 
number  of  what  appear  to  be  adiabatic  shear  bands.  These  bands  are 
detailed  in  Figure  12.  The  sheared  nature  of  the  a-uranlum  grains 
in  the  vicinity  of  the  bands-  Figure  12(b)  suggests  chat  a  large 
amount  of  localized  plastic  deformation  has  occurred.  The  featureless 
interior  of  the  bands  suggests  chat  rapid  heating  and  cooling  occurred 
during  this  shearing  process.  Aa  Figure  12(a)  indicates,  the  bands 
extend  for  large  distances  through  the  fragment.  Since  this  fragment 
was  separated  from  the  main  penetrator  early  in  the  penetration  process. 
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FIGURE  1 1  DETAILED  VIEWS  ALONG  THE  FLAT  SURFACE  ORIENTED  PERPENDICULARLY  TO  PENETRATOR  AXIS 
[As  seen  in  Figure  10(1))  1 


FIGURE  12  ADIABATIC  SHEAR  BANOS  IN  URANtUM-3/4  wl%  TITANIUM  ALLOY  PENETRATOR  FRAGMENT 
EMBEDDED  IN  STEEL  TARGET  PLATE 


Ic  appears  that,  at  lease  In  the  early  stages  of  penetration,  the 
penetrator  was  fractured  by  the  nucleaclon,  growth,  and  coalescence  of 
adiabatic  shear  bands.  The  failure  mechanisms  that  were  operative 
later  could  not  be  surmised  from  the  residual  penetrator  fragment. 

The  second  DU  fragment  examined  was  recovered  after  complete 
perforation  through  the  triple  armor.  Its  impact  velocity  was  985  m/s. 
Its  residual  velocity  was  S62  m/s,  and  Its  residual  mass  was  12.8  g. 

The  original  angle  of  obliquity  is  unknown.  A  cross-sectional  view  of 
the  penetrator  is  shown  in  Figure  13(a).  Detailed  examination  of  the 
penetrator  shaft  away  from  the  impact  end  (Figure  13b)  revealed  no 
evidence  of  any  activated  microfailure  mechanisms.  The  visible  porosity 
was  caused  by  the  pullout  of  uranium  oxide  inclusions  during  grinding 
and  electropollshing,  and  not  by  plastic  flow.  The  globular  inclusions 
are  crapped  inside  the  ingot  during  casting  and  are  then  strung  out  in 
lines  parallel  to  Che  rod  axis  during  extrusion. 

The  situation  is  different  near  the  impact  end  of  this  penetrator. 

Ic  is  evident  from  Figute  13(a)  that  extensive  localized  failure  occurred 
during  impact.  Figure  14  shows  chat  this  failure  was  caused  by  adiabatic 
shear  band  formation.  The  lens-shaped  a-uranlum  grains  are  noticeably 
sheared  in  the  vicinity  of  these  bands  (Figure  14b).  From  Figure  14(a), 
ic  is  also  apparent  chat  separation  along  the  shear  band  is  not  simply 
due  CO  brittle  crack  propagation.  The  surfaces  of  the  separated  shear 
band  are  quite  irregular,  suggesting  that  void  coalescence  contributed 
to  the  separation  process.  Irrespective  of  the  actual  separation  mech¬ 
anism,  it  is  clear  that  adiabatic  shear  band  nucleation,  growth,  and 
coalescence  played  a  dominant  role  in  Che  breakup  of  this  penetrator 
during  the  perforation  of  spaced  triple  armor. 

It  is  evident  from  the  observations  made  on  the  two  uranium-3/ 4  wt% 
titanium  penetrator  fragments  that  adiabatic  shear  banding  does  occur 
in  this  alloy  and,  furthermore,  is  the  dominant  failure  mechanism  for 
erosion  of  the  penetrator  nose  during  a  high  velocity  impact. 


28 


I 

i 

1 

I 


i 

I 


1 

i 

I 


I 


(a)  Ovarall  view  of  the  recovered  penetrator  showing 
how  the  impact  end  was  fractured  and  the  shaft 
bent  during  perforation. 


(b)  View  of  the  microstructura  of  the  penetrator 
shaft  back  from  the  impact  end.  The  holes 
most  likely  resulted  from  oxide  inclusions 
strung  out  along  the  extrusion  axis  and  then 
pulled  out  during  grinding  and  electropolishing. 


MP-7893-70 

FIGURE  13  RECOVERED  URANIUM-3/4  wt%  TITANIUM  PENETRATOR 
AFTER  PERFORATION  OF  SCALED  TRIPLE  ARMOR 
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(b)  Higher  magnification  view  of  the  tip  of  the  thear  band 
in  (a).  Shearing  of  the  a-uranium  platelets  near  the 
unseparated  part  of  the  shear  band  is  quite  visible.  Note 
the  offset  prior  7  grain  boundary  in  the  upper  left-hand 
corner  of  the  micrograph. 

MP-7893-71 

FIGURE  14  DETAILED  VIEWS  OF  THE  FRACTURED  IMPACT  END 

OF  THE  RECOVERED  PENETRATOR  SEEN  IN  FIGURE  13(a) 
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C.  Quarter-Scale  Impact  of  Steel  Into  RHA  at  Zero  Obliquity 

To  obtain  downrange  fragmentation  data,  we  fired  a  1.02-cm-dlamecer, 
blunt-nosed  A340  steel  (R^AO)  long-rod  projectile  at  a  normal  angle  and 
a  velocity  of  1.55  km/s  through  a  2. 5A-cm-thlck  plate  of  RHA  (cut  from 
a  lO-cn-thlck  slab).  This  experiment  Is  discussed  In  detail  In  Section 
IV.  In  addition,  we  examined  the  recovered  target,  residual  penetrator, 
and  some  downrange  fragments  to  add  to  our  phenomenology  data  base,  and 
In  particular  to  compare  the  failure  mechanisms  evident  In  RHA  subjected 
to  a  high  velocity  normal-angle  penetration  with  those  of  the  high  velocity 
hlgh-obllqulty  penetration  discussed  above  and  with  those  of  the  lower 
velocity,  normal-angle  tests  performed  last  year. 

Photographs  of  these  recovered  sanples  are  shown  In  Figure  15,  and 
polished  and  etched  cross  sections  of  the  residual  projectile  and  target 
place  are  shown  In  Figures  16  and  17,  respectively.  The  proximal  end 
of  the  residual  A3A0  steel  projectile  exhibits  shear  banding  failure 
In  planes  approximately  45*  to  the  Impact  axis,  very  similar  to  that 
shown  in  Figure  1  and  2  for  the  full-scale  oblique- impact  case.  This 
Implies  that  the  phenomenology  of  penetrator  erosion  at  high  Impact 
velocities  is  not  strongly  dependent  on  the  angle  of  impact. 

Figure  17  shows  a  highly  gouged  and  scalloped  inside  edge  of  the 
RHA  target  hole  and  a  very  high  density  of  shear  bands  on  planes  of 
weakness  in  this  highly  anisotropic  material.  Figure  17(c)  shows 
clearly  that  the  material  falls  along  these  bands.  This  failure  Is 
followed  by  increased  distortion  as  the  material  slides  along  the 
failed  planes  until  finally  the  highly  damage  material  is  sheared  out 
by  bands  cutting  across  the  planes  of  weakness  (as  shown  in  Figure  17b) . 

That  is  what  gives  the  surface  its  scalloped  appearance  and  is  the 

* 

same  process  shown  in  Figure  4  for  an  oblique  impact.  Furthermore, 

as  can  be  seen  in  the  bottom  right-hand  comer  of  Figure  17(a),  a 

* 

Note  that  for  the  oblique  Impact  case,  the  final  shearing  across  the 
planes  of  weakness  occurs  in  the  direction  opposite  to  that  of  the 
projectile  motion  because  of  the  backflowing  erosion.  In  the  normal 
impact  case,  the  final  shearing  is  in  the  direction  of  projectile 
motion  except  in  approximately  the  upper  quarter  of  the  target  place, 
where  projectile  backflow  is  again  dominant. 
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(b) 

MP-7893-1 12 


FIGURE  16  POLISHED  AND  ETCHED  CROSS  SECTION  (al  OF  THE  PROXIMAL 
END  OF  THE  RESIDUAL  4340  STEEL  PENETRATOR  FRAGMENT, 
WITH  DETAIL  (b)  OF  SHEAR  BANDED  REGION 
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FIGURE  17  POLISHED  AND  ETCHED  CROSS  SECTION  (a)  THROUGH  RHA  TARGET 
AND  DETAILED  VIEWS  (b|  AND  (c)  SHOWING  DENSE  SHEAR  BANDING 
IN  REGION  ADJACENT  TO  PENETRATION  HOLE 
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large  fragment  (extending  ~  120‘  of  circumferential  arc  and  weighing 
s  7  g)  was  removed  by  the  propagation  of  a  vertical  tensile  crack  from 
a  failed  plane  of  weakness  to  the  back  target  surface  in  a  location  far 
removed  from  the  penetration  axis.  Again  this  is  similar  to  what  is 
observed  for  an  oblique  impact  (see  Figure  7). 

Thus  It  appears  that  much  of  the  same  phenomenology  occurring  in 
a  high  velocity  oblique-angle  Impact  is  evident  in  a  high  velocity 
normal-angle  Inqjact.  As  a  contrast.  Figure  18  shows  the  results  of  a 
blunt-nosed  normal-angle  Impact  performed  last  year  at  about  half  the 
velocity  of  the  test  described  above  (0. 792  km/s,  or  just  above  the 
ballistic  limit),  using  the  same  RHA  target  material  and  a  somewhat 
harder  4340  steel  penetrator  (R^50).  In  the  lower  velocity  case,  no 
shear  banding  along  planes  of  weakness  is  observed,  and  the  smooth 
appearance  of  the  plug  walls  Is  Indicative  of  a  vertical  shear  band 
traversing  the  entire  thickness  of  the  target.  These  two  types  of 
penetration  behavior,  as  seen  in  Figures  17  and  18,  are  significantly 
different,  and  our  computer  code  must  be  able  to  predict  the  differences, 
using  the  same  failure  models  and  material  properties,  and  varying  only 
the  impact  velocity. 

Finally,  we  examined  metallographically  some  of  the  downrange 
fragments  that  were  recovered.  The  samples  were  chosen  to  span  the  size 
range  of  the  recovered  fragments  and  were  sliced  in  random  directions. 

The  polished  and  etched  surfaces  from  the  4340  steel  penetrator  fragments 
are  shown  in  Figure  19  and  those  from  the  RHA  target  fragments  are 
shown  in  Figure  20.  The  white-etching  bands  that  show  up  readily  on  the 
etched  4340  surfaces  occur  on  the  internal  or  not-yet-separated  surfaces, 
as  well  as  on  the  external,  already-separated  surfaces.  Shear  banding 
is  less  apparent  on  the  RHA  surfaces  becuase  the  lower  carbon  content 
decreases  the  tendency  to  form  martensite.  Nevertheless,  one  can 
easily  see  the  large  shear  deformations  in  the  grain  structure  and  the 
long  cracks  running  through  these  regions.  One  can  also  see  in  some 
of  the  RHA  fragments  (as  at  the  left  side  of  Figure  20a)  what  appears 
to  be  a  string  of  ductile  voids.  It  is  difficult  to  cell  by  metallo- 
graphlc  inspection  whether  this  is  subsequent  tensile  failure  along  an 
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(b) 

MP-7893-42A 

FIGURE  16  BLUNT-NOSED  4340  STEEL  PROJECTILE  EMBEDDED  IN  RHA  TARGET 
■•LUG  (a)  AND  CROSS  SECTION  THROUGH  IMPACT  AXIS  (b) 

FOLLOWING  IMPACT  AT  792  m/s 
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(b) 

MP-7893-119 


FIGURE  19  S  TIONS  THROUGH  RECOVERED  DOWNRANGE  4340  STEEL 
P6NETRATOR  FRAGMENTS,  SHOWING  EXTENSIVE  INTERNAL 
AND  EXTERNAL  SHEAR  BANDING 
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(b) 

MP-7893-120 

FIGURE  20  SECTIONS  THROUGH  RECOVERED  OOWNRANGE 
RHA  TARGET  FRAGMENTS.  SHOWING  EXTENSIVE 
SHEAR  DEFORMATION  AND  RESULTANT  BANDING 
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existing  shear  band  or  a  subsequent  shear  band  along  existing  tensile 
voids.  However,  since  our  computations  show  chat  large  shear  deformations 
precede  any  large  tensile  stresses  in  the  target  (except  perhaps  in  a 
narrow  region  adjacent  to  the  downrange  surface),  we  conclude  that  shear 
band  failure  is  the  prime  mechanism  of  importance  in  target  fragmentation. 

D.  Phenomenological  Conclusions  and  their  Significance  for  the 

Computational  Model 

The  following  five  conclusions  can  be  drawn  concerning  the  phenomenol¬ 
ogy  of  high  velocity  penetration  of  moderate-to-thick  armor  by  long-rod 
penecrators  on  the  basis  of  the  experiments  described  above: 

(1)  Shear  banding  is  the  principal  failure  mechanism  in  both 
the  penecracor  and  the  target,  and  it  is  the  coalescence  of 
shear  bands  that  forms  the  bulk  of  the  downrange  (and  uprange) 
fragments. 

(2)  Shear-band-induced  penecracor  backflow  and  shear  failure  along 
planes  of  weakness  are  significant  mechanisms  of  target  erosion. 

(3)  Projectiles  change  direction  within  the  target  (oblique  impacts 
only),  seeking  the  path  of  least  resistance. 

(4)  Tensile  failure  may  be  a  factor,  but  only  in  specific  instances 
such  as  at  the  rear  of  the  target  plate  in  the  final  stages  of 
penetration  for  a  projectile  just  above  the  ballistic  limit. 

(5)  Most  of  the  failure  mechanisms  occurring  in  a  highly  oblique 
impact  will  also  be  found  in  a  normal  Impact  of  the  same 
velocity  using  the  same  materials. 

Although  we  have  examined  only  a  few  target  and  penecracor  materials 
(RHA,  DU,  4340  steel  in  this  program  and  a  tungsten  alloy  in  a  previous 
program)^  in  a  velocity  range  up  to  1.5  km/s,  we  believe  chat  these 
phenomenological  results  are  applicable  to  a  more  general  range  of 
materials  and  conditions.  For  example,  we  expect  to  see  the  same 
phenomena  at  velocities  up  to  the  start  of  the  hypervelccicy  regime,  where 
strength  effects  are  no  longer  important.  Also,  because  of  the  high 
compressive  stresses  and  shear  strains  found  in  the  critical  regions  of 
long-rod  penetration,  we  expect  the  shear  failure  scenarios  to  be 
relevant  to  most  any  material,  particularly  metals,  that  may  be  considered 
for  armor  or  penecracor  use. 
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The  purpose  of  the  phenomenology  experiments,  in  addition  to 
providing  a  deeper  understanding  of  the  physical  processes  involved  in 
an  amor  penetration,  is  to  guide  us  in  selecting  the  types  of  failure 
models  and  other  features  to  be  put  in  our  predictive  computer  code  and 
the  types  of  materials  characterization  experiments  to  be  performed  on 
the  materials  of  interest.  From  each  of  the  five  conclusions  listed  above, 
we  can  elicit  a  guideline  to  aid  us  in  our  future  work  on  the  computational 
modeling  of  armor  penetration.  They  are  as  follows: 

(1)  Shear  banding  parameters  are  important  for  each  material 

of  interest,  including  both  penetrator  and  target  materials. 

(2)  An  anisotropic  shear  banding  model  is  needed,  and  the  computer 
code  needs  to  be  able  to  handle  the  extreme  deformation  and 
material  ejection  that  accompany  penetrator  backflow. 

(3)  A  self-directing  or  wandering  slide  line  is  needed  to  handle 
the  rotation  of  a  projectile  in  an  oblique  impact. 

(4)  Tensile  failure  parameters  may  be  relatively  unimportant, 
particularly  in  the  projectile. 

(5)  Simulations  of  normal  impacts  (which  of  course  are  much 
easier  to  perform  than  oblique  impacts  because  of  the  two- 
dimensionality  of  the  problem)  are  valid  tests  of  the  computer 
model's  ability  to  calculate  fragmentation  distributions 
because  they  Involve  most  of  the  phenomenology  (except  for 
projectile  rotation)  occurring  in  an  oblique  impact. 


40 


Ill  MATERIALS  CHARACTERIZATION  EXPERIMENTS 


Duiing  Che  second  year  of  this  program,  we  completed  a  series  of 
nine  contained  fragmenting  cylinder  (CFC)  experiments  using  RHA  specimen 
cubes.  A  schematic  diagram  of  the  CFC  experiments  is  shown  in  Figure  21. 
Figure  22  shows  the  shear  band  orientation  modes  possible.  A  detailed 
description  of  the  CFC  method  has  been  given  earlier  by  Erlich  ec  al.^ 

Table  1  shows  the  experimental  parameters  and  the  qualitative  damage 
results  for  these  experiments. 

Since  quasl-scaclc  tensile  tests  had  revealed  a  significant 
anisotropy  in  tensile  failure  in  RHA, ^  we  expected  a  similar  anisotropy 
to  exist  with  respect  to  resistance  to  shear  banding.  Early  experiments 
with  RHA  specimen  cubes  whose  axes  were  in  the  long  transverse  direction 
of  the  RHA  slab  (the  L-specimen,  as  shown  in  Figure  23)  confirmed  our 
expectations.  The  preponderance  of  shear  bands  in  planes  nearly  parallel 
CO  the  rolling  plane  In  specimens  L-2  and  L-4  (as  seen  in  Figures  24  and 
25,  for  example)  clearly  showed  that  shear  banding  will  occur  at  lower 
strains  in  the  rolling  plane  chan  in  planes  perpendicular  to  the 
rolling  plane. 

Once  it  was  clear  that  by  using  L-specimens  we  could  not  obtain 

many  shear  bands  in  planes  perpendicular  to  the  rolling  plane,  we  switched 

to  T-specimens  (whose  axes  were  in  the  chrough-the-thlckness  direction) 

for  which  all  the  Mode  1  shear  bands  would  lie  in  the  desired  planes. 

The  first  four  T-specimen  experiments  (T-1,  T-2,  T-3,  and  T-5),  using 

PETN  explosive  densities  similar  to  chose  used  in  all  previous  CFC 
6.7 

experiments,  yielded  only  incipient  Mode  1  shear  banding.  However, 
the  filial  three  experiments  (T-7,  T-9,  and  T-10) ,  using  much  higher 
explosive  densities,  yielded  significant  quantities  of  both  Incipient 
and  large  shear  bands,  as  well  as  some  fragmentation.  Photographs  of 
these  specimens  are  shown  in  Figures  26  and  27. 

Data  analysis  for  the  RHA  CFC  experiments  is  proceeding.  The  analysis 
consists  of  three  items:  (1)  determining  the  final  wall  thickness 
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Pressed  PETN  a-- -  Initiation  Point 
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ORIENTATION  MODE  1 


MODE  2 


B  —  Length  in  Direction 
Perpendicular  to 
Slip  Motion 

d  —  Depth  Along  45* 
Slip  Plane 

8  —  Shear  Oiiplacement 
Along  Slip  Plane 

Z  —  Axial  Direction 

r  —  Radial  Direction 

0  —  Circumferential 
Direction 


FIGURE  22  GEOMETRY  AND  NOMENCLATURE  FOR  SHEAR  BANDS 

IN  fragmenting  cylinder  experiments 
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Table  1 


PARAMETERS  AND  QUALITATIVE  RESULTS  FOR  RHA 
CONTAINED  FRAGMENTING  CYLINDER  EXPERIMENTS 


Specimen 

Luclte 

Buffer 

Thickness 

(cm) 

PETN 

Density 

(R/cm’) 

Calculated 

Detonation 

Pressure 

(kbar) 

Qualitative  Shear 

Banding  Results 

* 

L-4 

0.873 

1.2 

122 

Substantial  Mode  1  banding  con¬ 
centrated  in  rolling  direction. 

L-2 

0.873 

1.0 

80 

Same  as  L-4 . 

T-l'*’ 

0.873 

1.2 

122 

A  few  incipient  Mode  1  bands 
(some  tensile  cracking  on  outer 
surface  in  midaxial  region) . 

T-3 

0.873 

1.35 

162 

Incipient  Mode  1  banding  over 
most  of  inner  surface  (some 
tensile  cracking  as  in  T-1  above) • 

T-2 

1.600 

1.2 

122 

Few  incipient  Mode  1  bands. 

Moderate  Mode  2  banding  in  mid- 
axial  region. 

T-5 

1.600 

1.35 

162 

Extensive  incipient  Mode  1  banding 
Few  moderate  Mode  2  bands  in  mid- 
axial  region. 

T-7 

1.600 

1.47  to 
=  1.77 

202  to 
=  324 

Substantial  Mode  1  banding. 
Substantial  Mode  2  banding  in 
narrow  midaxial  region. 

T-9 

1.27 

1.54 

228 

Substantial  Mode  1  banding. 

Mode  2  bands  in  midaxial  region 
cuts  specimen  in  two. 

T-10 

1.27 

1.6 

250 

Same  as  T-9,  with  some 
fragmentation. 

L-specimens  have  axes  parallel  to  long  transverse  direction  of  RHA  block. 

T-specimens  have  axes  parallel  to  short  transverse  (through  the  thickness) 
direction. 
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FIGURE  23  CFC  SPECIMENS  MACHINED  FROM  RHA  SLAB.  SHOWING  MODE  1  SHEAR 
BANDS  IN  PLANES  PERPENDICULAR  TO  ROLLING  PLANE 


MP-7893-43 


FIGURE  24 


TWO  VIEWS  OF  HALF  OF  THE  RHA  SPECIMEN  TUBE  L-4  RECOVERED  FROM  CFG 
EXPERIMENTS 

Detonation  proceeds  from  top  to  bottom. 
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WP  7893-44 

FIGURE  25  polished  AND  ETCHED  SECTION  THROUGH  RHA  SPECIMEN  TUBE  L-4,  SHOWING 
PREPONDERANCE  OF  SHEAR  BANDING  DAMAGE  IN  ROLLING  DIRECTION 
Plane  is  cut  at  a  distance  of  2-3/8  inches  from  bottom  of  specimen  tube. 
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(and,  for  lower  damage  cases,  the  final  inside  radius)  of  che  recovered 
specimen  tube  as  a  function  of  axial  position,  (2)  determining  che  shear 
band  size  distribution  as  a  function  of  axial  position,  and  (3)  determining 
the  shape  or  outline  of  the  shear  band  in  its  slip  plane. 

Final  wall  thickness  measurements  are  made  after  slicing  the  specimen 
cube  parallel  to  the  axis.  They  are  first  used  for  comparison  with 
initial  no-damage  computer  simulations  of  the  CFC  experiments  to  check 
if  the  explosive  parameters  and  bulk  material  properties  (before  shear 
band  damage)  used  in  the  code  are  correct.  Later  che  measurements  are 
used  with  damage  simulations  to  determine  if  further  plastic  strain 
following  shear  band  nucleation  is  apportioned  correctly  between  bulk 
strain  and  localized  (shear  banding)  strain. 

The  final  wall  thickness  plots  for  all  nine  RHA  CFC  experiments 
are  shown  in  Figure  28,  along  with  qualitative  damage  values.  Note 
chat  che  qualitative  damage  results  shown  in  Figure  28  refer  only  to 
damage  by  Mode  1  banding.  The  presence  of  Mode  2  banding,  which  accounts 
for  much  of  che  equivalent  plastic  strain  in  che  peak  of  the  curve  T-1, 
for  example,  is  not  shown. 

For  determining  che  shear  band  size  distribution  as  a  function  of 
axial  position,  che  specimen  is  cut  axially  into  segments  no  greater 
than  about  60°  of  circumferential  arc.  (For  cases  of  higher  damage,  the 
existence  of  fragments  of  che  appropriate  size  makes  cutting  unnecessary.) 
Fach  segment  Is  first  divided  into  zones  by  scribing  lines  perpendicular 
CO  che  axis,  and  then  photographed.  Light  is  provided  at  a  grazing 
angle  from  alternate  sides  to  highlight  the  Mode  1  shear  bands  intersecting 
che  surface  at  the  two  different  45°  angles.  Figure  29  is  such  a  photo¬ 
graph  of  a  segment  from  specimen  T-7.  The  zones  are  =1/4  Inch  (6.35  ram) 
wide.  Then  for  each  zone,  a  shear  band  size  distribution  is  determined, 
such  as  che  one  shown  in  Figure  30  for  a  4340  steel  specimen  from  a 
previous  program.  This  type  of  plot  will  provide  the  data  base  for 
comparison  with  computer  simulations  performed  to  obtain  shear  band 
parameters  for  RHA. 

We  have  been  experiencing  difficulty  obtaining  che  shear  band  size 
distributions  for  RHA  from  the  optical  photographs  because  of  the 
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THICKNESS  —  inches 


MA-7893-B0A 


FIGURE  28  FINAL  WAUL  THICKNESS  AS  A  FUNCTION  OF  AXIAL  POSITION  FOR  ALL 
OF  THE  RHA  SPECIMEN  TUBES  RECOVERED  FROM  CFC  EXPERIMENTS 


31 


FIGURE  29  PHOTOGRAPHS  OF  INSIDE  SURFACE  OF  ONE  PIECE  OF  SPECIMEN  7-7 
Two  views  highlight  mode  I  shear  bands  intersecting  surface  at  different  45“  a 
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myriad  of  smaller  bands  and  r.he  poor  contrast  with  the  background 
metallic  surface.  We  are  currently  investigating  scanning  electron 
microscopy  as  well  as  ocher  optical  photographic  techniques  in  hopes 
of  improving  the  situation. 

For  determining  the  shape  or  outline  of  the  shear  band  in  its 
slip  plane,  the  specimen  is  sliced  perpendicular  to  the  axis,  polished, 
etched,  and  phocomicrographed  to  obtain  an  edge-on  view  of  the  shear 
bands,  as  in  Figure  31.  The  depth  (d)  and  displacement  (B)  of  each  band 
(see  Figure  22)  are  measured.  We  then  slice  or  grind  a  known  small 
amount  off  the  specimen,  repolisli,  etch,  and  measure  d  and  B  at  this  new 
poslclcn.  Finally,  we  relate  these  measured  values  of  d  and  B  with  the 
value  of  the  length  (i)  for  the  same  band  as  measured  on  the  inside 
surface  before  slicing  the  specimen.  Figure  32  shows  results  of  a 
similar  analysis  of  data  on  4340  steel,  which  exhibits  nearly  semi¬ 
circular  shear  bands. 

Again,  we  are  experiencing  difficulty  with  the  measurement  of  d 
because  of  the  high  density  of  bands  and  because  the  poor  contrast 
obscures  the  lower  end  of  the  bands  (see  Figure  31).  One  encouraging 
note,  however,  is  chat  all  the  bands  appear  to  intersect  the  inside 
surface  of  the  specimen  tube,  as  was  the  case  for  4340  steel. ^  There¬ 
fore,  measuring  a  surface  shear  band  size  distribution  will  yield  a 
complete  picture  of  the  damage.  We  need  to  make  only  a  few  slices 
through  the  tube  perpendicular  to  the  axis  to  determine  the  relationship 
between  d,  B,  and  £. 


54 


w 

CJ 

2 

< 

oc 

oc 

< 

UJ 

I 

(O 

u. 

o 

>• 


o 

1 

o 

X 

o 

2 

i 

o 

I 


2 

UJ 

s 

O  X 
UJ  H- 
cl  Cl 

</}  UJ 

<  O 

i  Q 

z 

u.  < 
O  CD 
w  o 

X  ~ 

< 
o 

h- 

ac 

< 

-j 
D 
O 


2 

S 

X 
UJ 
I- 

UJ 

n 

^  2 

§5 

CL  -J 
OC  ^ 

UJ  U 

,T 


2  5 
< 


CO 

UJ 

oc 

D 

(J 


55 


cm 


IV  DOWNRANGE  FRAGMENTATION  DATA  BASE  ACQUISITION 


J.  Barb  and  A-  J.  Rlchiazzl  of  BRL  have  provided  us  wich  a  substantial 
amount  of  fragmentation  data  on  Impact  of  DU  into  RHA  at  various  velocities 
and  obliquities.  We  expect  these  data  to  be  of  great  value  to  use  for 
comparison  with  our  computer  calculations.  However,  we  will  have  shear 
banding  model  parameters  for  RHA  and  4340  steel  before  we  expect  to  have 
them  for  DU,  and  our  code  will  handle  i  normal  Impact  before  it  is 
expected  to  be  able  to  handle  an  oblique  impact.  Therefore,  to  make  an 
early  comparison  between  downrange  fragmentation  data  and  our  computational 
model  predictions,  and  hence  determine  at  an  early  stage  the  ability 
of  our  phenomenologically  based  model  to  predict  fragmentation,  we 
performed,  at  the  quarter-scale  facility  at  BRL,  the  penetration  experiment 
described  below. 

A  l,C2-cm-diameter  by  10. 16-cm-long  right  circular  cylinder  of  4340 
steel  (R^40)  was  launched  at  zero  obliquity  and  a  velocity  of  1.55  km/s 
into  a  2. 54-cm-thlck  by  15.24-cm-diamecer  RHA  plate  (cut  from  a  10-cm- 
thick  RHA  slab).  The  experimental  arrangement  is  shown  in  Figure  33. 

The  target  plate  was  positioned  on  the  uprange  side  of  a  large  steel 
holder  plate  having  a  7 . 62-cm-diameter  hole  concentric  with  the  impact 
axis.  Celotex  recovery  packets  were  placed  both  uprange  and  downrange 
of  the  target  to  capture  fragments  emerging  from  the  two  sides  of  the 
target  plate.  Eight  x-rays  (three  vertical  and  five  horizontal)  were 
triggered  from  appropriately  placed  screens  to  record  preimpact 
projectile  velocity  and  yaw  as  sell  as  postimpact  fragmentation. 

The  first  two  pairs  of  orthogonal  radiographs  are  shown  in  Figure  34. 
Measurements  from  these  x-rays  reveal  an  impact  velocity  of  1.55  km/s 
and  a  final  yaw  of  less  than  2®.  (The  yaw  was  fortuitously  improved  by 
the  Impact  of  the  penetrator  with  the  pusher  plate  deflector;  however, 
this  unplanned  Impact  prevented  the  deflector  from  accomplishing  its 
original  purpose.) 
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Vertical  X-ray  Tubes 
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FIGURE  33  SCHEMATIC  OF  BRL  BALLISTICS  EXPERIMENT 


FIGURE  34  VERTICAL  (V,  AND  V^t  AND  HORIZONTAL  (H,  AND  RADIOGRAPHS 
SHOWING  PENETRATOR  TRAVEL  (from  right  to  left)  BEFORE  IMPACT 

Images  shown  are  100.6  ps  apart. 


The  three  postlmpect  horizontal  radiographs  are  shown  in  Figure  35. 
The  Isage  labeled  taken  t30  us  after  Impact,  shows  the  projectile 
more  than  halfway  into  the  front  of  the  target  and  many  fragments,  most 
of  them  probably  from  the  penetrator,  backscattered  from  the  target. 

(We  were  unable  to  recover  any  of  the  backscattered  fragments  because 
the  uprange  Celotex  recovery  pocket  was  destroyed  by  debris  from  the 
pusher  plate  deflector.) 

The  two  Images  labeled  and  taken  slA5  and  285  us,  respectively, 
after  impact,  show  the  fragments  and  residual  penetrator  traveling  dowr- 
range  at  a  velocity  of  1.257  km/a  or  slightly  more  than  80Z  of  the 
impact  velocity. 

The  white  numbers  adjacent  to  each  fragment  image  identify  the 
fragment  for  computational  analysis,  and  the  white  lines  show  possible 
fragment  trajectories  between  the  two  images.  Note  that  the  large  frag¬ 
ment  with  an  "X"  on  it  is  a  pice  of  the  pusher  plate  and  not  a  fragment 
of  interest. 

The  postimpact  vertical  radiograph  V^,  taken  at  the  same  time  as 
Hj,  is  shown  in  Figure  36,  along  with  an  entirely  different  set  of  identi¬ 
fying  numbers.  The  computational  analysis  consists  of  several  steps. 
First,  the  positional  coordinates  of  the  endpoints  of  each  fragement  in 
the  two  orthogonal  radiographs  are  determined.  Computer-drawn  pictures 
based  on  these  readings  are  shown  in  Figure  37.  In  each  radiograph, 
about  350  fragments  are  identified,  and  each  fragment  is  represented 
by  a  triangle  that  gives  an  accurate  length  in  the  direction  of  penetra¬ 
tion  and  an  approximate  width  parameter.  A  computer  program  then  matches 
the  fragments  images  from  the  two  radiographs  and  calculates  the  true 
length,  velocity,  and  trajectory  angles  for  each  fragment. 

The  next  step  in  the  analysis  is  to  determine  the  mass  distribution 
of  these  fragements.  The  total  mass  available  for  the  fragments  is 
roughly  98  grams:  the  3A  grams  lost  from  the  target  plus  the  39-gram 
residual  penetrator  plus  the  r25  grama  eroded  from  the  penetrator. 

However,  this  total  is  reduced  by  the  unknown  mass  of  the  fragments 
backscattered  uprange  (mostly  from  the  penetrator).  Our  method  was  to 
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FIGURE  35  HORIZONTAL  RADIOGRAPHS 

REVEALING  BACKSCATTER  FRAG¬ 
MENTATION  (H3  at  «  30  fii  after 
impact)  AND  DOWNRANGE  FRAG¬ 
MENTATION  (H4  and  H5  at  =«  14B 
and  285  p$,  respectively,  after  impact). 

White  numbers  and  lines  identify  frag¬ 
ments  and  their  paths  for  computa¬ 
tional  analysis. 


locate  the  downrange  fragments  captured  in  the  recovery  packet,  determine 
the  trajectory  of  each  fragment  from  the  position  coordinate  of  the  entrance 
hole  on  the  aluminum  witness  plate,  and  measure  each  fragment's  mass 
and  length.  By  comparing  trajectories,  lengths,  and  general  shapes 
from  the  witness  plate  data  with  those  from  the  x-ray  data,  we  might  be 
able  to  make  a  one-to-one  correlation  betweenmany  of  the  fragments  (mostly 
the  larger  ones)  recovered  and  their  images  on  the  radiographs.  In 
addition,  a  mass-versus-length  curve  can  be  obtained  and  used  to  estimate 
the  mass  for  those  radiographic  images  that  could  not  be  directly 
correlated  with  the  recovered  fragments. 

The  hole  pattern  on  the  0.5-mm-thick  aluminum  witness  plate  is 
shown  in  Figure  38  along  with  a  detailed  photograph  of  the  area  around 
the  axis  of  Impact.  The  triangles  in  Figure  38(b)  represent  dents  made 
by  penetrator  or  target  fragments  (or  perhaps  in  some  cases  other  debris) 
with  insufficient  energy  to  either  stick  to  or  perforate  the  witness 
plate.  The  300  circles  represent  the  fragments  that  did  stick  to  or 
perforate  the  plate;  about  240  of  these  fragments  were  recovered,  ranging 
in  length  from  less  than  0.6  mm  to  16.7  mm  and  in  mass  from  less  chan 
0.2  mg  Co  more  Chan  2  g. 

A  relationship  between  the  maximum  dimension  (2)  of  any  fragment  and 
the  fragment  mass  (m)  was  obtained  by  plotting  the  data  and  by  fitting 
them  to  a  fifth  order  least-squares  polynomial  of  the  form 

5  A,  . 

£n  m  ■  ^ 
i»o 

where  m  ■  mass  in  g,  f,  ■  length  in  mm,  and  through  A^  »  -9.62,  3.54, 
-0.790,  9.85  X  10  -5.89  x  10  and  1.32  x  10  respectively.  The 
mass  versus  length  data  and  the  polynomial  fit  are  shown  in  Figure  39. 

A  preliminary  calculation  of  the  downrange  fragment  mass  and  velocity 
distributions  was  performed  at  BRL,  using  the  above  mass-versus-length 
curve,  but  without  attempting  to  correlate  individual  fragments.  The 
results  are  shown  in  Figure  40.  These  data  were  transformed  into  the 
cummulative  mass  distribution  shown  in  Figure  S2  of  the  Summary  (and 
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FIGURE  36  VERTICAL  RADIOGRAPH  V5  REVEALING  DOWNRANGE  FRAGMENTATION 
AT  *  285 /is  AFTER  IMPACT. 

White  numbers  identify  fragments  for  computational  analysis. 


(b)  SIDE  VIEW 


FIGURE  37 


MP-7893-110 

COMPUTER  PLOTS  OF  FRAGMENTS  OBSERVED  ON  (•)  VERTICAL  RADIOGRAPH  Vj 
AND  (b)  HORIZONTAL  RADIOGRAPH  Hg.  AS  RECEIVED  FROM  BRL 


MP-7893-113 

FIGURE  38  OVERALL  VIEW  (a)  OF  WITNESS  PLATES  WITH  DETAIL  (b) 
OF  HEAVILY  IMPACTED  REGION  WITH  FRAGMENT  HOLES 
(Circles)  NUMBERED, 

Triangles  denote  impacts  that  did  not  perforate  aluminum. 


MA-7893-n7 


FIGURE  39  MEASURED  MASS-VERSUS-LENGTH  DATA  FOR  THE  RECOVERED 
FRAGMENTS  AND  A  LEAST-SQUARES  FIFTH-ORDER  POLYNOMIAL 
FIT  TO  THE  DATA 
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MA-7883-n8 


FIGURE  40  PRELIMINARY  MASS  (a)  AND  VELOCITY  (b|  DISTRIBUTIONS 
CALCULATED  FOR  QUARTER-SCALE  TEST.  AS  RECEIVED 
FROM  BRL 
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discussed  In  Section  VII)  to  compare  with  computational  predictions. 
Continued  analysis  of  this  experiment,  including  a  one-to-one  correlation 
between  the  larger  fragments  and  their  radiographic  Images,  is  proceeding. 
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V  SHEAR  BAND  MODEL  DEVELOPMENT 


The  shear  band  model  has  been  developed  co  represent  Che  primary 
fracture  process  of  ductile  materials  under  compression  and  large 
shear  strains.  The  model  represents  behavior  In  penetrator  and  armor 
during  long-rod  impacts  and  In  exploding  cylinders.  The  current  model, 
SHEAR3,  Is  an  extension  of  earlier  models  developed  for  fragmenting 
cylinders.^  ^  SHEAR3  la  based  on  a  careful  study  of  a  series  of  rontalned 
fragmenting  cylinders  of  43A0  steel.  The  model  has  now  been  fitted  to 
these  data  and  appears  Co  fit  all  aspects  of  the  data  fairly  well.  From 
this  agreement,  we  have  developed  a  general  procedure  for  fitting  the 
model  to  data  from  fragmenting  cylinders.  Three  main  topics  are 
covered  in  this  section:  results  of  fitting  SHEAR3  co  the  A340  data, 
modification  of  SHEAR3,  and  the  fitting  procedure. 

A.  Shear  Band  Parameters  for  4340  Steel,  R  40 

c 

The  4340  data  were  obtained  from  eight  contalned-fragmenting- 
cylinder  experiments.^  Shear  band  parameters  (listed  in  Table  2)  for 
the  SHEAR3  model  were  obtained  to  represent  these  data.  The  nature  of 
the  shear  band  data  and  a  comparison  of  the  simulated  results  with  these 
data  are  described  here. 

The  shear  band  data  were  obtained  by  counting  and  measuring  the 
shear  bands  observed  on  the  surfaces  of  six  cylinders  after  contained- 
fragmentlng-round  experiments.  (Negligible  numbers  of  bands  were 
observed  in  the  other  two  cylinders.) 

As  shown  in  Figure  21  In  Section  III,  the  experimental  configuration 
consists  of  the  test  cylinder  filled  with  explosive  and  surrounded  by 
three  cylinders:  an  acrylic  buffer  cylinder,  a  steel  confining  cylinder, 
and  an  outer  lead  cylinder.  The  buffer  thickness  and  Che  density  of 
Che  explosive  were  varied  to  provide  a  range  of  strain  races,  peak  strains 
and  levels  of  shear  banding.  The  explosive  at  either  end  of  the  cylinder  is 
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Table  2 


INPUT 

_ Symbol _ 

General 

K 

G 

r 

I 

Shear  band  array 

RFR(M,22) 

BFR(M,23) 

BFR(M,24) 

BFR(M,25) 

BFR(M.26) 

BFR(M,27) 

BFR(M,28) 

BFR(M,29) 

BFR(M,30) 

BFR(M,31) 

BFR(M,32) 

BFR(M,33) 

BFR(M,34) 

Yield  strength  array 
YS 

ES 


DATA  FOR  SHEAR  BANDING  IN  4340  STEEL,  R^40 


Value 


Definition 


7.85  g/cm^ 

1.59  X  10^^  dyn/cm^ 
8.19  X  lOl^  dyn/cm" 
1.69 


Density 

Bulk  modulus 
Shear  modulus 
Grilnelsen  ratio 


14.0 
35.0 
0.01  cm 

10  per  cm 
0.15 
0.07 
0.04  cm 


0.25 

0.577 

2 

2- 


10 


Growth  coefficient,  C 
Nucleation  coefficient,  A^ 
Nucleation  size  parameter,  R^ 

Number  at  first  nucleation 
Threshold  plastic  shear  strain 
Shear  strain  ratio, b 
Maximum  nucleation  band 
radius,  R^j^ 

Volume  effect  coefficient,  c 
Tan  <5.  friction  coefficient 
on  shear-banded  planes 
Storage  Indicator.  "2"  means 
store  all  plastic  strains 
Definition  of  t.  "2"  means 
T  -  TAUZ  (2) 

Work-hardening  modulus, 
unused  (see  yield  parameters 
below) 

Number  of  radii  in  each 
orientation 


1.03  X  10^*^  dyn/cm? 
1.07  X  10^®  dyn/cm- 
1.37  X  10^®  dyn/cm 


Yield  strength  at  three 
points  on  the  nonlinear 
work-hardening  curve 


0.,0.03,  1.0  Plastic  strains  at  three 

points  on  the  nonlinear 
work-hardening  curve. 


72 


freer  co  expand  chan  chac  In  che  center  so  the  central  portion  of  the  test 
cylinder  experiences  a  longer  duration  of  loading  and  higher  strains. 

As  a  consequence,  che  greatest  damage  occurs  In  the  center  of  che  cylinder, 
there  is  a  reduction  of  damage  with  distance  from  the  center,  and  no  damage 
occurs  at  the  ends  where  there  was  no  explosive  In  direct  contact. 

Following  counting  of  the  shear  bands,  the  cumulative  numbers  were 
plotted.  Then  each  distribution  was  fitted  to  che  exponential  relation 

n  ■  n  exp  (-r/r, )  (1) 

go  1 

2 

where  n  is  che  cumulative  number  of  bands/cm  with  a  radius  greater 
8 

chan  r 

2 

n  is  che  total  number  of  bands/cm 
o 

r^  is  a  size  parameter  for  the  distribution. 

With  the  distributions  reduced  to  two  characterizing  parameters,  n^  and 
rj^,  plots  were  made  to  study  the  initiation,  nucleation  rate,  and  growth 
rate  of  che  bands.  These  two  characterizing  parameters  are  exhibited  in 
Figures  41  and  42  as  functions  of  the  axial  position  along  che  cylinder. 
(Detonation  began  at  the  left  end.) 

To  match  these  data  with  a  shear  band  model,  we  modified  SHEAR3 
slightly  (as  described  later)  and  conducted  simulations  of  the  fragmenting 
cylinder  experiments  with  estimated  sets  of  the  shear  hand  parameters. 

The  simulation  results  were  compared  with  the  experimental  data,  revised 
values  of  che  parameters  were  estimated,  and  the  simulations  were  repeated 
until  reasonably  good  agreement  with  che  experimental  data  was  obtained. 
(This  multiple  simulation  process  1$  described  later.)  Most  of  these 
simulations  were  made  for  tests  2  and  8,  which  bracket  the  range  of  damage 
observed.  Figures  43  and  44  compare  che  measured  and  simulated  charac¬ 
terizing  parameters  n^  and  respectively,  for  tests  2  and  8. 

The  simulated  results  tend  to  be  within  the  range  of  scatter  In  the 
experimental  data.  Also  the  simulations  show  the  flat  top  (on  the  curves) 
Chat  is  characteristic  of  the  test  data.  The  simulations  show  some  Mode  2 
and  3  bands,  but  there  are.  many  more  Mode  1  bands;  this  also  agrees  with 
the  data. 
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PARAMETER  r,  OF  THE  MEASURED  SURFACE  DISTRIBUTIONS 


AXIAL  POSITION  ALONG  THE  CYLINDER  —  inches 

MA>7893-B3 

FIGURE  41  VARIATION  OF  THE  DISTRIBUTION  PARAMETER  r,  WITH  AXIAL 
POSITION  FOR  SIX  CONTAINEO-FRAGMENTING-CYLINDER 
EXPERIMENTS  IN  4340  STEEL  AT  R  ■  40 
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TOTAL  NUMBER  OF  SHEAR  BANDS  FROM  THE  SURFACE  DISTRIBUTIONS 


TOTAL  NUMBER  OF  SHEAR  BANDS  PER  UNIT  AREA  —  number/cm^ 


MA-7B93'9S 


FIGURE  43  COMPARISON  OF  THE  VARIATION  OF  THE  NUMBER  OF  SHEAR  BANDS 
ng  WITH  AXIAL  POSITION  FOR  CONTAINED-FRAGMENTING-CYLINOER 
EXPERIMENTS  IN  4340  STEEL 
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PARAMETER  r,  OF  THE  MEASURED  AND  COMPUTED  SURFACE  DISTRIBUTIONS  —  cm 


AXIAL  POSITION  ALONG  THE  CYLINDER  —  cm 

MA-78g3-96 

FIGURE  44  COMPARISON  OF  THE  VARIATION  OF  THE  SIZE  PARAMETER  r, 

WITH  AXIAL  POSITION  FOR  CONTAINED-FRAGMENTING-CYLINDER 
EXPERIMENTS  IN  4340  STEEL 


77 


B.  Modification  of  the  Shear  Band  Model  SHEAR3 

While  fitting  Che  shear  band  model  to  Che  4340  data,  we  have  been 
led  to  modify  the  nucleaclon  and  growth  processes  and  the  method  for 
representing  Che  shear  band  size  distribution.  Also,  for  listing  pur¬ 
poses,  the  subroutine  now  transforms  the  shear  band  distribution  to  a 
surface  distribution  and  fits  the  latter  distribution  to  an  exponential 
for  ease  In  comparing  with  experimental  data.  A  nonlinear  work-hardening 
process  has  been  Included  in  SHEAR3.  The  interaction  of  the  uucleation 
and  growth  processes  with  stress  relaxation  and  the  assignment  of  shear 
band  damage  to  planes  are  both  discussed. 

1.  Nucleatlon 

In  Che  current  model,  nucleatlon  begins  separately  on  each  of 
the  selected  planes  when  the  plastic  shear  strain  on  chat  plane  exceeds 
a  critical  value  e^.  During  nucleatlon,  bands  are  initiated  In  a  range 
of  sizes  with  a  distribution  given  by 


"g  ■  “♦e 


<2) 


where  Is  the  number  of  bands/cm  with  radii  greater  than  R 

3 

is  Che  total  number/cm  nucleated  on  the  i^6  plane 

Is  the  nucleatlon  size  parameter,  an  input  quantity  (cm). 

The  number  nucleated  Is  given  by  Che  following  relation 


where 


>9 


Vo 


N^q,  ar®  the  current  and  previous  numbers  of  bands/cm 

In  the  <1)9  orientation 


(3) 

3 


Ac 


<t)9 


Is  a  nucleatlon  coefficient 

is  the  plastic  shear  strain  (on  the  (j)9  plane) 
available  for  nucleatlon. 


Equation  (3),  which  represents  a  change  from  previous  versions  of  the 
model, ^  appeared  to  be  required  to  correctly  represent  the  4340  shear 
band  data  described  earlier. 
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Nucleaclon  and  growth  both  absorb  plastic  shear  strain.  For 
the  model,  it  is  presumed  that  growth  requires  less  effort  than  nucleation 
hence  when  strain  Is  applied,  growth  is  permitted  first.  If  not  all  the 
strain  is  taken  in  growth,  some  is  available  for  nucleation  and 

some  new  bands  will  be  formed  according  to  Eq.  (3).  If  nucleation  does 
not  absorb  the  available  strain,  the  remainder  is  taken  homogeneously. 
Because  of  this  interaction  of  nucleation  and  growth,  and  the  limit  on 
growth  rate,  nucleation  appears  to  be  strongly  rate-dependent  and  non- 
llnearly  related  to  the  amount  of  plastic  strain  in  the  contained  frag¬ 
menting  cylinder  calculations,  even  though  Eq.  (3)  is  actually  rate- 
independent. 

2.  Growth 

The  growth  of  shear  bands  follows  the  viscous  growth  law  used 
in  all  the  previous  models,  plus  the  restrictions  that  the  growth  velo¬ 
city  not  exceed  the  shear  velocity  and  that  the  strain  taken  in  growth 
not  exceed  the  total  strain  Imposed.  Now  an  additional  constraint  has 
been  added:  afcer  the  damage  T  exceeds  1.0,  the  bands  no  longer  grow. 
Thus,  it  is  presumed  that  further  straining  after  full  damage  (t  ■  1.0) 
is  taken  by  sliding  of  fragments,  and  not  by  growth.  (No  nucleation  or 
homogeneous  plastic  strain  occurs  either.)  This  change  in  the  growth 
process  leads  to  the  flat  tops  on  the  n^  versus  x  and  r^^  versus  x  curves 
in  the  contalned-fragmenting-cy Under  simulations  (Figures  43  and  44). 

3.  Representation  of  the  Distribution 

Shear  band  size  distributions  tend  to  have  the  approximately 
exponential  form  shown  in  Figure  30  in  Section  III,  which  shows  the  4340 
steel  data  from  Experiment  2.  In  various  versions  of  the  shear  band 
model,  we  have  represented  these  distributions  either  by  a  single  expo¬ 
nential  or  by  a  series  of  points.  Neither  was  very  satisfactory,  so  the 
problem  of  representing  the  distribution  was  reexamined.  The  results  are 
a  new  multiple-size-group  method  and  a  single  exponential  method.  For 
the  multiple-size  group  representation,  there  is  also  a  provision  for 
transforming  to  an  exponential  for  displaying  the  results. 
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The  new  multiple-size  group  version,  described  in  Appendix  A, 
provides  for  an  exponential  between  points  of  the  distribution.  This 
representation  provides  high  accuracy  in  defining  the  shape  of  the  dis¬ 
tribution  and  also  In  the  computation  of  the  damage,  T.  (This  latter 
feature  is  an  improvement  over  the  previous  use  of  a  multiple-point 
distribution. ) 

Because  the  multiple-point  distributions  require  much  more 
computer  storage  than  the  single  exponential  versions,  we  made  several 
efforts  to  condense  the  distribution  to  a  single  exponential.  In  the 
most  recent  attempt,  described  in  Appendix  B,  the  requirements  for  the 
distribution  are  to  conserve  the  damage  T  and  to  represent  a  midpoint 
slope  of  the  true  distribution.  The  new  condensation  procedure  is  superior 
to  previous  ones,  but  it  still  tends  to  underestimate  the  damage.  However, 
the  procedure  should  be  useful  in  many  situations  where  high  fidelity  is 
not  required. 

To  aid  in  comparing  the  multiple-point  distribution  results 
with  experimental  data,  two  new  features  were  added: 

•  The  volume  distributions  are  transformed  to  surface 
distributions. 

•  The  surface  distributions  are  fitted  to  an  exponential 
form. 

The  volume  distributions  (number  per  unit  volume)  of  shear  bands  are  used 
in  all  the  shear  band  calculations.  However,  for  printing  the  damage,  a 
conversion  is  made  to  the  numbers  of  bands  that  would  Intersect  a  plane 
passed  through  the  material.  This  number  per  unit  area  is  obtained  by  a 
statistical  transformation  described  in  Appendix  C. 

Following  the  transformation  to  a  surface  distribution,  this 
distribution  is  fitted  to  an  exponential,  using  the  damage  T  as  a  weighting 
factor.  The  least  squares  fit  procedure  is  outlined  in  Appendix  C.  The 
fitted  distribution  parameters  are  then  readily  compared  with  the  experi¬ 
mental  data,  as  in  Figures  43  and  44. 

4.  Work  Hardening 

A  general  model  for  work  hardening  was  incorporated  into  SHEARS 
for  treating  homogeneous  straining  in  the  intact  material  around  shear 
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bands.  The  work-hardening  yield  strength  is  taken  to  be  a  function  of 
Che  scalar  plastic  strain,  in  the  conunon  power-law  form 

n 

Y  -  Aj(EP)  ^  (4) 

where  Y  is  the  current  yield  strength 

and  are  constants  pertaining  to  the  i^^  interval  of  the 
work-hardening  curve 

e  is  the  plastic  strain. 

The  yield  strength  versus  plastic  strain  function  is  inserted  as  a  series 
of  points,  and  it  is  presumed  chat  the  curve  between  the  points  has  the 
form  of  Eq.  (4).  The  model  is  described  in  detail  in  Appendix  D. 

Tills  work-hardening  model  has  been  Inserted  Into  SHEAR3  for 
shear  banding  calculations.  The  yield  curve  parameters  are  read  into 
SHEAR3  at  IniciallzaClon  with  the  ocher  special  shear  banding  input  data. 
The  yield  stress  calculation  is  made  in  SHEAR3  immediately  following  the 
stress-relaxation  calculation  (which  computes  C*^) .  Because  of  this  order 
of  calculations,  the  stress-relaxation  algorithm  is  always  using  a  yield 
strength  appropriate  to  the  beginning  of  a  time  increment  instead  of  the 
end.  This  lack  of  correct  centering  of  the  yield  strength  should  not 
cause  significant  errors,  however,  because  the  subroutine  subcycles  to 
avoid  strain  Increments  exceeding  about  0,22  and  very  little  change  occurs 
in  the  yield  strength  for  such  a  strain  increment. 

5.  Damage  Assignment 

When  shear  band  damage  occurs,  strengths  and  stresses  in  several 
directions  are  reduced.  The  current  shear  band  model  SHEAR3  contains  a 
complex  stress  relaxation  process  (outlined  in  Appendix  A  of  Ref.  1)  for 
relating  the  damage  to  the  stress  tensor  and  reducing  the  strengths  and 
stresses.  In  this  section,  we  review  the  procedure  for  assigning  the 
damage  to  the  stress  orientation  and  clarify  the  nomenclature  somewhat. 

The  purpose  of  the  darige  assignment  procedure  is  to  determine 
the  effect  of  a  band  or  distribution  of  bands  on  stresses  in  any  arbitrary 
orientation.  The  following  requirements  for  the  procedure  were  developed 


81 


for  compucaclonal  purposes  and  do  not  necessarily  reflect  physical 
requirements: 

•  The  bands  remain  circular.  Hence,  shear  strain  in 
one  direction  in  the  plane  causes  uniform  growth  in 
all  directions  in  the  plane,  and  strength  is  reduced 
equally  in  all  directions  in  the  plane. 

•  Only  the  size  and  orientation  of  the  band  have  an 
effect  on  stress,  not  the  growth  in  the  last  time 
Increment.  Hence,  only  Che  residual  strength  Is 
determined.  This  method  should  be  satisfactory  if 
the  work-hardening  modulus  is  much  less  than  the 
elastic  modulus. 

•  When  damage  is  assigned  to  three  orthogonal  planes, 
the  total  damage  on  all  three  planes  is  constant  and 
independent  of  the  orientation  of  the  damage.  This 
assignment  criterion  tends  to  make  the  damage  invariant 
under  coordinate  rotations. 

•  Shear  strength  and  stress  are  reduced  for  shear  bands 

such  that  S  2«  S22,  and  $-2  ii'  Figure  45  are 

reduced. 


These  four  requirements  have  been  met  by  considering  three  orthogonal 
planes,  such  as  the  1-2  plane  in  Figure  45,  one  at  a  time.  The  damage 
in  each  shear  band  orientation  is  assigned  to  each  of  these  three  planes 
using  Che  proposed  relation 


(5) 


where 


is  the  damage  on  any  plane  such  as  1-2  in  Figure  45 


T  is  the  damage  in  the 
P 


th 


shear  band  orientation 


is  the  angle  between  the  normal  to  the  Ij  plane  and 
Che  normal  to  the  shear  band  orientation. 


With  these  definitions,  iji  •  90®  for  the  angle  between  the  band  shown  and 
the  1-2  plane  (and  also  for  the  2-3  plane).  That  is,  T,-  “  1  and  T..  ® 

p  Zj 

However,  for  the  1-3  plane,  ■  0,  and  the  band  has  not  effect  on  the 
stresses. 

The  sine-squared  factor  In  Eq.  (5)  is  not  rigorously  derived 
but  proposed  by  similarity  Co  the  more  usual  cosine-squared  relations 
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FIGURE  45 


MA.7893-97 


SHEAR  (S)  AND  NORMAL  <a)  STRESSES  CN  A  BLOCK 
CONTAINING  A  SHEAR  BAND 
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and  by  its  ability  to  meet  the  four  requirements.  With  the  sine-squared 
factor,  the  sun  of  Che  damage  is 


13 


+  V,,  +  T,o  -  2t 


‘23 


12 


(6) 


This  Invariance  can  be  verified  by  considering  an  arbitrary  vector  with 
cosines  Sf,  m,  and  n  to  three  orthogonal  axes.  These  axes  can  be  the 
normals  to  our  three  planes  of  interest.  Then  the  sum  of  the  sines 
squared  is 


I  sin  -  (1  -  ih  +  (1  -  m^)  +  (1  -  n^)  -  2 
2  2  2 

because  £.  -t*  m  -t  n  ■  1.  This  invariance  can  also  be  verified  by 

examining  Table  A-3  of  Ref.  1  and  summing  rows  of  the  weighting  factor, 

.  2, 

sin  ip. 

Thus,  Che  damage  assignment  procedure  used  in  SHEAR3  appears 
CO  meet  qualitative  physical  requirements  and  also  an  invariance  criterion. 

C.  Procedure  for  Fitting  Shear  Band  Model  to  Data 

Based  on  the  analysis  of  the  4340  fragmenting  cylinder  data,  the 
following  tentative  procedure  has  been  developed  for  determining  shear 
band  model  parameters  from  the  data. 

The  cylinder  data  should  represent  a  range  of  damage  levels  and 
strain  races.  The  cylinders  are  sectioned,  and  then  shear  bands  are 
counted  at  several  positions  along  Che  cylinders  and  assembled  to  cumu¬ 
lative  size  distributions  as  in  Figure  46.  These  surface  counts 
(Figure  46b}  are  Chen  transformed  to  volume  counts  (Figure  46c)  usinp, 
a  standard  metallurgical  transformation  process  such  as  that  of  Schell^^*'* 
or  of  Seaman  et  al.^2  or  the  approximate  procedure  described  in  Appendix  E. 
The  volume  size  distributions  are  each  fitted  to  the  exponential  function 

Ng  -  exp  (-R/R^)  (7) 
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Shear  Bands 
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t 


I 


(c)  VOLUME  DATA 

MA-7893-98  I 

I 

FIGURE  46  STEPS  IN  OBTAINING  CUMULATIVE  SHEAR  BAND  DISTRIBUTIONS  ! 

FROM  CONTAINED  FRAGMENTING  CYLINDER  DATA 

i 
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parameters  of  the  distribution.  The  and  parameters  are  then  plotted 
versus  their  position  along  the  cylinder.  (The  plots  are  like  those  in 
Figures  41  and  42  for  n^  and  r^.} 

To  obtain  estimates  of  the  strains  in  the  cylinder,  a  computer  simu¬ 
lation  (without  damage)  is  made  of  each  cylinder  test.  The  final  radii 
and  elongation  of  the  simulated  cylinder  should  be  compared  with  the  ex¬ 
perimental  values  to  verify  that  the  explosion  process  and  material  pro¬ 
perties  are  being  represented  correctly  in  the  simulation.  From  these 
simulations,  the  values  of  plastic  shear  strain  on  the  k-th  shear 
band  plane,  are  derived  for  the  locations  along  the  cylinders  where  the 
counts  were  made.  Such  computer-generated  strain  histories  are  shown  in 
Figure  47.  Then  Che  experimental  parameters  and  are  plotted  versus 
the  computed  shear  strains  to  obtain  estimates  of  the  nucleation  and 
growth  race  parameters  for  the  shear  band  model,  as  shown  in  Figures  48 
and  49. 

With  these  shear  band  parameters,  the  cylinder  tests  are  simulated 
again  using  SHEAR.1  to  obtain  computed  values  of  and  R^  as  functions 
of  axial  location.  These  computed  values  of  N  and  R,  are  compared  with 

O  1 

Che  measured  values  at  each  axial  location  for  each  cylinder.  This  com¬ 
parison  may  suggest  revisions  in  the  growth  or  nucleation  parameters. 

If  so,  the  parameters  are  changed,  the  simulations  are  repeated  and  the 
results  are  again  compared.  Alternatively,  the  surface  distribution 
parameters  n^  and  r^  may  be  compared  to  determine  how  to  modify  the 
nucleation  and  growth  parameters. 

The  foregoing  procedure  was  used  in  a  more  expanded  and  less  direct 
form  for  the  4340  data  because  the  nucleation  and  growth  functions  as 
well  as  the  parameters  were  being  changed,  and  Che  shear  band  model  was 
being  developed  at  the  some  time. 

From  Che  nucleation  and  growth  equations,  it  might  appear  chat  there 
would  be  no  need  for  repeated  simulations  of  the  experiments  to  obtain 
the  shear  band  parameters.  However,  Che  model  is  written  with  an  inter¬ 
action  between  nucleation,  growth,  and  stress-relaxation  that  causes  some 
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NUMBER  DENSITY  N.  FROM  THE  MEASURED  VOLUMETRIC  DISTRIBUTIONS  —  number/cm^ 
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FIGURE  48 


COMPUTED  PLASTIC  STRAIN  ON  THE  MODE  1  PLANE 

MA-7a93-Cr9 


MEASURED  NUMBER  OF  SHEAR  BANDS  VERSUS  COMPUTED  PLASTIC 
STRAIN  FOR  DETERMINING  AN  INITIAL  ESTIMATE  OF  THE 
NUCLEATION  RATE  PARAMETER 
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PARAMETER  R,  FROM  THE  MEASURED  VOLUMETRIC  DISTRIBUTIONS 


nonlinearity.  As  nancloned  in  the  previous  section  under  nucleaclon, 
the  plestic  strain  is  first  taken  in  growth.  The  renalning  strain  causes 
nucleaclon  of  new  bands.  If  some  strain  has  not  been  absorbed  in  growth 
and  nucleaclon,  chat  remainder  is  taken  homogeneously.  The  developing 
damage  causes  stress  relaxation  chat  may  alter  adjacent  damage  growth. 
Therefore,  multiple  aimulacions  of  the  experiments  are  required  to  obtain 
a  reasonable  match  to  the  experimental  shear  band  data. 

D.  Plans  for  the  Shear  Band  Model  Refinement 

The  observations  of  damage  In  the  long~rod  penetration  into  RHA 
reported  in  Section  II  and  our  experience  in  the  simulations  reported 
below  in  Section  VII  have  led  us  to  examine  some  areas  of  the  shear  band 
model  chat  should  be  modified.  RHA  shows  significant  anisotropy  in  the 
production  of  shear  bands:  bands  form  much  more  readily  in  the  plane 
normal  to  the  short  transverse  direction  than  in  the  other  planes.  Also 
we  notice  Chat  the  coalescense  and  fragmentation  processes  have  not  been 
treated  in  sufficient  detail.  Here  we  describe  these  problem  areas 
further  and  outline  approaches  coward  their  solution. 

1.  Anisotropy 

The  anisotropy  of  RHA  probably  plays  an  important  role  in  the 
penetration  process.  Shear  bands  appear  to  form  normal  to  the  short 
transverse  direction  (in  the  plane  of  the  place)  at  about  20Z  plastic 
strain,  but  do  not  appear  in  the  plugging  direction  until  about  60% 
plastic  strain  occurs.  The  type  of  anisotropy  observed  might  be  tenr.ed 
orchocropy  because  only  one  plane  differs  from  all  others.  This  anisotropy 
affects  the  production  of  bands,  but  does  not  directly  affect  the  elasticity 
relations  or  the  yielding  process  in  RHA.  The  introduction  of  this  aniso¬ 
tropy  into  SHEAR3  will  require  several  changes: 

•  Addition  of  damage  orientations  to  account  separately 
for  bands  in  the  weak  plane  and  normal  to  it  (these 
two  sets  of  bands  are  treated  as  one  set  now). 

•  Provision  for  different  nucleatlon  and  growth  processes 
for  each  band  orientation. 


90 


•  A  revision  of  Che  band  rotation  process  to  account 
for  Che  fact  that  the  angular  inclination  between 
bands  initially  is  90^,  but  gradually  approaches 
smaller  values  during  large  deformation. 

•  A  complete  change  of  Che  yield  and  stress  relaxation 
calculations. 

The  first  three  changes  are  considered  minor  because  they  add  only  a 
special  treatment  of  one  shear  band  orientation.  The  changes  in  the 
stress  calculation,  however,  are  major.  We  envision  a  complete  change 
of  Che  stress  computation  portion  of  Che  subroutine.  Instead  of  deter¬ 
mining  Che  principal  stress  directions  and  using  these  for  the  yield 
calculations,  we  will  use  a  stress  orientation  chat  has  one  axis  lying 
in  the  weak  plane.  The  yield  calculations  will  be  performed  using  three 
Mohr  stress  diagrams:  one  lying  in  the  weak  plane  and  two  in  planes 
orthogonal  to  that  plane  and  to  each  ocher.  These  Mohr  di?  -.rams  will 
not  necessarily  contain  the  principal  stresses  and  therefo’"'  provide  only 
an  approximation  to  the  yield  process.  To  improve  thu  accuracy  of  yielding 
calculation,  we  will  add  one  or  more  intermediate  planes  and  construct 
Mohr  diagrams  on  these  planes  also.  Initial  tests  of  this  approximation 
technique  indicate  chat,  with  four  Mohr  planes,  the  stresses  will  be 
accurate  within  5X  of  the  yield  strength  in  the  worst  case. 

We  expect  the  anisotropic  treatment  of  yielding  to  requir  ^  a 
somewhat  smaller  amount  of  code  Chan  Che  previous  treatment,  but  one  addi¬ 
tional  orientation  of  shear  bands.  Hence  the  next  version  of  SHEAR3  will 
not  be  significantly  larger  or  harder  to  use. 

2.  Coalescence  and  Fragmentation 

The  coalescence  and  fragmentation  processes  in  SHEAR3  will 
require  some  additional  effort  to  complete  them  in  a  physically  realistic 
manner.  Currently,  the  subroutine  ends  with  a  size  distribution  of  shear 
bands.  What  is  needed  is  either  an  added  section  in  SHEAR3  or  a  post¬ 
processor  (a  program  such  as  MOTT,  which  is  run  after  the  penetration 
calculation)  to  monitor  the  intersection  of  bands  to  form  fragments  and 
to  produce  a  fragment  size  distribution.  The  first  step  is  to  borrow 
features  from  the  MOTT  program;  however,  further  effort  should  be  directed 
toward  understanding  and  quantifying  the  intersection  process  leading  to 
fragmentation. 
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VI  SYNTHESIS  OF  A  COMPOSITE  CODE  FOR  PENETRATOR  SIMULATIONS 


To  predict  penetration  and  back-surface  fragmentation  in  a  penetrator- 
target  encounter,  it  is  necessary  to  have  a  computational  procedure  for 
simulating  the  materials  and  geometry  of  a  penetration  event  in  consider¬ 
able  detail.  From  phenomenological  studies  of  penetrators  and  targets 
(Section  II)  and  from  preliminary  simulations  of  projectile  penetrations 
(to  be  discussed  in  Section  VII)  we  have  generated  a  set  of  requirements 
for  a  computational  procedure.  Here  we  list  these  requirements  and  review 
the  capabilities  of  the  available  codes  that  have  been  used.  A  choice  is 
made  of  a  basic  code  (HEMP)  plus  elements  from  other  codes  to  enhance  its 
capabilities.  Some  current  results  with  the  Composite  HEMP  are  shown  to 
Indicate  the  state  of  development. 

A.  Requirements 

The  basic  problem  considered  is  the  penetration  of  a  long-rod  missile 
into  one  or  a  series  of  armor  plates.  The  target  may  be  defeated  by 
formation  of  a  plug,  excessive  plastic  flow  allowing  perforation,  rear- 
surface  spall,  or  a  combination  of  these  effects.  The  rod  may  be  defeated 
by  excessive  plastic  flow,  shear  banding  and  fracture  near  the  tip,  or  by 
fracture  along  its  length.  The  initial  problem  we  are  considering  is  two- 
dimensional:  either  a  normal  Impact  of  an  axisymmetric  missile  onto  a 
plate  or  an  oblique  impact  of  a  slab  onto  a  slab.  Later,  the  possibility 
of  simulating  oblique  Impacts  of  symmetric  missiles  (a  three-dimensional 
problem)  should  be  investigated. 

Our  attempts  and  chose  of  others  to  simulate  the  penetration  encounter 
have  led  us  to  outline  several  basic  requirements  for  a  computer  program 
for  the  simulation; 

«  Slide  lines 

*  Rezonlng 

•  Material  model  for  fracture  damage. 
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The  slide  lines  permit  materials  to  slip  past  each  other  or  to  impact  and 
separate.  The  penetrator  nose  and  sides  must  be  able  to  slide  along  the 
upper  surface  of  the  target  and  on  the  sides  of  the  crater  as  it  forms. 

In  addition,  the  target  may  be  severed  to  form  a  plug:  the  Interface 
between  the  plug  and  the  remainder  of  the  target  requires  a  slide  line 
if  the  relative  motion  is  significant.  The  length  and  direction  of  the 
severing  line  is  determined  during  the  penetration,  so  the  slide  lines 
must  be  able  to  unzip  gradually  and  to  wander  through  the  finite  difference 
or  finite  element  grid.  The  large  distortions  experienced  by  both  pene- 
trator  and  target  during  a  simulation  indicate  a  capability  to  rezone  or 
reconfigure  the  computational  grid  at  Intervals  during  the  calculation. 

The  material  property  model  is  essential  for  determining  the  ballistic 
limit  and  for  predicting  the  fragment  spray  at  the  end  of  the  event. 

Because  such  a  model  requires  many  variables,  the  penetration  code  must 
be  configured  to  permit  large  data  arrays  for  cells  in  which  the  material 
is  undergoing  damage. 

Besides  these  basic  requirements,  there  are  a  number  of  nominal 
requirements  that  make  a  computer  program  readily  usable: 

•  Input  in  a  free-field  format,  like  the  AMMRC  version  of  HEMP 

•  Finite  element  numbering,  as  in  EPIC 

e  Usable  on  available  computers 

•  Documentation. 

The  convenience  of  a  simple  input  form  is  important  because  it  may 
greatly  reduce  the  time  required  to  lay  out  a  new  problem  and  it  makes 
the  problem  definition  more  apparent  to  other  users. 

The  finite  element  numbering  method  allows  the  cells  to  be  configured 
in  a  more  general  way  than  in  rows  and  columns.  This  numbering  is 
convenient  for  discretizing  nonrectangular  or  multiply  connected  bodies. 

It  seems  essential  for  providing  slide  lines  that  extend  in  both  x  and  y 
directions,  as  along  the  nose  of  a  projectile  and  then  along  the  side. 

The  computational  logic  of  a  self-directed  slide  line  seems  more  straight¬ 
forward  with  finite  element  numbering  because  new  nodes  can  be  readily 
added  and  the  relationships  between  these  nodes  and  the  cells  can  be 
easily  updated. 
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The  code  must  be  In  a  standard  form  that  can  be  used  on  CDC,  CRAY, 
Univac,  Digital  (VAX),  and  other  computers.  This  transportability  means 
that  the  code  must  be  in  standard  FORTRAN  (without  the  special  and  different 
advantages  of  each  computer  vendor),  use  a  fairly  small  memory,  and  not 
depend  on  availability  of  special  peripheral  devices. 

The  code  must  be  well  documented.  Documentation  Includes  a  manual 
with  derivations  and  sample  input  and  a  code  with  ample  comment  statements 
and  equations  written  with  the  same  variables  as  the  manual. 

B.  Review  of  Available  Codes 

A  review  was  made  of  many  codes  that  have  been  used  for  penetration 
calculations.  The  review  was  limited  to  general-purpose  wave  propagation 
codes  written  for  large-deformation  problems  in  solids  treated  by  complex 
material  models.  Here  only  some  of  the  special  advantages  or  features  of 
these  codes  are  described. 

The  Eulerian  codes  HELP  and  HULL  are  able  to  treat  large  distortion 
problems,  although  this  advantage  leads  to  difficulties  in  treating 
sharply  defined  boundaries  and  slide  lines.  Also  there  tends  to  be 
smearing  of  cell  properties  such  as  the  fracture  quantities. 

The  Lagrangian  finite  element  code  EPIC  has  both  two-and  three- 
dimensional  versions  and  thus  is  particularly  suited  to  oblique  impacts. 

The  triangular  cells  resist  large  distortions,  but  falsify  the  stresses 
computed  from  pressure-sensitive  stress-strain  relations.  A  simple, 
frictionless  slide-line  treatment  is  available  In  the  code. 

STEALTH  is  also  Lagrangian  and  has  two—  and  three-dimensional  versions. 
It  features  automatic  rezonlng  and  slide-line  capabilities. 

The  Lagrangian  code  TOODY  has  been  used  in  many  penetration  calcu¬ 
lations.  It  has  a  discrete  rezoner:  wave  propagation  calculations  are 
halted,  the  entire  mesh  is  reconfigured,  and  then  the  calculation  continues. 
The  slide  line  treatment  is  particularly  accurate.  The  cells  on  each  side 
of  the  slide  line  are  accounted  for  in  detail  so  that,  if  no  slip  occurs, 
the  presence  of  the  line  should  be  undetected.  Stick  or  frictionless 
conditions  and  Coulomb  friction  properties  are  provided  on  the  slide  line. 
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The  slide  line  Is  tied  to  the  usual  row-and-column  schene.  The  number  of 
variables  per  cell  can  be  readily  Increased,  but  all  cells  are  assigned 
the  same  number  of  variables. 

WAVEL  Is  a  Lagranglan  code  with  advanced  slide  line  and  continuous 
or  automatic  rezoning.  It  has  been  developed  for  penetration  problems 
and  extended  to  a  quasl-three-dlmenslonal  form  for  oblique  impacts. 

HQIP  Is  the  forerunner  of  many  two-and  three-dimensional,  finite- 
difference  Lagranglan  codes  used  In  penetration  calculations.  It  has  a 
discrete  rezoner,  as  TOODY  has,  and  a  slide  line  capability  that  Is  tied 
to  the  rov-and-column  cell  numbering.  A  recent  version  is  HEMPDS,  which 
has  a  double  slide  line  (intersecting  slide  line)  capability,  but  no 
rezoner.  The  double  slide  provision  required  a  major  restructuring  of 
HEMP,  so  It  is  not  a  minor  addition  nor  an  element  that  can  be  readily 
transferred  to  another  code.  The  HEMP  code  has  many  versions,  each  with 
different  capabilities.  The  version  available  at  AMMRC  has  an  input 
provision  termed  a  NAMELIST  READ  statement.  This  provision  permits  all 
numerical  quantities  to  be  Identified  by  name  and  all  data  to  occur  in 
a  free-fleld  format  (no  fixed  positions  on  a  line).  This  greatly  simplified 
input  provision  speeds  the  laying  out  of  new  problems  and  minimizes  errors. 

DYNA2D  is  a  large  finite-element  code  that  has  been  used  for  simu¬ 
lating  self-forging  fragments  and  penetrations.  It  has  a  capability  for 
a  slide  line  that  can  extend  both  horizontally  and  vertically  at  the  same 
time.  Elements  of  higher  order  than  the  constant  strain  elements  in  HEMP 
are  available  for  use  on  elastic  or  other  small-distortion  problems. 

The  small  two-dimensional  Lagranglan  code  TROTT,  which  we  have  used 
in  many  of  our  applications,  has  the  advantage  of  simplicity  In  adding 
new  material  models  with  large  numbers  of  variables.  The  new  model  is 
added  by  the  Insertion  of  one  or  two  call  statements.  Extra  variables 
for  the  cells  treated  by  special  material  models  are  available,  and  the 
number  of  extra  variables  Is  designated  In  the  input .  The  number  of 
variables  per  cell  varies  from  cell  to  cell  so  that  relatively  Inactive 
cells  can  have  a  mlnumum  number  of  variables,  while  cells  in  the  region 
of  fracturing,  for  example,  can  have  50  or  100  extra  variables.  TROTT 
also  has  automatic  rezonlng,  like  that  In  WAVEL  and  STEALTH.  A  simple 
frictionless  slide  line  treatment  Is  available. 
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The  review  of  available  codes  indicaced  to  us  that  there  are  many 
good  candidates.  Although  no  one  code  has  all  the  desired  features,  all 
the  desired  features  are  available.  Therefore,  a  synthesis  from  existing 
codes  seemed  desirable. 

C.  Synthesis  and  Testing  of  a  Code 

The  HEMP  code  was  chosen  as  the  basic  code  to  work  with.  Onto  this 
has  been  added  finite-element  numbering,  automatic  rezoning,  a  TOODY-like 
slide  provision,  free-field  input,  and  provision  of  large  numbers  of 
variables  for  cells  requiring  them.  These  features  are  outlined  in  this 
section. 

The  status  of  each  portion  of  the  code  is  noted  in  Table  3,  with  the 
source  of  the  portion.  The  subroutine  containing  that  portion  of  the 
code  is  listed  in  parentheses.  The  remaining  portions  to  be  debugged  are 
the  slide  and  rezoning  parts.  The  rezoner  is  a  translation  of  the  TROTT 
subroutine  REZONT,  which  has  been  tested  on  Impact  problems.  The  slide 
line  treatment  is  based  largely  on  the  HEMP  logic,  but  with  the  use  of 
some  concepts  from  TOODY  and  from  STEALTH. 

Three  test  cases  have  been  used  in  debugging  composite  HEMP.  Test 
problem  1  was  a  one-dimensional  Impact  of  a  steel  projectile  onto  a  steel 
target  at  1000  m/s.  The  Impact  stress  achieved  was  about  225  kbar, 
attended  by  the  usual  oscillations.  These  results  were  compared  with 
analytic  values  to  verify  the  accuracy.  For  test  problem  2,  a  two- 
dimensional  planar  Impact  was  simulated.  A  detail  of  the  mesh  distortion 
near  the  plane  of  Impact  is  shown  in  Figure  50.  The  results  were  examined 
for  correct  stress  levels  under  the  impactor  (like  the  one-dimensional 
case)  and  physically  reasonable  motions  elsewhere.  Test  problem  3  repre¬ 
sents  a  blunt-nosed  projectile  fired  into  a  ligid  circular  trough  as 
shown  In  Figure  51.  This  problem  was  constructed  to  test  the  boundary 
conditions  and  the  time  step  control  under  severe  cell  distortion.  A 
corresponding  TROTT  calculation  showed  very  similar  results. 

Test  problem  3  is  being  used  now  to  test  rezonlng.  Comparable 
rezonlng  runs  have  also  been  made  with  TROTT,  so  we  have  an  indication 
of  the  results  to  expect. 
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Table  3 


STATUS  OF  COMPOSITE  HEMP 


Portion 
of  Code 

Purnose 

Source 

Status 

Input 
(in  GENR) 

Read  the  input  data 

HEMP  (BRL  and 
AMMRC  versions) 

In  use 

Setup 
(in  GENR) 

Lay  out  and  initialize 
arrays 

HEMP,  with 
finite  element 
logic 

In  use 

CYCLE 

Wave  propagation 
calculations 

HEMP 

In  use 

VQP 

Call  stress-strain 
models 

HEMP,  TROTT 

In  use 

Material 

Models  (VQP) 

Compute  stresses  at 
each  cycle 

TROTT 

All  models  used 
in  TROTT  are  in 
but  not  tested 

Slide 

(CYCLE) 

Permit  slide  lines 

TOODY,  HEMP 

Written,  not 
tested 

Double  slide 

Intersecting  slide 
lines 

HEMPDS 

Planning 

Re zone 
(REZONH) 

Modify  mesh  to 
minimize  cell 
distortion 

TROTT,  STEALTH 

Written,  in 
debugging 

Grid  listing 
(CYCLE) 

List  status  of  all 
cell  and  node 
quantities 

HEMP 

In  use 

Histories 

(EDIT) 

List  histories  of 
all  cell  and  node 
variables 

HEMP,  TROTT 

Written  and 
tested  for 
some  variables 
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y-distance. 


Projectile 


MA  .7893-1 36 


FIGURE  50  CONFIGURATION  AND  DETAIL  OF  RESULTS  FOR  TEST 
PROBLEM  2  USING  COMPI^ITE  HEMP 
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FIGURE  51  CONFIGURATION  FOR  TEST  PROBLEM  3  AND  COMPUTED  RESULTS 
AT  37  lU  (300  cyciwl.  USING  COMPOSITE  HEMP 
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Following  the  rezonlng  test,  we  plan  a  slide  line  test  with  a  geometry 
like  that  In  test  problem  2,  hut  the  projectile  block  will  have  a  tangential 
velocity  component.  The  next  test  will  be  a  normal  Impact  near  the 
ballistic  limit  and  with  the  material  represented  by  the  shear  band  model. 
This  fifth  test  will  Initially  Include  only  a  slide  line  between  the 
projectile  and  target,  but  will  later  include  rezonlng. 

Five  special  features  of  the  composite  HEMP  (Input,  finite-element 
numbering,  time  step  control,  rezonlng,  and  slide  lines)  are  described 
In  greater  detail  In  Appendix  F. 
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VII  SIMULATIONS  OF  PROJECTILE  PENETRATIONS 


Preliminary  projectile  penetration  calculations  using  the  shear  band 
model  were  made  with  the  TROTT  two-dimensional  wave  propagation  computer 
program,  an  SRI  program  that  is  similar  to  HEMP.  These  computations  were 
conducted  to  study  the  feasibility  of  performing  such  calculations  with 
a  wave-propagation  code,  locate  computational  problem  areas,  and  aid  in 
planning  an  augmented  code  with  the  capability  to  readily  handle  the 
simulations.  In  addition  we  wished  to  assess  the  possibility  of  pre¬ 
dicting  the  observed  fracture  r  td  fragmentation  using  a  shear  band  or 
ocher  fracture  model. 

Two  penetration  calculations  were  made  with  blunt-nosed  cylindrical 
projectiles  of  4340  steel  impacting  normally  onto  rolled  homogeneous 
armor  (RHA)  targets.  The  penetrators  were  10.16  cm  long  and  1.016  cm 
in  diameter  and  the  targets  were  2. 54  cm  thick.  In  both  cases  the  4340 
was  simulated  with  shear  band  parameters  obtained  from  our  fit  to  the 
fragmenting  round  data,  as  described  in  Section  V.  The  RHA  was  also 
simulated  with  the  shear  band  model  using  a  sec  of  parameters  derived 
mainly  from  the  4340  set,  but  with  some  modifications  to  fit  obvious 
characteristics  of  Che  RHA.  The  strain  threshold  was  raised  from  15!(^ 
for  4340  to  20Z  for  RHA  to  represent  the  initiation  of  shear  banding  on 
the  weakest  plane  (normal  to  Che  short  transverse  direction) .  The  initial 
yield  strength  for  4340  was  1.03  GPa  and  for  RHA  it  was  0.53  GPa.  Work- 
hardening  functions  derived  trom  quasi-static  tensile  data  were  provided 
for  each  steel. 

The  penetration  calculation  without  a  slide  line  was  stopped  at  17  us 
because  excessive  cell  distortion  had  led  to  infinitesimally  small  time 
steps  in  the  calculational  cycle.  The  final  configuration  of  the  target 
and  projectile  is  shown  in  Figure  52.  The  figure  also  contains  contours 
of  equivalent  plastic  strain  e^.  Although  the  shear  band  development  on 
a  particular  plane  is  a  function  of  individual  tensor  components  of  strain, 


103 


Che  equlvalenc  plastic  strain  Is  a  good  taeasure  of  the  general  level  of 
shear  band  damage.  In  the  penecrator  there  la  a  concentration  of  plastic 
strain  near  the  impact  plane  and  essentially  no  distortion  in  the  rest  of 
the  rod.  The  target  shows  large  amounts  of  plastic  strain  at  the  edge  of 
Che  crater  and  just  beneath  the  penetrator.  In  the  target  under  the 
penecrator,  there  is  little  radial  variation  of  strain  until  a  radius 
equal  to  Chat  of  Che  penetrator  is  reached;  Chen  Che  contours  change 
direction.  From  this  calculation,  we  can  note  that  the  front  3S%  of  the 
penetrator  and  considerable  material  in  the  target  around  Che  crater  would 
be  fragmented.  Because  the  projectile  still  has  considerable  forward 
momentum,  the  final  state  of  penecrator  and  target  cannot  yet  be  estimated 
from  the  calculation. 

In  the  second  penetrator  calculation,  a  slide  line  was  provided  within 
the  target  at  a  radius  equal  to  that  of  the  penetrator.  As  the  shear  band 
criterion  for  initiation  was  reached  in  adjacent  cells,  the  predetermined 
slide  line  was  activated  and  sliding  was  permitted.  In  this  calculation, 
it  was  necessary  to  Increase  the  penecrator  velocity  to  2000  m/s  to  avoid 
the  impact  of  the  sides  of  Che  penetrator  onto  the  upper  surface  of  Che 
target.  Figures  S3  and  54  show  Che  configurations  of  penetrator  and 
target  at  17  us  (to  compare  with  Figure  52) and  at  23  us,  respectively. 

The  calculation  was  halted  at  23  us  because  of  cell  distortion.  At  this 
time  Che  slide  line  had  propagated  from  Che  impact  point  to  the  last  node 
at  Che  back  of  the  target.  s 

The  geometry  and  strain  levels  in  Figures  52  and  53  can  be  compared 
directly  because  the  conditions  and  time  differ  only  by  the  presence  of 
the  slide  line  in  the  second  calculation.  The  amount  of  penetrator  damaged 
in  the  two  cases  and  Che  level  of  damage  are  nearly  identical.  However, 
with  the  slide  line,  the  penecrator  has  gone  nearly  twice  as  deep  and 
retained  a  nearly  flat  front.  In  contrast  the  penetrator  in  Figure  52 
is  heavily  distorted  and  the  front  is  very  curved.  For  the  calculation 
with  a  slide  line,  the  target  crater  is  much  larger  in  radius  as  well  as 
in  depth.  Also  there  is  a  discontinuity  in  strain  contours  along  the 
slide  line.  The  strain  contours  in  the  target  in  Figure  53  also  differ 
from  chose  in  Figure  52  because  of  the  greater  depth  of  the  penetrator. 
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Impact  Velocity  =  2000 


[♦I? 


At  23  MS  (Figure  54),  the  penecrator  is  0.3  cm  deeper  into  the  target 
than  it  was  at  17  jis.  The  plastic  strain  contours  have  extended  deeper 
into  the  target,  but  are  essentially  unchanged  in  the  penecrator. 

The  major  problems  that  arose  during  the  calculations  were 

•  Cell  distortion,  which  led  to  very  small  time 
steps  and  hence  effectively  halted  the  calculation. 

•  Lack  of  an  appropriate  criterion  for  propagation 
of  the  slide  line. 

•  No  provision  for  contact  between  the  target  and  i:he 
sides  of  Che  penecrator. 

The  problem  of  cell  distortion  can  probably  be  treated  by  a  discrete  or 
continuous  (automatic)  rezoner,  which  will  reconfigure  the  cells  to  keep 
them  approximately  square.  The  slide  line  propagation  problem  became 
apparent  because  there  was  not  a  strong  concentration  of  plastic  strain 
or  shear  band  damage  along  the  path  Che  slide  line  was  expected  to  take. 
This  problem  may  be  diminished  by  using  smaller  ceils  so  that  steeper 
strain  gradients  can  be  maintained.  Resolution  of  this  problem  will 
probably  require  careful  study  of  both  the  slide  line  process  and  the 
shear  banding,  as  well  as  Che  coupling  between  the  two. 

The  third  problem  was  Che  tendency  for  the  free  surface  of  the  target 
to  rise  and  eventually  coincide  with  the  bulging  sides  of  the  penecrator. 
Hence  there  is  a  need  for  a  slide  line  between  the  front  and  sides  of  the 

penecrator  and  the  top  of  the  target.  In  our  two  calculations,  the 
penecrator  had  a  flat  front  and  so  the  computation  began  with  contact 
between  penecrator  and  target.  No  separation  or  sliding  between  the  two 
was  permitted.  In  the  penetration  calculation  with  a  slide  line  in  the 
target,  the  impact  velocity  was  artificially  raised  to  avoid  Impact  of  the 
projectile  sides  onto  the  target. 

Thus,  a  complete  calculation  appears  to  require  a  rezoner,  a  slide 
line  between  the  sides  and  front  of  the  penecrator  and  the  top  of  the 
target,  and  a  propagating  slide  line  in  the  target.  All  these  features 
are  available  in  some  wave  propagation  codes,  although  not  in  the  same 
code,  so  the  performance  of  a  full  penetration  calculation  seems  feasible. 
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To  gee  an  Indication  of  Che  fragment  size  distribution  chat  would  be 
computed  if  the  penetrator  completely  perforated  the  target,  we  examined 
Che  state  of  shear  banding  at  23  ys.  The  MOTT  program  described  elsewhere^ 
was  used  for  Che  transformation  from  shear  bands  to  fragments.  For  this 
calculation,  the  damage,  T^,  is  described  by 


T 


f 


By^  Tp  Z  Z 
*  m  n 


N  R  ^ 
mn  mn 


where  B  is  the  ratio  of  number  of  fragments  to  number  of  bands,  y  is  the 
ratio  of  fragment  radius  to  band  radius,  T^  Is  a  constant  factor  related 
to  fragment  aspect  ratio,  N  is  the  number  of  shear  bands/cm^,  and  R  is  the 
radius  of  a  shear  band.  The  subscripts  m  and  n  refer  to  the  orientations 
and  the  sizes,  respectively  in  a  shear  band  distribution.  When  the  material 
is  fully  fragmented,  ■  1.0.  For  partial  fragmentation  (T^  <  1.0),  we 
use  as  Che  fraction  of  cell  mass  that  Is  fragmented.  When  these  rules 
were  used  In  Che  MOTT  program,  total  masses  of  fragments  and  size  distri¬ 
butions  of  fragments  were  obtained  separately  for  Che  penetrator  and  target. 

The  computed  mass  fraction  of  the  penetrator  that  was  completely 
fragmented  at  the  time  shown  in  Figure  54  is  31%,  in  fair  agreement  with 
the  value  of  40%  obtained  by  measuring  the  mass  of  the  recovered  unfrag¬ 
mented  distal  portion.  Furthermore,  the  computed  mass  lost  from  the  target 
in  the  form  of  the  plug  and  collaterally  fragmented  material  is  35  grams, 
in  good  agreement  with  the  34  grams  obtained  by  weighing  the  armor  plate 
before  and  after  the  experiment. 

Comparison  of  these  computations  with  experimental  observations 
(transformed  from  Figure  40(a)  of  Section  IV)  in  Figure  55  shows  that  we 
compute  Coo  many  small  fragments  and  too  few  large  ones.  This  discrepancy 
can  have  several  causes.  First,  a  large  but  unknown  number  nf  small 
fragments  were  not  recovered  in  the  experiment,  as  evidenced  by  the  large 
number  of  nonperf orating  dents  in  Che  witness  plate  shown  in  Figure  38(b) 
of  Section  IV.  Second,  since  the  SHEARS  parameters  for  RHA  have  not  yet 
been  derived  from  Che  CFC  data,  in  the  computations  we  used  the  measured 
nucleatlon  threshold  strain  for  RHA,  but  the  nucleaCion  rate  function 
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FIGURE  54  CONFIGURATION  OF  TARGET  AND  PENETRATION  AT  23  IN  TROTT  SIMULATION  WITH  UNZIPPiivib  slIDE 
LINE.  INCLUDING  CONTOURS  OF  EQUIVALENT  PLASTIC  STRAIN 


CUMULATIVE  NUMBER  OF  FRAGMENTS  WITH  MASS  GREATER  THAN  M 
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FIGURE  SB  COMPARISON  OF  MEASURED  OOWNRANGE  FRAGMENT  SIZE  DISTRIBUTION 
FROM  EXPERIMENT  SHOWN  IN  FIGURE  $1  WITH  DISTRIBUTIONS 
CALCULATED  FROM  COMPUTER  SIMULATION 


appropriate  for  the  4340  penetrator  steel.  On  the  basis  of  our  experience 
with  another  ductile  steel  (Armco  iron),  we  expect  the  nucleation  rate  to 
be  less  for  the  RhA  than  for  the  4340  steel.  Furthermore,  the  calculations 
of  Figure  54  used  a  somewhat  higher  impact  velocity  than  was  achieved  in 
the  experiment.  Accounting  for  all  three  of  these  factors  would  reduce 
Che  differences  between  observation  and  computation.  An  additional  cause 
for  computational  error  is  the  as-yet  unmodeled  effect  of  the  significant 
anisotropy  of  the  RHA. 

The  above  uncertainties  mean  that  these  comparisons  of  computed  and 
observed  fragment  size  discributiuns  can  only  be  semiquantitatlve.  Never¬ 
theless,  the  computations  provided  the  Information  desired.  The 
overall  approach  is  shown  to  be  promising;  SHEAR3  parameters  obtained  from 
laboratory  CFC  experiments  appear  to  produce  encouraging  results  in  the 
penetration  computations.  At  an  impact  velocity  well  above  the  ballistic 
limit,  SHEAR3  will  produce  a  spray  of  target  and  penetrator  fragments  with 
a  reasonable  mass  or  size  distribution.  Work  is  continuing  to  perform  a 
similar  computational  simulation  of  the  case  in  which  a  plug  formed  at  a 
velocity  just  above  the  ballistic  limit  (Figure  S6  of  the  Summary). 

Further  calculations  must  then  await  the  determination  of  the  SHEAR3 
parameters  for  RHA  material. 


Ill 
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PUNS  FOR  THE  THIRD  YUR 


For  the  experimental  phase  of  this  program,  our  tentative  plans  for 
the  third  year  are  as  follows: 

(1)  Phenomenology  experiments.  No  new  phenomenological  studies 
are  planned.  We  feel  that  we  have  a  substantial  data  base 
for  understanding  the  mechanisms  of  penetrator  nose  erosion 
and  RHA  target  plate  penetration  at  the  velocities  and 
obliquities  of  most  interest.  Studies  by  Netherwood  of 
BRL^^  have  shown  that  breakup  of  the  distal  portions  of 
Che  penetrator  during  a  multiple-plate  impact  do  not  appear 
to  be  as  consnon  as  once  thought.  So  we  have  cancelled 
plans  for  tests  at  BRL  Co  study  this  problem. 

(2)  Material  characterization  experiments.  We  plan  to  conclude 
the  analysis  of  the  RHA  CFC  experiments  and  to  obtain  shear 
band  size  distributions  for  use  in  conjunction  with  our 
computer  simulations.  Ue  also  plan  to  begin  some  charac¬ 
terization  experiments  on  DU.  In  particular,  we  plan  to 
perform  quasi-static  tensile  tests  and  CFC  tests  on  the  DU 
specimens  Chat  have  been  machined  for  us.  The  tensile 
specimens  include  axial  specimens  and  radial  specimens 
with  and  without  centerline  porosity.  The  CFC  specimens 
include  only  axial  specimens;  that  is,  their  axes  coincide 
with  the  axes  of  the  cylinders  from  which  they  were  machined. 

(3)  Downrange  fragmentation  data  base  experiments.  No  new 
experiments  are  planned,  but  we  will  conclude  (with  the 
help  of  J.  Barb  of  BRL)  the  analysis  of  the  quarter-scale 
normal- Impact  test  and  obtain  final  fragment  size,  velocity, 
and  angle  distributions. 

The  plans  for  the  computational  portion  of  the  program  next  year 
Include  completing  the  development  of  a  method  for  predicting  fracture 
and  fragmentation  during  a  long-rod  projectile  impact.  This  method  is 
based  on  a  shear  band  model  that  suitably  represents  the  fracture  and 
fragmentation  process  and  on  a  wave  propagation  code  that  can  handle  the 
large  distortions  occurring  during  the  penetration. 

The  shear  band  model  will  require  development  of  an  anisotropic 
treatment  of  strength  and  a  more  detailed  coalescence  and  fragmentation 
process.  These  plans  were  given  in  detail  in  Section  VD.  In  addition. 
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we  plan  Co  fit  the  model  to  Che  shear  band  data  obtained  from  a  series  of 
RUA  contained  fragmenting  cylinders. 

In  the  development  of  Composite  HEMP,  the  next  step  is  testing  the 
rezoning  subroutine.  Probably  some  minor  adjustments  to  the  rezoner  will 
be  required  later  as  we  g-ln  experience  with  it.  The  slide  line  develop¬ 
ment  will  proceed  as  follows:  test  the  single  slide  treatment,  add  an 
unzipping  or  propagating  slide  line,  and  add  the  intersecting  or  double¬ 
slide  treatment.  Continued  effort  will  be  made  to  minimize  the  user  effort 
by  introducing  automatic  layout  features. 

Simulations  of  normal  impacts  will  be  continued  to  test  the  development 
of  code  features,  guide  the  insertion  of  new  features,  and  test  the  suit¬ 
ability  of  the  shear  band  model.  Oblique  impact  simulations  will  be  made 
i^ter ,  We  expect  these  simulations  to  guide  in  further  code  and  shear  band 
model  developments. 
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APPENDIX  A 

REPRESENTATION  OF  THE  SHEAR  BAND  DISTRIBUTION  BY  MULTIPLE  SIZE  GROUPS 
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The  shear  band  size  distribution  may  be  represented  by  a  series  of 
exponential  line  segments  as  shown  in  Figure  Al.  The  user  may  select  Che 
number  of  segments  or  nodes  and  thereby  determine  Che  precision  with 
which  Che  dlscrlbucion  is  represented.  This  section  describes  the  method 
for  specifying  the  distribution  and  treatment  during  nucleation.  growth, 
and  the  computation  of  damage. 

At  the  first  nucleation,  a  shear  band  distribution  is  created  with 
a  series  of  radii  Rj  as  abscissae  and  Mj  as  ordinates  at  the  nodes.  The 
number  of  Intervals  in  BFR(M, 34),  and  the  size  of  the  largest  R  ■ 
is  BFR(M, 28) .  The  smaller  radii  are  initialized  by  the  formula 

Rj  -  ^H1  (Al) 

^RATIO 

where  provided  in  the  initialization  portion  of  SHEAR3  as  1.4. 

The  individual  Nj  values  are  Initialized  from  a  nucleated  distribution  of 
the  form 

N  »  N  e”*^^*^n  (A2) 

n 

where  N  and  R  are  the  number  and  size  parameters  of  the  nucleated 
n  n 

distribution.  Then  Nj  Is  given  by 

N  -  (A3) 

j  n 

The  numbers  and  radii  are  stored  in  the  CN  array  in  the  SHEAR3  subroutine 
in  the  order 


N  ,  R-,  N, ,  R.,  N-. 
o’  1’  1  2’  2 


NS 


where  NS  ■  BFR(M, 34),  the  number  of  intervals. 

For  nucleation  on  cycles  following  initialization  of  the  disttibution, 
the  numbers  at  each  node  are  merely  incremented  as  shown  in  Figure  A2  (a) . 

Nj  -  N^o+  exp(-Rj/R^)  (A4) 
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Cn  N  —  8n  numbef/cm^ 


1 


No 


RADIUS  —  cm 
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FIGURE  A1  SHEAR  BAND  SIZE  DISTRIBUTION  REPRESENTED  BY  A  SERIES 
OF  POINTS  AND  EXPONENTIAL  LINE  SEGMENTS 
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where  N.  is  the  number  before  Che  time  increment.  Thus  the  nucleated 

Jo 

distribution  is  added  to  Che  old  distribution  at  Che  nodes  of  the  old 
distribuCicn.  Between  these  nodes  the  distribution  is  assumed  to  remain 
a  series  of  straight  line  segments.  This  segmented  distribution  is  only 
an  approximation  to  the  actual  sum  because  Che  sum  of  two  exponentials 
Is  not  in  general  an  exponential. 


Growth  of  the  bands  is  represented  by  an  increase  in  the  at  all 
Che  nodes  as  shown  in  Fig.  Alb.  The  rate  of  growth  is  given  by  the 
following  viscous  law.^ 

■‘"l  ‘"^^6 

J  -  C  R  — ^ 

dt  dt  (AS) 


where  C_  is  a  growth  coefficient  and  is  the  plastic  strain  on  the 

G  90 

plane  associated  with  the  shear  bands  under  consideration.  Eq.  (A5)  is 
used  with  Che  constraint  chat  the  growth  velocity  cannot  exceed  the  shear 
wave  velocity.  When  an  exponential  distribution  grows  according  to  Eq. 

AS,  it  remains  an  exponential.  Thus  the  straight  line  segments  exactly 
represent  the  growth  of  the  distribution  {see  Fig.  Alb).  However, 
if  the  limit  velocity  is  reached  over  some  part  of  the  distribution,  the 
new  distribution  is  no  longer  exponential  and  Che  line  segments  are  only 
a  discrete  approximation  to  the  grown  distribution  (see  Fig.  Ale). 

Thus  both  nucleation  and  growth  are  readily  treated  with  the  distribution 
represented  by  multiple  size  groups,  although  approximationsoccur  in  both 
computations. 

The  computation  of  the  damage  x  is  performed  exactly  under  the 
assumption  of  an  exponential  curve  between  nodes  of  the  distribution. 

Here  the  expression  for  this  curve  is  obtained  in  terms  of  the  coordinates 
at  the  adjacent  nodes  and  Chen  an  integration  is  performed  to  obtain  T. 
Within  each  Interval  the  distribution  is  given  by 

N  -  N  ,  exp  (A  ,R)  (A6) 

mj  mj 

Then  the  number  at  the  nodes  on  either  end  of  the  interval  are 

*^j“l  “  ”mj  ^*mj^j-l^ 
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(a)  SUMMATION  OF  A  PREVIOUS  (b)  GROWTH  OF  A  MULTIPLE 
DISTRIBUTION  AND  A  POINT  DISTRIBUTION 

NUCLEATED  DISTRIBUTION 


(c!  GROWTH  OF  A  DISTRIBUTION  NEAR  THE 
LIMITING  VELOCITY 

MA-7S93-102 

FIGURE  A2  NUCLEATION  AND  GROWTH  PROCESSES  FOR  A  SHEAR  BAND 
DISTRIBUTION  GIVEN  BY  A  SERIES  OF  POINTS 
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N.  =  N  ,  exp  (A  .R.) 
J  nj  mj  J 


Equations  (A7)  and  (A8)  can  be  used  to  solve  for  and  A^^ 


„  R .  inN  .  1  -  ,  in  N  . 

inN  .  -  1  1-1  1-1  1 

mj  d - ± 

"j  -  ^j-1 


in  (Nj_j/Nj) 


with  the  expression  for  the  curve  N  versus  R  defined  (Eq.  A6  with  constants 
from  A9  and  AlO),  the  integral  for  T  can  be  evaluated. 

The  damage  factor  T  in  any  orientation  is  defined  as  a  sum  over  the 
entire  distribution 

T  -  EAN  -cR  •  nR^  (All) 

where  AN  is  the  number  of  shear  bands  per  unit  volume  with  radius  R,  c  is 

a  proportionality  factor  (in  the  range  0.1  to  IrO;  input  as  BFR(M,29)], 

2 

and  ffR  is  the  area  of  the  band.  The  factor  cR  is  introduced  to  form  a 

volume  of  material  ttcR^,  which  is  associated  with  the  band.  For  an 

interval  from  node  j-1  to  node  j,  the  damage  t  is 


f  3 

Tj  ■  TTC  I  R'^ 


The  limits  on  the  integral  are  reversed  to  make  the  sum  positive.  To 
prepare  the  integral,  dN  is  obtained  as  a  function  of  R  from  Eq .  (A6) 
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r' 


Let  £  “A  ^R.  Then  the  integral  becomes 
mj 


J-1 

Tj  -  ^  I 


’’j 

frcN  ,  r 
mJ  I 

mj 


(AlA) 


This  Integral  has  the  solution 


Tj  -  e''  -  3?^  +  65  -  6] 

mj 


'j-1 


(A15) 


For  convenience,  a  new  quantity  Z  is  defined 

Z  -  -  (5^  -  35^  +  65  -  6)  -  e^(l  -  5  +  1/25^  -  1/65^)  (A16) 

Then  T  makes  the  form 


6itc  N 


mj 


A^ 

mj 


"j- 


(A17) 


The  Z  quantity  is  approximately  equal  to  one  and  5  Is  always  negative. 
The  first  and  last  Intervals  are  treated  somewhat  differently.  For  the 


first  interval  R 


,  ,  -  R  is  always  taken  as  zero  and  is  therefore  not 
j-1  o 

Z  ■  1.0.  For  the  last  interval  the 


stored  in  the  array.  Then  Z.  , 

J-1  o 

original  expression  for  T  (Eq.  All)  is  reexamined.  The  distribution 
shown  in  Fig.  A1  contains  a  number  of  bands  N 


NS 


at  R 


R^  and  then 


a  continuous  distribution  for  all  smaller  radii.  The  damage  quantity  is 
therefore 


‘NS 


TIC 


;^s-n 


NS 


R^dN 


Ss 


(A18J 
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TtC 


With  appropriate  transformations,  these  terras  on  the  right  can  be  combined. 


^NS  "  'Sis  'SlS^ 


^NS  N  .  exp  (f  „) 
~Y-  mj  NS 

A  . 
mj 


with  j  -  NS.  Now  Tj^g  becomes 


Sttc  N 


m1 


NS 


mj 


”  ^j-1  ■*"  ^NS 


Stic  N  . 

(Y.  -Z..,) 


mj 


j  j-1 


(A19) 


where  the  new  function  Y  has  the  definition 


Y  -  c^d  -  C  +  1/2C^) 


(a20) 


With  the  use  of  Eqs .  (A17)  and  (A19),  the  damage  T  can  be  computed  for 
the  multiple  size  distribution. 

The  approximations  noted  above  in  the  assumption  of  an  exponential 
curve  in  each  interval  of  the  distribution  gives  some  indication  of  the 
number  of  intervals  required.  A  single  interval  with  three  data  items — 
N^,  Rj^,  and  N^^ — can  be  used  but  will  lead  to  too  large  values  of  N  and  T 
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Three  or  four  intervals  (7  or  9  variables)  are  probably  sufficient  if 
the  ranges  in  radii  are  selected  so  that  they  span  the  distribution  well 
(RATIO  may  be  modified  to  assist  in  the  representation)  •  In  our  simu¬ 
lations  of  the  fragmenting  round  data  to  determine  the  shear  band 
parameters,  we  used  ten  Intervals. 
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APPENDIX  B 

COMBINATION  PROCESS  TO  MAINTAIN  THE  DISTRIBUTION  AS  A  SINGLE  EXPONENTIAL 
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The  experlmencal  size  dlscrlbucion  for  shear  bands  has  approxlmacely 
the  shape  of  a  single  exponential.  Therefore,  Ic  should  be  possible  to 
use  a  single  exponential  Instead  of  the  multiple  size  groups  to  represent 
the  distribution  In  the  SHEAR3  computations.  The  purpose  of  this  appendix 
is  to  develop  a  procedure  for  computing  the  parameters  of  the  single 
exponential  In  a  way  that  provides  reasonable  fidelity  both  to  the  shape 
of  the  true  N-R  distribution  and  to  the  value  of  t.  Such  a  procedure 
has  the  advantage  of  requiring  only  3  variables  to  represent  the  distribu¬ 
tion  Instead  of  the  7  to  21  required  for  the  multiple  size  group  repre¬ 
sentation. 

In  these  calculations,  we  first  consider  the  distribution  changes 
associated  with  growth  and  nucleation,  and  then  the  means  for  combining 
the  growing  and  nucleated  distributions  into  a  single  exponential.  The 
damage  t  is  computed  and  comparisons  are  made  between  results  from  the 
single  exponential  and  from  a  raultiplc-size-group  calculation. 

The  distribution  is  represented  by  three  parameters:  (the  inter¬ 

cept  at  R  “  0) ,  (the  shape  parameter),  and  R^  (the  maximum  radius). 

N  -  exp  (-R/Rj^)  1  'Sj 

Growth  of  this  distribution  is  treated  by  allowing  both  Rj^  and  Rj^  to  grow 
according  to  the  viscous  growth  law  (Eq.  A5) ,  with  the  shear  wave 
velocity  limit  on  the  growth  rate. 

The  first  step  for  nucleation  is  to  compute  a  nucleated  distribution 
with  the  form  of  Eq.  (Bl) ,  but  with  the  nucleation  parameters  and  . 

N  -  AN  exp  (-R/R  )  R  <  (B2) 

n  n  —  bn 
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where  is  computed  according  to  the  nucleation  process,  and  and 
are  input  parameters. 

The  nucleated  and  growing  distributions  are  then  combined  to  form 
one  exponential  distribution  that  has  the  same  damage,  t,  and  represents 
the  general  slope  of  the  sum  of  the  two  distributions.  This  second 
requirement  was  made  specific  by  noting  chat  in  some  central  portion  of 
the  N-R  distribution  there  lS  a  major  contribution  to  the  r  integral. 

Hence  there  is  a  central  point  where  di/dR  is  maximum.  The  function  in 
the  vicinity  of  this  point  must  be  well-represented  for  the  combined 
distribution  to  reasonably  match  the  sum  of  the  nucleated  and  growing 
distributions.  The  two  requirements  are  then: 

•  The  slopes  dr/dR  in  the  conbined  distribution  and  the  sum 
reach  maxima  at  the  same  R  value. 

•  T  of  the  conbined  distribution  equals  t  of  the  sum. 

2  2 

For  a  single  exponential,  dr/dR  reaches  a  maximum  when  d  r/dR  equals  zero. 

dR  “  i  ”  • 

.2  N  R^ 

exp  (-R/P,)  (3  -  |-)  (B4) 

dR'^  h  ^  h 

Setting  d^x/dR^  -  0  gives  R  »  3R,  (B5) 


when  dx/dR  is  maximum. 


The  corresponding  expressions  for  the  sum  are 


.  5  N  AN 

10  n 


d^x  2  R^  ^10 

2  "  R —  ^~R/Rio) 

dR^  *^10  ‘'lO 


2 

+  (-3R^  +  |_)  ^  exp  (-R/R„) 
n  n 
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where  and  N  _  are  the  shape  parameter  and  intercept  number,  respectively, 

Xv  XU  2  2 

of  the  growing  distribution.  Setting  d  t/dR  to  zero  in  Eq.  (B7),  we  can 

obtain  a  value  of  R  ■  R  at  the  point  of  maximum  dt/dR  (the  solution 

c 

procedure  is  outlined  below).  Then  from  Eq.  (B5), 


Hence  the  combined  distribution  is  set  to  have  its  maximum  dr/dR  at  the 
same  radius  as  for  the  maximum  of  the  sum. 

The  condition  to  conserve  t  is  approximated  as 


N  R,' 
o  1 


^10  '^lO 


AN  R  ■ 
n  n 


in  which  the  cutoff  of  the  distribution  at  R  »  has  been  neglected. 

With  R^  known  from  Eq.  (B8) ,  Eq.  (B9)  can  be  used  to  determine  N^.  Thus 
the  two  parameters  and  R^  of  the  combined  distribution  are  determined. 

The  solution  of  Eq.  (B7)  for  R  »  R^  is  obtained  starting  with  a  change 
of  variable.  Let 


so  that  Ht  the  new  unknown,  is  a  small,  dimensionless  quantity.  Eq.  (B7) 
then  cakes  the  form: 


(1-n)—  -  1 
n 


AN  R,- 
n  10 

Nin^ 
10  n 


exp  -3(l-n)(^  -  1) 
n 


(BID 


This  relation  is  readily  solved  by  iteration  and  the  equation  is  cast 
in  the  form  appropriate  for  the  solution  method  used.  The  initial 
estimate  of  n  is  zero.  This  value  is  used  on  the  right-hand  side  and 
the  second  estimate  is  obtained  as  the  value  on  the  left-hand  side.  This 
process  is  repeated  until  successive  values  of  n  differ  by  less  chan 


Usually  only  one  to  four  Iterations  are  required  b'cause  n  is 


small. 
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Three  test  calculations  were  made  to  compare  results  of  the  current 
combination  procedure  with  those  from  a  multiple-point  size  distribution 
and  from  a  previous  combination  procedure.  A  shear  strain  of  43X  was 
imposed  in  37  equal  steps  to  reach  the  size  distribution  shown  in 
Figure  Bl.  The  variation  of  t  values  during  the  calculations  are  shown 
in  Figure  B2.  The  multiple-point  results  are  probably  the  most  accurate. 
The  current  combination  procedure  gives  a  distribution  (Fig.  Bl  that  is 
tangent  to  the  multiple-point  curve  in  the  vicinity  of  the  maximum 
dx/dR,  but  provides  a  smaller  t  (Fig.  32).  The  previous  combination 
procedure  gives  a  very  poor  indication  of  both  the  shape  of  the  distribu¬ 
tion  and  the  value  of  T. 

Comparisons  of  the  conbination  and  multiple-point  procedures  have 
also  been  made  for  complete  fragmenting-cylinder  simulations.  The 
combination  procedure  gave  significantly  lower  amounts  of  damage.  These 
comparisons  of  results  with  the  current  combination  procedure  and  with 
the  multiple-point  procedure  suggest  that  the  combination  procedure  should 
be  used  only  for  rough  approximations.  Probably  the  shear  band  parameters 
should  be  determined  with  the  multiple-size  version. 
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CUMULATIVE  NUMBER  OF  BANDS  GREATER  THAN  R  —  number  hix? 


MA-7893-103 


FIGURE  B1  SHEAR  BAND  SIZE  DISTRIBUTIONS  COMPUTED  BY  SHEARS 
AT  43%  PLASTIC  STRAIN  USING  DIFFERENT  PROCEDURES 
FOR  representing  THE  BAND  SIZE  DISTRIBUTION 
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DAMAGE,  T  “  jr/R-^dN 


APPENDIX  C 

VOLUME -TO -SURFACE  TRANSFORMATION  FOR  SHEAR  BANDS 
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In  our  model  calculations  of  shear  bands,  the  bands  occur  as  number 
densities  or  number  per  unit  volume.  However,  the  shear  bands  are  observed 
experimentally  as  a  number  that  intersects  a  surface  and  are  therefore 
counted  as  a  number  per  unit  area.  To  make  a  comparison  between  the 
computed  and  measured  values,  a  transformation  must  be  made.  In  our 
study  it  is  convenient  to  transform  the  volume  quantities  to  surface 
quantities  as  part  of  the  operation  of  the  SHEAR3  shear  band  subroutine. 

The  following  analysis  includes  a  derivation  of  the  equations  and 
the  solution  procedure  used.  The  method  differs  from  Scheil ' s^*^  in  that 
the  radius  Intervals  are  of  variable  slxe  in  the  present  analysis. 

Following  the  transformation,  the  surface  distribution  is  fitted  to  an 
exponential  for  ease  in  comparing  with  experimental  data. 

The  desired  transformation  relates  the  areal  counts  N  to  the  volume 

a 

counts  N^.  For  the  transformation  both  functions  are  discretized  into 
counts  that  pertain  to  a  radiu.s  interval.  The  relation  between  the 
counts  is  then 

J 

N^(i)  (Cl) 

j-1 

where  N^(i)  is  the  number  of  circles  per  unit  area  with  radii  between 
r^  and  is  ih®  number  of  bands  per  unit  volume  with  radii 

between  and  is  the  transformation  matrix  to  be  derived,  and 

J  is  the  number  assigned  to  the  upper  limit  of  observed  radii. 

In  addition  to  the  discrete  distribution  quantities  N  and  N  ,  we 

a  V 

assume  the  existence  of  continuous  distribution  functions  C(r)  and  P(R), 
the  areal  and  volume  densities,  respectively.  These  are  defined  such  that 
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{C2) 


N  (1) 
a 


N^(l) 


r 


i 


i+l 

p(R)dR 


R 


1 


(C3) 


For  Che  derivation  we  define  C'(r,R)drdR  as  the  number  of  bands  of  radius 
R  per  unit  area  that  intersect  the  plane  of  polish  with  a  circle  of 
radius  r,  and  C"(R)dR  as  the  total  number  of  bands  of  radius  R  to  inter¬ 
sect  the  plane  of  polish.  The  densities  S'(r,R)  and  C"(R)  are  related 
by  a  probability  function  P; 

5'(r,R)  -  P(r,R)  C"(R)  (C4) 


where  P  is  the  probability  of  r,  given  R.  This  probability  can  be  readily 
derived* from  Figure  Cl.  The  differential  probability  for  a  slice  with  a 
thickness  dh  is 

dP  -  P(h)dh  ■  P(r.R)dr  -  ^  - 

R 

R 


P(r,R) 


r  <  R 


-  0 


for  r  >  R 


(C5) 


The  areal  density  5(r)  is  related  to  the  volume  density  of  bands  p(R)  by 
noting  chat  Che  plane  of  polish  intersects  all  circles  of  radius  R  whose 
centers  lie  within  R  of  the  plane,  that  is,  in  the  volume  2R»A  (A  »  area) . 
For  a  unit  area 

C"(R)  -  2R  p(R)  (C6) 
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Conbine  this  with  Eqs.  (C4)  and  (C5)  to  get 

CD  00 

5(r)  -  J  5'(r,R)dR  -  J  P(r.R)C'(R)dR 

o  o 


„  <C7) 

J  P(r,R)2Rp(R)dR 
o 


The  areal  counts  are  numbers  between  two  radii;  that  is,  from  Eq .  (C7), 

2  J  j  P(r.l 


/.i+l 

N^(l)  -  J  C(r)dr 


.R)Rp(R)dRdr  (C8) 


"i  ° 


The  order  of  integration  is  reversed  and  the  integral  is  separated  into 
four  components  as  follows: 

li+l  Ji  5l+l  ^i+1  ?  lft+1  ri+1  ^  k+1 

J  J  V  J  J  J  J 


r^  o  o  ri 


ri  r^  r^  R 


’^i+l  ’^i 


(r)  (R)  (R)  (r)  (R)  (r)  (R)  (r)  (R)  (r) 

The  first  and  third  of  the  integrals  on  the  right  are  zero  according  to 
the  second  part  of  Eq.  (C5) .  The  resultant  equation  is 


N  (!)  -  2 

a 


r/ 


==7  p(R)drdK  (R  <  .  .  . , ) 
2  2 
R  -  r 


‘i+1 


^i+1  ^ 


— — 

i  J  -  I 


p(R)drdR  (R  > 


-  2 R^  -  r^  p(R)dR  +  2  J  (  Jr^-  r^  R^  -  r^^^^  p(R)dR 


i+1 


(C9) 


The  procedure  for  computing  the  is  essentially  the  same  as  Scheil ' s 

Rewrite  Eq.  (C9)  in  ths  form  of  Eq.  (Cl),  using  Eq.  {C9)  to  define  the 
a's.  Then  using  Eq.  C3  t  define  N^(l),  we  obtain 
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(CIO) 


and 


’'j+1 


/ 


P(R)  (^/  ~  dR 


(Cll) 


Ij 


J+l 

J  p(R)dR 


r . 

3 

To  make  the  compucation  of  the  oi's  tractable,  we  use  the  approximation 
of  Scheil  and  let  p(R)  be  constant  in  each  radius  interval.  Then  the 
a's  become 


“ii 


Vo2  2 
R 


dR 


’^i+l  "  '’i 


and 


“ij  ■  ‘J  (V'"  -  'i  -  V"'  - 


’'j+1  ■  “"j 


(C12) 


(C13) 


These  integrals  can  be  evaluated  to  obtain 
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Vl  ^j+1.1  “  ""iVl-i  "  ""j  "  Vl  Vl.  1+1 


2  2 
’^i+l  ’’j+l.i+l  *  ^j,i+l  ‘  *’i+l  ^j,i+l 


where 


-  '^j> 


C  M  2 

=j.i-  V'l-'i 
■'j.i  •  ‘"*'J  *  Y'J  ■  '■^ 

The  corresponding  expression  for  Is 


(C14) 


“ii  ■  'i*i  'i 


(CIS) 


The  system  of  equations  specified  by  Eq.  (Cl)  has  the  form 


N  (J)  »  a,  -N  (J) 

a  J » J  V 


N  (J~l)  -  a,  ,  ,  ,N  (J-1)  +  a,  ,  ,N  (J) 

a  J— 1,J— 1  V  J-itJ  V 


(C16) 


N  (1)  -  a,  ,N  (1)  +  a,  ,N  (2)  .  .  .  a,  ,N  (J) 

a  l.lv  1.2v  l,jv 


where  J  is  the  number  of  radius  intervals.  When  the  coefficients  a  are 
computed  and  the  volume  quantities  are  known,  the  values  are  easily 
evaluated  from  Eq .  C16,  From  the  values,  the  curaulative  numbers  C 
larger  than  the  given  radius  are  determined. 

J 

C(r^)  - 

The  C(rj^)  values  then  form  the  surface  distribution,  which  may  be  compared 
to  experimental  data. 
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Because  the  surface  data  in  the  form  of  C(r)  versus  r  tend  to  be 
approximately  exponential,  it  is  convenient  to  fit  the  points  to  an 
exponential  of  the  form 

C(r)  «  n  exp  (-r/r  )  (CIS) 

o  1 


where  n  and  r,  are  the  parameters  of  the  expoticntial .  The  fitting  uses 
o  1 

the  standard  least  squares  procedure  with  a  weighting  factor.  Let  in  n(r) 
be  approximated  by  the  series  k 


y(r) 


(C19) 


The  coefficients  a^^  are  obtained  by  minimizing  the  square  of  the  difference 

»  2,n  C(r^)  -  y(r^) 

at  the  nodes  r^  of  the  distribution.  That  is,  minimize 


J 

E 

i+1 


2 

“i^i 


in  n(r.) 

X 


I 


2 


(C20) 


Differentiation  of  this  expression  with  respect  to  and  a^  and  setting 
the  derivatives  to  zero  leads  to 


a  • 

o 


2 

Zw.in  c^r .)Iv .T .  -  Ew.r.in  C(r.)Ew.r^ 

X  X  11 _ IX  1  1  1 

Ew^Ew^r^  -  (Ew^r^)^ 


(C21) 


Iw  Iw.r.Zn  C(r.)  -  Lw.r.Ew  £n  C(r.) 
1  11 _ i _ 11  i _ 1 

Iw  Iw  r .  -  (Ew  r.)"^ 

Ill  1  1 


Taking  the  logarithm  of  Eq.  (C18)  and  comparing  the  result  with  Eq .  (C19) , 

we  can  identify  a  and  a,  as 
o  1 


a  »  in  n  or  n  -  exp(a  ) 

O  O  0^0 


3]^  “  *"1  *  ■ 


(C22) 
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To  emphasize  the  portion  of  the  distribution  that  contributes  most  to 
damage  in  computing  the  fit,  the  weighting  factor  was  chosen  as 


w 


i 


C(r^)r 


2 

i 


(C23) 


Thus  the  weighting  factor  is  effectively  T^. 

The  preceding  transformation  and  least  squares  fit  procedures  were 

inserted  into  the  printing  section  at  the  end  of  SHEAR! .  The  printing 

is  requested  as  part  of  edit  listings  on  some  computational  cycles.  Then 

the  volume  density  quantities  are  listed  as  usual,  plus  the  surface  densities, 

and  n  and  r,  . 
o  1 
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APPENDIX  D 

WORK-HARDEN'ING  MODEL 
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A  general  model  for  elastic-plastic  work-hardening  is  derived  here 
for  use  in  the  intact  material  in  the  shear  banding  model  SHEAR3.  Yield 
strength  is  taken  as  a  function  of  the  scalar  plastic  strain  e  only. 

For  the  model,  the  yield  versus  curve  is  divided  into  a  series  of 
strain  intervals.  Within  each  interval  (except  the  first),  the  yield 
strength  varies  with  the  power  law 

y  « 


In  the  first  interval  the  y  -  relation  is  linear. 


The  model  derivations  are  made  considering  the  elastic  loading  and 
unloading  paths  and  the  yield  function  shown  in  Figure  Dl.  The  ordinate 
0  is  the  equivalent  stress,  a  scalar  measure  of  the  distortional  stresses. 


a 


+  3a 


12 


(Dl) 


where  a'j  are  deviator  stresses.  The  abscissa  is  an  equivalent  strain 
£,  a  scalar  measure  of  the  distortional  strains. 


£ 


dE 


(D2) 


(D3) 


The  definitions  of  0  and  dE  are  written  so  that  o  is  the  axial  stress  and 
E  is  the  axial  plastic  strain  for  a  uniaxial  test. 
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e,  EQUIVALENT  STRAIN 
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FIGURE  D1  DEFINITION  OF  QUANTITIES  FOR  MULTISEGMENTED 
WORK-HARDENING  CURVE 


The  data  for  the  model  are  Inserted  as  a  series  of  yield  strain 
points  (Y^,  E^).  In  the  first  interval  the  yield  curve  has  the  linear 


Y  -  Y^  +  A^E 


In  all  subsequent  intervals  the  power  law  form  is  used 


Y  -  A^(£^) 


For  convenience  in  using  the  program,  only  the  yield  curve  points  are 
input  data:  the  parameters  and  n^  are  initialized  in  the  subroutine. 

During  initialization  the  points  in  the  yield  curve  are  read  in  and 
then  A^  and  n^  are  computed.  For  the  first  interval,  n.  *  1.0,  and 


Y  -  Y 

‘2  n 


and  is  required  to  be  zero.  In  subsequent  strain  intervals,  n^.  and 
A^  are  Initialized  by  comparing  the  yield  strengths  at  each  end  of  the 
strain  interval. 


1-1  ■  ''i(vi)  * 


Hence 


1 
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and 


A 


i 


(D9) 


With  these  preparatory  calculations,  the  computations  for  each  cell  are 
made  very  quickly  using  r.qs .  (D4)  or  (D5). 
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APPENDIX  E 


AN  ANALYTICAL  TRANSFORMATION  OF  SURFACE  COUNTS 
TO  A  VOLUME  SIZE  DfSTRI BUTTON 
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Surface  counts  of  voids  or  cracks  in  a  single  plane  often  lead  to 
cumulative  size  distributions  with  an  exponential  form.  For  this  case 
an  approximate  analytical  transformation  can  be  made  to  a  volume  size 
distribution  that  Is  also  exponential.  This  analytical  transformation 
Is  outlined  here. 

The  surface  crack  or  void  size  distribution  is  assumed  to  be  in  the 

form 


n(r) 


-r/r, 
n  e  1 
o 


(El) 


where 

n(r)  is  the  cumulative  number  per  unit  area  with  observed  radii 
greater  chan  r 

n  is  the  total  number  per  unit  area  in  the  distribution 
o 

r^  is  the  shape  parameter  of  the  distribution. 

Hence,  it  is  assumed  chat  the  count  data  have  been  fitted  to  Eq.  (El) 
to  obtain  the  parameters  n^  and  r^. 

The  desired  volume  size  distribution  has  Che  form 


N(R)  -  N  e“^'^^l 
o 


(E2) 


where 

N(R)  is  the  number  per  unit  volume  with  true  radii  greater  than  R 
is  the  number  per  unit  volume 
is  the  shape  parameter  of  the  distribution. 
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The  transformation  of  Eq.  (El)  was  performed  numerically  with  BABSl, 
a  computer  program  based  on  the  transformation  analysis  of  Scheil  ’ 

Then  the  resulting  N(R)  values  were  nondimensionalized  and  plotted.  These 
distributions  were  then  fitted  to  the  exponential  function  in  Eq.  (E2)  to 
determine  N  and  R, . 

O  1 

The  transformed  volumetric  distributions  were  nondimensionalized 

and  plotted  as  r,N(R)/n  versus  R/r, .  The  results  are  shown  in  Fig.  El 
1  o  1 

for  r^  •  2  and  r^  •  4;  The  two  sets  of  points  were  indistinguishable 
except  at  large  R/r^^.  The  differences  at  large  R/v^  are  probably  caused 
by  inaccuracies  in  the  BABSl  transformation. 

Having  the  computed  points  of  a  volumetric  size  distribution,  we 
next  fit  these  points  to  the  exponential  form  in  Eq.  (E2).  Because  of 
the  nonlinearity  of  the  volumetric  distribution,  subsets  of  the  distri¬ 
bution  were  selected.  For  a  fit  over  one  decade  in  N,  the  parameters 
are 


N  ■  0.94  n  /r, 
o  o  1 

Rj^  ■  0.50  r^, 

and  for  two  decades  in  N,  the  parameters  are 

N  ■  0.82  n  /r, 
o  o  1 

R^  -  0.655  r^ 

With  these  fitted  results,  an  approximate  transformation  to  a  volumetric 
size  distribution  can  be  made  very  quickly. 


154 


r/r^  or  R/r^  —  dimensionless  crack  size 

MA-7893-105 


FIGURE  El  VOLUME  COUNTS  N.H  COMPUTED  FROM  THE  AREA  COUNT  n  •  n^e  ‘""i 
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APPEMDIX  F 


FEATURES  OF  COMPOSITE  HEMP 


1S7 


Progress  lias  been  made  coward  developing  a  two-dimensional  wave 
propagation  code  for  penetration  calculations  by  combining  the  best 
features  of  several  codes.  The  basic  structure  is  that  of  the  HEMP 
code,  but  a  free-form  data  input  technique  like  that  available  with 
a  few  versions  of  HEMP  has  been  incorporated.  Other  additions  Include 
automatic  rezoning,  an  advanced  slide  line  treatment,  a  new  time  step  control, 
and  finite-element  numbering.  These  five  features  are  outlined  in  this 
appendix.  The  code  is  still  under  development,  so  it  is  anticipated 
chat  some  of  these  features  will  be  modified  later.  As  a  background 
for  these  features,  the  physical  and  numerical  basis  for  wave  pr.ipagation 
calculations  is  given  first. 

Wave  Propagation  Calculations 

The  HEMP  programs  are  all  based  on  the  solution  of  the  Lagrangian 
equations  governing  motion  of  a  continuous  medium.  The  solution  technique 
is  called  the  method  of  artificial  viscosity  because  of  the  introduction 
of  viscous  forces  to  permit  a  continuous-flow  computation  in  regions  of 
high  stress  gradients.  Such  regions  are  interpreted  as  locations  of 
shock  fronts,  although  no  discontinuities  occur  in  the  computed  flow 
field.  With  this  artificial  viscosity  method,  the  equations  of  continuous 
flow  can  be  used  ev<»ij'where  and  no  special  equations  are  required  at 
shock  fronts.  HEMP  uses  the  leapfrog  method  of  von  Neumann  and  Richtmyer^^ 
to  integrate  the  flow  equations.  This  method  provides  second-order 
precision  in  time  and  space  without  the  use  of  second  differences. 

In  the  following  paragraphs  we  first  introduce  the  Lagrangian  set 
of  partial  differential  equations  governing  one-dimensional  planar  flow. 

These  are  Chen  changed  to  integral  and  numerical  form,  and  those 
requiring  it  are  extended  to  multidimensional  form.  Then  the  solution 
procedure  is  outlined. 
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The  Lagrangian  partial  differencial  equations  to  be  solved  for 
one-dimensional  planar  flow  are 


3X^ 

-  U 

(velocity) 

(F-1) 

3U  ^ 

3^/h 

■-t  < 

^3t\ 

(momentum) 

(F-2) 

3e\ 

i 

1 

H 

QjJ 

0. 

(energy) 

(F-3) 

3X^) 
.3h/  , 

n 

(mass) 

(F-4) 

T  ■  P  +  o'  +  Q  (constitutive  relations)  (F-b) 

-  f (E,D, . . . ) 

where  X  “  Eulerian  or  current  position 
t  •  time 

h  ■  Lagrangian  or  initial  position 
U  “  particle  velocity 
C,  •  current  and  original  density 

T,  o'  ■  total  mechanical  stress  and  deviacor  stress  in  the  direction 
of  propagation 
V  »  1/p,  specific  volume 
P,  Q  “  pressure  and  artificial  viscous  stress 
E  ■  internal  energy. 

These  equations  relate  velocity  to  coordinate  motion;  provide  for  conser¬ 
vation  of  nooentum,  energy,  and  mass;  and  give  the  stress-strain  or 
constitutive  relations. 

Now  each  of  these  five  equations  is  excunined,  extended  to  two  or 
three  dimensions,  and  put  into  numerical  form.  The  velocity  relation, 

Eq.  (F-1),  is  expanded  to  a  vector  form 
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(F-6) 


where  1  runs  over  the  three  coordinate  directions  and  h  refers  to  all 
three  directions.  In  numerical  form  the  integral  of  Eq.  (F-6)  is  simply 


+ 


U 


n  +  1/2 

ij 


(F-7) 


Here  n  and  n+1  refer  to  times  t*'  and  t”**^  and  j  refers  to  the  node  number. 
Velocity  U  and  time  increment  At  are  labeled  to  indicate  that  they  are 
centered  midway  between  t”  and 


The  momentum  equation  (F-2)  takes  the  following  three-dimensional 

form: 


(F-8) 


Here  1  again  Cakes  the  values  x,  y,  and  z,  and  h  ,  h 

X  y 

Several  strategies  such  as  those  in  HEMP,  TOODY,  and 


h^  are  Lagrangian 


coordinates. 

STEALTH  are  available  in  the  literature  for  putting  Eq.  (F-8)  into 
suitable  numerical  form.  Instead  of  using  any  of  these,  we  return  to 
a  simple  integral  form  of  (F-8) : 


(F-9) 


which  is  clearly  F  »  Ma.  The  mass  M  is  one-half  of  the  sum  of  the  masses 
of  Che  cells  surrounding  the  node  j  as  indicated  in  Figure  FI.  The  summation 
over  F  is  constructed  from  the  1-th  components  of  the  products  of  the 
stresses  in  each  cell  times  the  area  (of  the  mass  M)  on  which  they  act. 

These  l-th  components  of  Che  forces  can  be  written  as  a  stress  times  a 
projected  area 
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FIGURE  FI  CONFIGURATION  OF  FOUR  CELLS  SURROUNDING 
NODE  j  AND  ITS  ASSOCIATED  MASS  M 


j  -  1  i  j  +  1 

LAGRANGIAN  POSITION 


MA-7893-134 

FIGURE  F2  GRID  FOR  LEAPFROG  CALCULATIONS  SHOWING  INCREMENTS 
IN  SPACE  AND  TIME 


Hence  we  have  effectively  found  a  discretization  for  the  derivatives  in 


Eq.  (F-8) 


j.  ik 


(F-11) 


Hence  In  the  indiclal  notation  of  Eq.  (F-7),  Eq.  (F-8)  becomes 


U  y  n-1/2  X 

^ij  “ij  M  i.k  ^ik»^^£k 


(F-12) 


Here  the  subscript  i  refers  to  the  four  cells  (in  two  dimensions) 
surrounding  the  j-th  node.  In  planar  two-dimensional  flow,  the  sum  will 
contain  eight  terms  for  each  direction  i.  The  time  increment  At  is 


n+1/2  n-1/2  +  At”"^^^ 

t  -  t  -  - 2 - 


(F-13) 


Although  (F-12)  is  obtained  from  a  finite  element  viewpoint,  it  agrees 
exactly  with  the  results  of  finite-difference  derivations  for  planar 
two-dimensional  flow.  Differences  arise  only  in  axisyrametric  flow. 

The  energy  equation  (F-3)  is  rewritten  in  tensor  form  as 


ii  ,  vi:  T  ^ 

3t  i.k  ik  3t 


(F-14) 


where  is  the  strain  tensor.  For  each  cell  this  expression  takes 
the  numerical  form 


“"ik"'' 


(F-15) 


Here  the  barred  quantities  are  averages 

r.n  v"  +  .  =n  t"  +  t"”^ 

V  -  - - -  and  T  -  - ^ - 


The  sum  in  Eq.  (F-15)  runs  over  four  terms  la  two  dimensions. 

The  expression  for  the  conservation  of  mass  (F-4)  takes  the  simple 
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n  ^  ^  _M_ 

P  ”  ‘^o  AV^  ”  AV^ 

where  Av”  is  the  cell  volume  at  time  c**. 

The  constitutive  relation  (F-5)  is  expanded  to 


(F-16) 


ik 


P  + 


(F-17) 


where  P  and  Q  remain  scalar  pressures. 

14 

In  the  leapfrog  scheme  of  von  Neumann  and  Richtmyer  ,  the  five 
fundamental  equations  are  not  solved  simultaneously,  but  in  a  prescribed 
order.  This  order  is  suggested  in  the  diagram  in  Figure  F2.  First  there  are 
the  cell  calculations  at  the  time  t”  and  mid-node  positions  j  -  1/2  and 
j  +  1/2: 


Ae 


These  are  followed  by  the  node  point  calculations: 

,^+1/2  yH+l 

j  '  J 

These  steps  are  undertaken  for  each  cell  and  node.  Then  the  time  is 
incremented  and  the  process  is  repeated.  In  composite  HEMP,  tnese 
calculations  are  performed  in  two  separate  loops:  one  over  all  cells 
and  one  over  all  nodes.  This  sepaiatlon  of  the  loops  was  dictated  by 
the  finite  element  numbering  procedure.  In  the  section  on  time-step 
control,  this  separation  of  cell  and  node  loops  will  be  seen  to  have  a 
significant  advantage  over  the  usual  combined  loop. 


Data  Input 

The  data  are  provided  to  composite  HEMP  in  a  free-field  form  like 
that  required  for  the  HEMP  version  at  AMMRC  and  the  version  sent  to  BRL 
from  SRI  in  1979.  In  this  form  each  number  is  provided  with  an  alphanumeric 
label.  The  data  items  can  be  liiserced  in  a  nearly  arbitrary  order. 

Position  of  the  labels  and  numbers  on  a  line  are  also  fairly  arbitrary. 

Some  sample  input  decks  are  presented  and  described  below,  the  rules  for 
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constructing  the  input  are  given,  and  the  method  used  in  the  code  for 
interpreting  this  input  is  outlined. 

Sample  Input 

The  data  input  is  grouped  under  several  type  labels:  TITLE,  GENERAL, 
MATERIAL,  BLOCK,  BOUNDARY,  HISTORY,  AND  REZONE.  Specific  input  parameters 
and  labels  are  provided  under  oach  of  these  type  labels  as  shown  in 
Table  F'-l.  An  example  of  this  input  is  shown  in  Figure  F3.  The  first  line  is 
the  TITLE  line  and  need  not  contain  the  word  TITLE.  All  other  data  must 
be  preceded  on  each  line  by  one  of  the  type  labels.  For  example,  on  the 
line  beginning  with  GENERAL  in  Figure  F3,  AMAT  is  the  label  for  the  material 
name,  RHO  is  the  label  for  the  density,  7.85.  The  next  lines  begin 
with  $,  the  continuation  symbol.  Some  labels,  such  as  XBC  under 
BOUNDARY,  can  precede  a  series  of  numbers. 

The  problem  described  by  the  data  in  Figure  F3  is  also  illustrated 
by  a  sketch  there.  A  steel  plate  is  impacting  a  circular  trough  from 
the  left.  The  BLOCK  data  provide  Che  configuration  as  one  quadrilateral 
contained  within  four  points  given  in  counterclockwise  order.  BLOCK 
also  specifies  the  initial  velocities,  number  of  cells  in  each  direction, 
and  the  material.  The  BOUNDARY  specifications  provide  the  line  of 
symmetry  (Y  ■  0)  and  the  configuration  of  the  trough  (given  by  10  points). 
Histories  of  several  quantities  are  requested  at  a  number  of  points  in 
the  plate.  The  data  Include  a  request  for  a  complete  rezcning  every  five 
cycles,  leaving  the  points  on  the  line  of  symmetry  and  the  rear  of  the 
projectile  untouched. 

Data  Input  Rules 

The  following  general  rules  have  been  developed  to  govern  all  input. 

•  Each  line  is  in  Che  form:  "typelabel  label  number  label 
number"  or  "typelabel  label  number  number  number" 

•  Labels  and  numbers  must  be  separated  by  one  o'*  more  of 
the  three  characters:  blank,  equal  sign,  and  comma. 

•  Numbers  may  be  provided  in  integer,  decimal,  or  exponent 
form  and  may  be  positive  or  negative.  Repeated  numbers 
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Table  FI 


Parameter 

INPUT  PARAMETERS 

Block 

(Type  Label) 

AND  DEFINITIONS 

Def initions 

A 

HISTORY 

.  2 

Area,  cm 

AM 

HISTORY 

Mass  assigned  to  a  nodal  point,  g 
or  g/cm 

AMAT 

MATERIAL 

Name  of  material  (alphanumeric) 

BC.BCN 

HISTORY 

Boundary  condition  indicators 

CLIN 

GENERAL 

Coefficient  of  linear  artificial 
viscosity,  nominally  0.05 

COMC 

HISTORY 

Array  for  cell  variables 

CQSQ 

GENERAL 

Coefficient  of  quadratic  artificial 
viscosity,  nominally  4.0 

D 

HISTORY 

Density,  g/cm^ 

DEXX, 

DEYY.DEZZ 

HISTORY 

Strain  increments  in  the  X,  Y,  or  Z 
directions 

DEXY 

HISTORY 

Tensor  shear  strain  increment  in 
the  X-Y  plane 

E 

HISTORY 

Internal  energy,  erg/g 

EQSTC 

MATERIAL 

2 

Bulk  modulus,  dyn/cm 

EQSTD 

MATERIAL 

Second  coefficient  in  series  for 
bulk  modulus,  dyn/cm^ 

EQSTE 

MATERIAL 

Sublimation  energy,  erg/g 

EQSTG 

MATERIAL 

Grvlneisen  ratio 

EQSTH 

MATERIAL 

Gas  Grvlneisen  ratio 

EQSTN 

MATERUL 

Exponent  in  expansion  equation  of 
state 
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Table  FI  (continued) 


Parameter 

Block 

(Type  Label) 

Def Initlons 

EQSTS 

MATERIAL 

Third  coefficient  in  series  for 
bulk  modulus,  dyn/cm^ 

ESC 

MATERIAL 

Material  property  array 

EXX, 

EYY.EZZ 

HISTORY 

Strain  in  the  X,  Y,  or  Z  direction 

EXY 

HISTORY 

Tensor  shear  strain  in  the  X,  Y 
plane 

FX.FY 

HISTORY 

Forces  in  the  X,  Y  directions 
assigned  to  a  nodal  point,  dyn  or 
dyn/cm 

G2 

MATERIAL 

2 

Twice  the  shear  modulus,  dyn/cm 

IDIR 

BOUNDARY 

Direction  of  definition  for  a  wall 
or  slide  line:  positive  for  master 
material  or  wall  on  left  when  facing 
in  direction  of  Increasing  node 
point  numbers 

IPRINT 

GENERAL 

Cell  and  node  listings  are  provided 
at  time  steps  that  are  multiples 
of  IPRINT,  nominally  10 

MAT 

BLOCK 

Number  of  the  material  in  the  quadri 
laterial 

MU 

MATERIAL 

2 

Shear  modulus,  dyn/cm 

NBCPT 

BOUNDARY 

Number  of  points  in  the  XBC,  YBC 
array  for  a  wall  condition 

NaCTYP 

BOUNDARY 

Boundary  condition  type: 

1  -  simple  control  on  X  or  Y 

3  -  slide  on  wall 

-  slide  on  wall,  can  leave  wall 

and  impac  t 

5  -  sticks  to  wall 

6  ~  sticks  to  wall,  can  leave  wall 

and  impact 
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Tabl'?  FI  iccntinucJ 


traaeter 

Block 

(Tvpe  Label* 

fc.  e  I  a  c.  A  K  « ■.  r*.  > 

NC12. 

NC23 

BLOCK 

Suat«r  or  b«ti.«:er.  Is;  *-.4 

2:».;  *n<j  b*tw««n  2'C  }r4 

poir.rs  tr.t  s-air  eral 

NFR£Q 

RE  CON 

Thv  reiv'fiin^  pattern  c-o  tour 
coRseciit  ;ve  c>-;l«s  begins  a;  t  ist 
steps  that  arc  ault  .ples  o:  NFRI  . 

NCEOM 

GENERAL 

Problest  geoaetry:  1  :or  planj:. 

2  tor  axis>-a»et r  : t 

NMAX 

GENERAL 

Maximuo  nuaber  o:  t  me  steps  ;r. 
the  calculation 

F 

HISTORY 

•i 

Pressure,  dyn/cm",  positive  in 
compi essicn 

PMIN 

MATERIAL 

Minimum  (tensile)  pressure  permitted 
a  spall  cutoff,  dyn/cm“ 

REZMIN 

REZCN 

Minimum  amount  of  node  motion  that 
is  considered  in  rezoning,  nominallv 
0.02 

RE2WT 

REZON 

Fraction  determining  the  amount  or 
computed  rezone  motion  to  apply  to 
a  node,  nominally  0.5 

RHO 

MATERIAL 

Initial  density  (distinguished  from 
RHOS  for  3  porous  material),  g/ctn^ 

RHOS 

MATERIAL 

3 

Initial  density  of  a  solid,  g/cm 

SBAR 

HISTORY 

«  2 

Equivalent  stress,  o,  dyn/cm 

SPRINl, 

SPRIN2, 

SPRIN3 

HISTORY 

Principal  stresses  in  first,  second, 
and  third  directions,  dyn/cm-, 
positive  in  tension 

SP2 

HISTORY 

2  2 

Sound  speed  squared,  cm  /sec 
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Table  FI  (continued) 


Parameter 

Block 

(Type  Label) 

Definitions 

STOPT 

GENERAL 

Maximum  problem  time  in  the  calcul* 
ation,  sec 

SXX, 

SYY.SZZ 

HISTORY 

Thermodynamic  stress  in  the  X,  Y, 
or  Z  directions,  dyn/cm^ 

TH 

HISTORY 

Rotation,  radians,  positive  counter¬ 
clockwise 

TRIQ 

GENERAL 

Coefficient  of  triangular  artificial 
viscosity,  nominally  0.02 

TXX, 

TYY.T22 

HISTORY 

Total  stresses  in  the  X,  Y,  or  2 
directions,  dyn/cm^,  positive  in 
tension 

TXY 

HISTORY 

Shear  stress  on  X  or  Y  planes, 
dyn/cm^,  tensile  positive  sign 
conversion 

WBFR 

MATERIAL 

Array  of  special  fracture  properties 

WCMP 

MATERIAI- 

Indicator  for  a  composite  material 
model 

WDS 

MATERIAL 

Indicator  for  a  devlatot  stress 
model 

WEQS 

MATERIAL 

Number  of  material  property  para¬ 
meters  used,  usually  at  least  30 

WFR 

MATERIAL 

Indicator  for  a  fracture  model 

WH 

HISTORY 

Condition  Indicator 

WMAT 

HISTORY 

Material  number 

WPOR 

MATERIAL 

Indicator  for  a  porous  material 
model 

WTSR 

MATERIAL 

Array  of  special  material  properties 
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Table  FI  (continued) 


Block 

Parameter 

(Type  Label) 

Def initions 

WVAR  MATKRIAL  Number  of  extra  variables  requested 

for  each  cell 


X,  Y 

HISTORY 

Nodal  positions,  cm 

XA.YA 

BLOCK 

Four  coordinates  of  a  quadrilateral 
listed  in  counterclockwise  order,  cm 

XAFECT, 

YAFECT 

REZON 

Nonzero  values  indicate  that  the  X 
or  Y  value  should  not  be  rezoned 

XAFECT 

BOUNDARY 

Type  of  boundary  condition  control 
used  with  NBCTYP  -  1: 

1  X  greater  than  or  equal  to  XBC 

2  X  less  than  or  equal  to  XbC 

3  X  equal  to  XBC 

4  or  0  X  free 

XBC, 

YBC 

BOUNDARY 

Point  or  points  defining  a  boundary 
condition,  cm 

XCONST, 

YCONST 

BOUNDARY 

X  and  Y  values  used  in  defining 
boundary  points  to  be  affected  by 
a  particular  boundary  condition,  cm 

XCONST, 

YCONST 

REZON 

X  and  Y  values  used  in  defining 
points  to  be  affected  by  a  rezoning 
restriction,  cm 

XD,  YD 

HISTORY 

Nodal  velocities,  cm/s 

XDNOT, 

YDNOT 

BLOCK 

Initial  velocities  of  a  quadrilateral 
cm/ a 

XO,  YO 

HISTORY 

Initial  nodal  positions,  cm 

YAD 

MATERUL 

2 

Work  hardening  modulus,  dyn/cni 
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Table  FI  (concluded) 


Block 

Parameters  (Type  Label) 


Definitions 


YAFECT  BOUNDARY  Type  of  boundary  condition  control 

used  with  NBCTYP  -  1: 


0  or  1  Y  greater  than  or  equal  to  YBC 

2  Y  less  than  or  equal  to  YBC 

3  Y  equal  to  YBC 

4  Y  free 


YC 

MATERIAL 

Initial  yield  strength,  dyn/cm 

YY 

HISTORY 

2 

Yield  stress,  dyn/cm 

Z 

HISTORY 

Mass,  g  or  g/cm 

2 


171 


c.  nCHf  TrsT  NJ.  3  pnojecTiue  ihpact  at  hoc  n/sec  into  a  circucan  trough 


CtNtRAC  lPRjNT*iO 
S  CLIn«0«05 


H'tAXalO  STOPTat.O 
COSQaa.O  TAtfiaO.aS 


lATCnlAU  A>l«TaSTEEI. 

eiSTC«USRE*l2 

TCauo3E«tO 


'E9SBSO.O 

EOSTOaS.TTE^lS 

iUaa.m^U 


HyAMauO 

eOSTSaS.17E«IS 


RHOaT.ai 

EOSTSat,** 


ILOCk  TO.NOTaO.O  X3N0Tal,5E»5 
i  XA  -«,aiTT«  13.9622*  13. 0*22*  •A.aSTT* 
t*  O.a  6.4  *.9  ••0 

:Kl2*i3  NC23a5  maT*i 


iOUNOABY  ubCTYPai 

TCONSTaO.O  XAPECTaO.O 


YAPBCTaS.O 


aUNOARY 


x^CTYPan  lOlRal 
XdC  3<6  5.0  9.0762  9.S 
X6C  13.6*33  16. 6 


NBCPTalO 

10.7BSI  11,8664  13.6*226  13. *572 


y8C  1 0.0  9.3  9.4972  7.7942  6.8944  9. 7Ml  4.0  3.0782 
Y8C  1.5628  0.0 


xC0*9r*l3.6*22*  XAFECTaO.6  YAPCCTaO.O 


lUUMOARY  NacTYP^a  lOtPai  leCPTaiO 

X8C  5.0  5.6  6.0762  9.5  10.7651  11.8940  13.0*226  13.6572 
XbC  13. 8633  14.0 

vec  10.0  9.0  6.4572  7.7962  *.6944  5.7651  4,0  3.0762 
Y6C  1.5*26  0,0 

YCONOTao.O  XAPCCTaO.O  YAPECTaO.O 


BOUnOaRV  NSCTYP*!  X8C  14.0 

0  XC0NS74t3.0«22*  YC0N8Ta»»0  XAPECTa2.0  TAPCCT*3.0 


HtaTORT  TlX  0  P  13.0  0  13.0  u.O  EXX  C0NC2  12.0  1.0  CUNNl*  XD  TO 
S  10.0  3.5  14.0  1.5 

REZOnE  NPRE0a5  REZNlNaO.Ol  RE2NTa0.5 

S  VC0M6TS0.0  Y*F€cTa2, 

*  XC0M5T*-4. 93774  XAFECT*!.  YAPECTaZ. 


MA-7893-1 12 

FIGURE  F3  DATA  INPUT  AND  GEOMETRY  FOR  A  STEEL  PROJECTILE 
IMPACTING  A  RIGID  CIRCULAR  TROUGH 
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may  be  inserted  with  a  star  as  the  repeat  symbol.  Thus 
3*2.6  means  three  values  of  2.6. 

•  A  dollar  sign  can  be  used  to  replace  a  repeated  typelabel. 

•  Abbreviated  labels  and  typelabels  may  be  used  if  they  are 
unique:  M  may  be  used  for  MATERIAL  but  BL  instead  of  B 
is  the  shortest  abbreviation  for  BLOCK  because  BOUNDARY 
also  starts  with  B. 

•  Expanded  labels  such  as  BLOCKTARGETRHA  may  also  be  used 
if  they  contain  no  separators. 

In  addition  to  these  general  rules,  there  are  special  interpretations 
used  for  input  under  each  typelabel. 

The  BLOCK  input  is  read  by  the  subroutine  BLOCK,  which  also  initializes 
the  X,  Y  arrays  for  cells  and  nodes.  This  layout  occurs  whenever  all  the 
data  for  a  block  have  been  read.  BLOCK  determines  chat  a  block  is  complete 
when  another  typelabel  BLOCK  is  encountered  (not  merely  a  repeat  symbol) 
or  when  additional  values  of  any  of  the  quantities  XA,  YA,  XDNOT,  YDNOT, 

MAT,  NC12,  or  NC23  is  encountered.  The  input  "SLIDE  1  TO  3"  permits  a 
disconnection  between  the  nodes  from  Che  first  Co  third  corners  of  a 
quadrilateral  block.  Thus  a  slide  line  is  provided  between  this  block 
and  a  previously  described  block. 

Two  kinds  of  data  are  provided  under  the  typelabel  BOUNDARY:  those 
governing  the  nature  of  the  boundary  condition  and  chose  describing  Che 
nodal  points  affected.  Boundary  types  are  designated  by  NBCTYF  with  the 
meanings  given  in  Table  FI.  NBCTYF  must  precede  all  other  BOUNDARY  data 
for  each  type.  The  labels  XAFECT  and  YAFECT  are  used  with  NBCTYF  1  to 
specify  Che  control  desired  for  X  and/or  Y  quantities.  XBC  and  YBC  are 
arrays  containing  coordinates  of  a  boundary  line.  The  individual  nodes 
or  lines  of  nodes  affected  by  the  boundary  conditions  are  provided  by 
X,  Y  values  labeled  XCONST  and  YCONST.  The  boundary  condition  indicators 
at  all  nodes  are  updated  following  the  reading  of  a  line  containing 
XCONST  and/or  YCONST.  Hence  the  other  boundary  data  such  as  XBC,  YBC,  and 
IDIR  should  precede  XCONST  and  YCONST.  Also  only  one  XCONST  and/or  YCONST 
sec  should  be  inserted  on  a  line.  Because  of  the  immediate  effect  on 
the  layout,  the  BLOCK  data  must  precede  BOUNDARY  data. 
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The  HISTORY  deca  permit  a  request  for  histories  of  cell  or  node 
quantities  at  a  series  of  X,  Y  locations.  Following  reading  of  this 
input,  the  aubroucln:t  HISTORY  selects  Che  cell  center  or  node  point 
nearest  the  designated  X  and  Y  point  for  each  history.  Note  that  the 
requested  quantities  need  not  be  array  variables.  The  requested 
quantities  are  scored  in  arrays  during  the  cell  and  node  loops  in 
subroutine  CYCLE  and  then  printed  as  histories  by  EDIT  at  the  end  of 
a  calculation. 

For  Che  MATERIAL  input,  AMAT  must  be  the  first  label  encountered 
for  each  new  saterial:  AMAT  triggers  the  counter  that  keeps  track  of 
the  material  number  and  storage  locations  allotted  for  the  material,  and 
it  initializes  some  indicators.  In  Figure  F4  the  input  deck  contains 
some  special  labels:  SHEAR  and  EXPLODE.  These  labels  call  for  the 
material  model  subroutines  by  the  same  names.  These  subroutines  read 
the  data  preceded  by  these  labels. 

In  the  MATERIAL  arrays,  the  two  iudicators  WEQS  and  WVAR  designate 
the  space  to  be  allocated  for  the  material  data.  WEQS  indicates  the 
number  of  equation  of  state  constants  (such  as  density,  moduli,  sound 
speed)  and  indicators  required  to  describe  the  material.  These  quantities 
are  scored  sequentially  in  Che  C0H2  array.  WVAR  gives  the  number  of 
additional  variables  required  for  each  cell.  Usually  there  are  19 
variables  per  cell  in  the  COMl  array:  these  include  a  yield  strength, 
cell  rotation,  and  a  condition  indicator,  as  well  as  mass,  stress,  density, 
and  other  quantities.  If  more  variables  are  required  to  describe 
quantities  such  as  porosity,  fracture  processes,  or  dislocations,  these 
are  provided  by  setting  WVAR  to  the  number  needed. 

Two  types  of  data  are  provided  to  govern  rezonlng:  chose  governing 
the  nature  of  the  rezonlng  and  those  designating  the  nodes  to  be  rezoned. 
The  REZONE  input  must  follow  the  block  data  because  WREZ,  a  nodal 
quantity,  is  set  following  reading  of  each  line  containing  the  XCONST 
and/or  YCONST  values.  Several  lines  of  XCONST  and  YCONST  values  can  be 
inserted  to  control  rezonlng  at  nodes  or  lines  of  nodes.  The  labels 
XAFECT  and  YAFECT  mean  that  X  and  Y,  respectively,  are  not  to  be  changed 
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C.  HCH*  FUtCMCNTINC  KOUNO  OF  FIFO  tTCCL*  IMOT  a 

SCNCRAL  IFKINTlJOe  NMAXaTOO  $TOFTal.E-« 

•  CLlNaO.l  COSQaa.C  TKtOaO.02 

a  NceoHai 

MAtCftlAL  AMATaFCTN  HCQSa|2  IIHOaO.514 

a  aPNat  aVAMa;  COaTCaO.ai 

a  CXFLOOE  NTYFlait  QCXPtal.ClO  XOCTaS.O  DlSTaO.S 

MATERIAL  AHATaaSaOSTCCL  MCOBaSa  RHOaT.BS 

a  wPRa)  avARa*) 

a  CQSTCal.smS  lOSTOaS.lTCia  tOSTSaS.lStlS 

a  COSTEaT.Saeia  lOSTCat.**  CQSThbO.ZS 

a  SHCAH  aPRa2t*o.*  la..  sa«,  a.ot.  i«.*  o.is.  o.or,  o.o« 

a  aiiEAR  aPR2R«.2s«  .aT7«  2.,  2.*  0..  lo. 
a  aHEAR  FNuc«i.*  i.«  i.f  I. 

a  aHEAR  NTTR£a2  NRANCta) 
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a  aHEAR  E3a0,.  a.O}«  i.o 

a  VCal.OSClO  MUaa.lREU 
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MATERIAL  AHATaatQOaTEEL  HCQaalO  RHOaT.BS 

S  HV«Ra2 

a  EOSTCai.S4El2  EOatOaS.|7Ei2  f OSTSaS. I7E 1 S 

«  CQBTEaT.JoElO  EOBTCat.oo 

a  TCao.EB  MuaR.tOEIl 
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a  CQSTCaS.OIEll  EaBTOaa.ooen  E2STSa2.02EI2 

»  EOSTEaO.lofo  fOSTCaE.E 

BLOCXEXP  XAao.aSS  IS.33S  li.3S5  0.0)5  haT*1 

a  TAao.  0.  2.2«  2.20 

a  NCUaao  NC2)ao 

SLUCKASao  XAbO.  15. 2«  15. 2a  0.  HATa2 

a  TAbE.EB  2.26  a.«2  ).«2 

a  MCl2a25  NC23al 

a  elide  I  TO  2 

OLOCkPmma  XAaO.  15.24  15.24  0.  HATa) 

a  rAa).42  1.02  4.37  4.37 

a  NCt2a25  nC23bi 

BLOCkETCEL  XAbO.  15.24  15.24  0.  MATaa 

a  TAaa.ST  4.37  10. lo  10.10 

a  NCt2a25  NC23a2 
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NC23a2 

15.28  0.  MATB5 

12,7  12.7 

boundary 

NBCTYPa] 

VCONETaO. 

YAPEC7a3 

boundary 

NBCTYPan 
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history 

a 

a 

C0MC22f 

COHCol. 

0*f  )•« 

C0MC24. 
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12..  3. 

C4NC25.  COMC40, 

COMC40  5..  3. 

MA-78B3-133 

FIGURE  F4  DATA  INPUT  FOR  A  RUNNING  DETONATION  IN  A  CONTAINED 
FRAGMENTING  CYLINDER 
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by  razonlng.  The  nuaerlcal  values  provided  for  XAFECT  and  YAFECT  have 
no  significance  here. 

Two  more  sample  input  dacha  are  shown  in  Figures  ^5  and  F6.  These 
are  for  the  first  two  te>.'t  cases  run  with  Composite  HEMP. 

Subroutines  to  Handle  Input 

The  free  field  input  form  is  handled  by  four  subroutines  (POP, 

FIND,  BESTMATCH,  and  UNPAGC)  plus  special  data- insertion  statements  In 
the  initialization  subroutines.  POP  reads  the  data  lines  one  character 
at  a  time,  examines  these  characters,  and  assembles  them  into  arrays  of 
labels  (LABEL)  and  numbers  (RLABEL).  FIND  is  called  to  compare  labels 
in  the  data  line  with  the  standard  labels  of  each  subroutine  to  determine 
which  variable  is  being  initialized.  If  no  exact  match  is  found,  FIND 
calls  BESTMATCH  to  examine  parts  of  the  labels.  BESTMATCH  first  calls 
UNPACK  to  separate  the  input  and  standard  labels  into  separate  characters. 
Then  BESTMATCH  counts  how  many  matches  there  are  between  characters  in 
the  input  label  and  each  of  the  standard  labels.  The  Initialization 
routine  controls  reading  and  calls  appropriate  routines  such  as 

BLOCK  to  accept  the  input  data.  BLOCK  examines  the  LABEL  values  with 
the  aid  of  FIND  and  assigns  the  numbers  in  RLABEL  to  the  correct 
constants  or  arrays. 

Discussion 

The  free  field  input  form  sketched  above  has  several  advantages 
over  the  traditional  fixed  location  format  used  with  most  wave  propaga¬ 
tion  codes.  This  format  is  expecially  convenient  to  use  on  a  video 
screen  where  column  positions  are  not  provided.  Only  the  required  data 
are  insertsd;  standard  data  are  provided  by  default  and  unneeded  variables 
are  omitted.  Becauai  all  data  are  labeled,  the  input  r^nds  to  be  easily 
understood  and  confusion  is  minimized.  The  use  of  abbreviations  and 
extended  names  provides  added  convenience  to  the  programmer.  When 
new  variables  are  needed,  they  can  be  readily  added  by  appending  new 
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C.  HEMP  PIB8T  1-0  RUN  IMPACT  STEEL  ON  STEEL 


GENEhal  IPRInT*!  nhaX*50  STOPT«1,0 
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FIGURE  FS  DATA  INPUT  AND  GEOMETRY  FOR  A  ONE-DIMENSIONAL  IMPACT 
OF  TWO  STEEL  PLATES 
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C.  hFHP  first  2*0  RUN  Impact  steel  on  steel 
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FIGURE  F6  DATA  INPUT  AND  GEOMETRY  FOR  A  TWO-DIMENSIONAL 
PLANE-STRAIN  IMPACT  OF  TWO  STEEL  PLATES 
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A4btli  to  ch«  Hats  available  and  adding  replacement  statements  [such 
aa  VARIABLE  ■  RLABEL  (LM3)1  In  the  appropriate  Initialization  subroutine. 

Finite-Element  Numbering 

With  finite-element  numbering,  each  cell  and  node  is  numbered  inde¬ 
pendently  of  Ita  position  in  the  grid.  This  contrasts  with  the  1,J- 
nuoberlng  common  to  most  finite  difference  codes  in  which  a  fixed  numbering 
relationship  exists  between  each  node  or  cell  and  its  neighbors.  This 
finite-element  numbering  allows  a  somewhat  more  flexible  layout  of  cells 
In  problems  with  complex  shapes.  It  is  essential  for  problems  in  which 
new  cells  or  nodes  are  created  as  by  wandering  slide  lines  or  by  general 
rexonlng  (grid  reconfiguration).  The  disadvantage  is  that  the  nodes 
around  each  cell  and  cells  around  each  node  must  be  Individually 
designated  because  there  Is  no  automatic  node-cell  relationship  as  In 
the  1, j-nufflberlng  system.  Here  we  describe  the  Indicator  arrays 
required  to  show  these  relationships  and  then  exhibit  some  cell  layouts 
that  are  possible  with  the  finite-element  numbering. 

Adjacency  Arrays 

During  layout,  the  nodes  and  cells  are  numbered  In  the  order  In 
which  they  are  constructed,  A  partial  node  and  cell  numbering  for  a 
quadrilateral  is  shown  In  Figure  F7.  During  layout  three  special  arrays 
are  constructed  to  provide  the  needed  adjacency  relationships.  The 
UNOD  array  of  four  integers  designates  the  nodes  surrounding  each  cell. 

For  cell  10  in  Figure  F7,  the  WNOD  array  is  12,  13,  8,  and  7.  The  nodes 
are  provided  In  counterclrckwlse  order. 

The  nodes  and  the  cells  around  each  node  must  be  available  to  the 
resoner:  these  Integers  are  In  the  WCELL  and  WNDND  arrays.  For  node  7, 
the  surrounding  cells  are  WCELL  -  1,  9,  10,  2,  and  the  nodes  are 
WNDND  ■  2,  6,  12,  8.  These  sets  of  cell  and  node  numbers  are  both  in 
counterclockwise  order.  Also  WNDND(i)  and  WNDND(1+1)  are  on  either  side 
of  WCELL(l).  For  node  6  only  two  cells  are  adjacent  so  these  two  arrays 
are 
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FIGURE  F7  A  PORTION  OF  A  GRID  SHOWING  INDEPENDENT 
NUMBERING  OF  NODES  AND  CELLS 
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WCELL  •  9,  1 
WNDND  -  11,  7,  1. 


In  BLOCK  these  two  arrays  are  filled  with  call  and  node  numbers  In 
ascending  order.  Then  ORDER  la  called  by  CENR  to  rearrange  the  elements 
of  each  array  to  produce  the  special  relationships  referred  to  above. 

Grids  Using  Finite-Element  Numbering 

The  finite  element  numbering  provides  some  added  flexibility  In 
laying  out  the  grid  or  In  changing  it  during  a  calculation.  Here  we 
describe  three  situations:  a  problem  where  a  varying  grid  size  is 
required,  a  computation  with  a  nonrectangular  object,  and  a  calculation 
with  a  moving  slide  line 

Often  It  la  of  Interest  to  provide  small  cells  In  one  region  of  a 
problem  and  large  cells  elsewhere.  For  accuracy,  the  cells  should  vary 
gradually  from  small  to  large  and  should  remain  nearly  square.  In  general, 
these  requirements  cannot  be  met  with  an  i,J -numbering  procedure  without 
Introducing  slide-line  discontinuities  In  the  grid.  Even  ch*hn  the  slide 
line  Is  restricted  to  one  direction  (1  or  j,  not  both)  so  the  cell  size 
can  be  varied  In  only  one  direction.  Figure  F8  shows  four  cell  configur¬ 
ations  that  permit  gradual  changes  In  cell  sizes  in  one  or  two  directions 
In  a  finite-element  numbering  scheme.  Because  of  the  Irregular  adjacency 
relations,  these  schemes  cannot  be  used  with  l,J-numbering  procedures. 

They  all  show  regions  In  which  the  cell  size  varies  lOZ  to  SOX:  this 
variation  may  be  too  rapid  for  wave  propagation  problems,  but  Is 
satisfactory  for  Intermediate  rate  problems  such  as  projectile  penetrations. 

Odd-shaped  objects  such  as  projectiles  and  spheres  are  difficult  to 
construct  with  nearly  square  cells  with  an  i,J-numberlng  technique. 

Figure  F9  shows  a  sphere  (modeling  a  human  head)  with  a  finite-element 
numbering  technique.  Note  that  all  cells  are  approximately  equal  In  area 
and  have  aspect  ratios  near  one. 

A  third  feature  of  codes  that  shows  the  need  for  finite-element 
numbering  is  the  slide  line,  especially  a  slide  line  that  propagates 
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(a)  A  One-Dir«ctional 
Call  Size  Variation 


(b)  One  Possible  Two-Directional 
Cell  Size  Variation  Derived 
From  the  One-Directional 
Variation 


(c)  Pattern  for  a  25%  Cell  Size 
Increase  Per  Cell  in  Two 
Directions 


(d)  Pattern  for  a  10%  Cell  Size 
Increase  Par  Cell  in  Two 
Directions 


MA-7893'127 

FIGURE  F8  CELL  SIZE  VARIATION  PATTERNS  THAT  MAY  BE  USED 
WITH  A  FINITE-ELEMENT  NUMBERING  TECHNIQUE 
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GA-31 4583-1 6A 


FIGURE  F9 


GEOMETRY  OF  HUMAN  HEAD  FOR  STUDY  OF  IMPACT  AND  CONCUSSION 
SHOWING  FINITE-ELEMENT  LAYOUT  FOR  CIRCULAR  OBJECT 
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during  the  calculation  along  an  irregular  and  unpredetermined  direction. 
Such  a  general  slide  line  is  illustrated  in  Figure  10.  As  the  slide 
line  moves  through  the  grid,  we  presume  that  it  may  move  to  any  adjacent 
node  along  its  general  direction  of  travel.  In  some  moves  it  merely 
separates  cells.  Then  one  node  such  as  A  must  be  divided  into  two  to 
permit  separate  motions  of  the  adjacent  cells.  However,  if  the  slide 
line  divides  a  quadrilateral  as  along  BC  into  two  triangles,  additional 
nodes  are  required  at  B,  C,  D,  and  E.  The  nodes  at  D  and  E  are  needed 
to  turn  the  new  triangles  back  Into  quadrilaterals.  Following  the 
construc',.lon  of  new  nodes,  the  WNOD,  WCELL,  and  WNDND  arrays  of  adjacent 
nodes  and  cells  are  revised  to  account  for  the  new  node.  The  new  node 
(or  cell)  Is  merely  added  to  the  node  (or  cell)  list  and  provided  with 
Che  usual  properties.  None  of  these  steps  is  trivial,  but  they  can  be 
accomplished  without  modifying  the  basic  computational  algorithm:  the 
only  changes  are  those  required  at  Che  time  of  cell  and  node  creation. 

Time  Step  Control 

A  new  time  step  control  procedure  was  derived  to  prevent  the 
excessive  distortion  that  results  in  quadrilateral  cells  with  negative 
area.  The  procedure  is  based  on  the  Courant  condition  as  usual,  but 
Includes  several  beneficial  new  features. 

The  Courant  condition  requires  that  the  computational  time  step  not 
exceed  the  time  for  stress  waves  to  cross  a  cell.  Thus  the  time  step  is 

At  <  -^  (F-18) 

where  AX  Is  the  cell  dimension  and  C  is  the  sound  speed.  This  condition 
is  required  for  stability  of  wave  propagation  calculations  in  which  an 
explicit  solution  method  is  used,  as  in  HEMP. 

The  time  step  control  procedure  is  derived  below,  including 

•  Expansion  of  the  Courant  stability  condition  to  allow 
treatment  of  general,  time-dependent,  stress-strain 
relations  and  the  artificial  viscous  stress. 
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•  Nocin  added  to  permit  separation  of  cells  at  nodes 
X  Nodes  added  to  maintain  quadrilateral  cells 

MA-7893-12B 

FIGURE  F10  POSSIBLE  CONFIGURATION  OF  CELLS  AND  NODES  ALONG  A  WANDERING 
SLIDE  LINE  REPRESENTING  A  GROWING  CRACK  OR  SHEAR  BAND 


•  Addition  of  a  convective  term,  which  Is  especially  Important 
In  regions  of  large  distortion  and  wherever  particle  velocity 
Is  comparable  to  sound  speed. 

•  Extension  of  the  Courant  condition  and  convective  correction 
from  the  usual  one-dlmenslonal  form  to  a  two-dimensional  form 

•  Placement  of  the  time  step  calculation  at  a  optlmlum  location 
in  the  flow  of  the  hydrodynamic  computations  to  make  the  time 
step  sensitive  to  the  current  cell  configuration. 

One  computational  result  with  the  new  time  step  control  Is  given  below 
to  Illustrate  its  effectiveness. 


Courant  Stability  Condition 

For  one-dlmenslonal  hydrodynamic  calculations,  a  fairly  full 
investigation  of  the  Courant  condition  has  been  undertaken  by  Rlchtmyer 
and  Morton^^  and  by  Herrmann  et  al.^^  However,  a  further  development 
is  required  for  a  situation  In  which  more  general,  time-dependent, 
stress-strain  relations  are  Included  in  one-  or  two-dimensional  calculations. 

Here  a  new  concept  for  stability  is  developed  based  on  the  hypothesis 
that  all  forms  of  stability  depend  simply  on  an  effective  sound  speed. 

The  addition  of  artificial  viscosity  or  other  types  of  viscosity  tend 
to  require  shorter  time  steps  In  the  calculations  because  they  Increase 
the  effective  sound  speed.  This  effective  sound  speed  is  related  to  an 
effective  modulus  for  the  material  as  follows: 


M 

e 


2  .  , 

AP  +  2Q  3  ^°1 
■  V  ^i~  3V 


(F-19) 


where 

AF  •  the  change  In  pressure  at  a  given  cell  during  the 
current  time  Increment 

Q  ■  the  artificial  viscous  stress 

AV 

~  «  the  relative  change  in  volume 
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Aa'  ■  the  change  In  devlator  stress 

•  the  strain  In  the  direction. 

First  a  calculation  Is  made  for  the  one-dlmenslonal  problem  to  com¬ 
pare  the  results  from  the  hypothesis  above  with  that  of  Che  more  usual 
approach  followed  by  Herrmann  et  al.^^.  The  linear  and  quadratic  forms 
for  Che  artificial  viscosity  are 

Q  .  C^41!C||  .  pcj^i  (F-20) 

where 

C, ,C.  *  coefficients  of  Che  linear  and  quadratic  terms  for 

^  artificial  viscosity 

AX  «  cell  thickness 
C  ■  sound  speed  for  elastic  waves 
p  ■  density 
c  •  time 


This  expression  Is  used  only  if  Q  Is  positive,  that  is,  for  compression; 
Q  Is  set  to  0  for  rarefaction  waves.  Equation  (F-20)  is  transformed  by 
using  the  Identities  relating  density  and  specific  volume,  and  Che  fol¬ 
lowing  equation  describing  mass  conservation 


_1  ^  9u 

p  3t  3x 


(F-21) 


Equation  (F-20)  Chen  Cakes  the  following  form  in  which  the  absolute  value 
sign  Is  used  with  the  change  In  particle  velocity  so  that  Q  will  be  posi¬ 
tive  In  compression. 

Q  -  -  T  P  If  tc^l^ul  +  C^C]  (F-22) 
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where 


V  «  the  specific  volume 

Ac  «  Che  time  Increment  used  In  the  computaclons. 

From  Eq.  (F-22),  a  modulus  can  be  defined  as  In  Eq.  (F-19),  chat  Is, 

"q  ■  —k'  ■  *  V 

2V 


The  factor  2  Is  Inserted  because  the  Q  Is  computed  at  Ac/2,  providing 
twice  the  stiffness  that  would  occur  If  were  averaged  over  the 
Interval.  Then  the  effective  sound  speed  which  will  govern  stability 
of  the  calculations.  Is  from  Eqs.  (F-19)  and  (F-23) 
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p  p 
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^  AV  «  AV  “7  AV\ 

V  V  ‘^(^l-3v) 
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+  2(CqiAuI  +  C^C)  ^ 


(F-24a) 

(F-24b) 

(F-24c) 


where 

K  ■  bulk  modulus 
y  -  shear  modulus 

C  ~  the  usual  sound  speed  based  on  these  moduli. 


For  stable  calculations,  the  Courant-Frledrlchs-Lewy  condition  as 
follows  Is  used  to  determine  the  time  Increment  permitted 


where  is  usuelly  Interpreted  to  be  the  local  sound  speedy  but  here  It 
will  be  taken  as  the  effective  sound  speed  C  given  by  Eq.  (F-~2A).  In¬ 
sertion  of  Eq.  (F-2A)  into  relation  (F-25)  gives 
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+  (C^lAul  +  C^C)  II 


(F-26) 


When  this  quadratic  Is  solved  for  the  maxlniun  time  increment,  the 
criterion  is  found  to  be 
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+  c,V  -  1 
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V  J 

(F-27) 


This  expression  is  identical  to  that  found  on  page  73  of  reference  16. 
Hence  the  new  concept  does  appear  to  give  the  correct  stability  criterion 
for  one-dimensional  flow. 


For  multidimensional  problems,  the  last  term  in  Eq.  (F-24a)  is  taken 
as  the  maximum  of  the  values  in  the  three  principal  directions.  An  al¬ 
ternate,  and  simpler,  course  is  to  use  the  last  term  from  Eq.  (F-2Ab) 
with  y  equal  to  the  largest  shear  modulus.  Then  the  effective  sound 
speed  is 
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(F-28) 


This  criterion  is  sufficiently  general  for  stress  relaxation  and  other 
time-dependent  stress-strain  relations. 


Particle  Velocity  (Convective)  Effect  in  One  Dimension 

In  the  time  step  algorithm  used  here,  we  Introduce  a  convective 
term  to  account  for  the  fact  Chat  AX  changes  size  during  the  time  step. 

In  small-distortion  problems,  this  correction  is  unnecessarily  meticulous, 
but  it  is  vital  in  preventing  the  cell  tangling  that  occurs  in  large- 
distortion  problems.  Cell  tangling  occurs  when  any  node  of  a  quadrilateral 
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moves  through  an  opposite  side.  Such  motion  can  have  no  physical 
meaning  in  a  continuum,  so  the  physical  relevance  of  the  calculation 
ceases  at  that  time. 

The  convection  calculation  begins  with  the  Courant  condition 
(Eq.  F-18) ■  Here  is  rewritten  as  the  average  cell  dimension  during 
the  time  step.  That  is, 

«n^  1  ,  v^  ^n^  1 

dX  *2  2  ■  *1  ■  *1  (F-29) 

^t  -  -  -  - - 

The  cell  extends  from  Xj^  to  X2.  The  superscripts  n  and  n+1  refer  to 
time  increments  such  that  t*'^^  -  t”  ■  At.  The  locations  at  the  n+1 
time  are  replaced  by  using  the  velocity  relation: 


X""^^  -  X"  +  u”  *  At  (f-30) 

where  u”  ^  ^  is  the  particle  velocity  at  t*'  +  At/2.  Then  Eq.  (F-29) 
becomes 


At 


X^  -  xj  +  I  (u;  -  up  i)  At 


Combining  Che  At  terms,  we  obtain 


At  ■ 


X°  +  x” 
2  ^1 


(F-31) 


In  the  denominator  this  expression  for  Ac  contains  Che  convective 
correction  AU/C,  which  characterizes  Che  cell  compression  during  the 
time  interval. 


Crossing  Time  and  Convective  Effect  in  Two  Dimensions 

A  two-dimensional  geometry  is  examined  here  to  determine  the 
effective  cell  dimension  AX  for  Eqs.  (F-18)  or  (F-31)  and  the  convective 
correction  AU/C.  For  a  quadrilateral  cell,  Che  minimum  distances  for  waves 
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to  cr«v«l  4re  frcm  a  node  to  an  opposite  side.  The  eight  possible 
mlnlfflua  paths  for  a  general  quadrilateral  are  shown  In  Figure  Fll. 

Some  perpendiculars  from  the  nodes  to  the  sides,  such  as  a  and  b,  lie 
outside  the  cell;  In  this  case  the  minimum  path  Is  the  side  1-2.  If 
the  Interior  angle  at  a  node  such  as  1  Is  greater  than  90”,  then 
perpendiculars  from  adjacent  nodes  <2  and  A)  to  sides  intersecting  at 
1  will  lie  outside  Che  cell  (perpendiculars  b  and  h) .  The  paths  are 
therefore  computed  by  calculating  the  Interior  angle  for  each  node 
and  Chen  the  two  paths  pertaining  to  chat  node.  After  the  paths  are 
determined,  the  relative  velocity  along  each  path  Is  found.  Then  the 
time  step  la  computed  from  Eq.  <F-31). 

The  computation  of  the  minimum  crossing  distance  AX  and  Che 
convective  correction  AU/C  requires  calculation  of  the  following 
quantities: 

•  The  lengths  of  Che  sides  of  the  quadrilateral. 

•  The  angles  between  Che  sides  of  the  quadrilateral  at 
each  Joint. 

•  The  relative  velocity  In  Che  direction  of  a  side. 

a  The  perpendicular  distance  from  a  node  to  a  line. 

•  The  relative  velocity  along  the  perpendicular  distance 
between  a  node  and  a  line. 

The  derivations,  which  are  not  Included  here,  are  all  readily  performed 
by  cresting  the  quadrilateral  sides  and  other  lines  as  vectors  and  using 
dot  or  cross  products. 

Placement  of  the  Time  Step  Calculation 

The  location  of  the  time  step  calculation  In  the  flow  of  the  wave 
propagation  calculations  Is  critical  for  obtaining  a  time  Increment 
chat  Is  a  good  estimate  for  the  next  computational  cycle.  Here,  an 
optimum  location  Is  found,  and  a  test  showing  the  importance  of  this 
location  la  described. 
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Th«  six  basic  scaps  In  Cha  wava  propagation  calculations  are  assembled 
hare  for  datarmlnlng  whars  to  put  the  time  step  calculation. 

(call  loop) 


1.  o"  -  m/av" 
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x""^, 

u"' 

-1/2 
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(F-16) 

(F-32) 

(F-15) 

(F-17) 


(node  loop) 


ikl  ^""ik 


5  u  .  n-1/2  At^i:  ^ 

“ij  ij  l.k 

hi  *  '‘d  * 


Tl^l  /? 

The  desired  time  step  is  At  from  Eq.  (F-31),  or 

yn+3./2j 


(F-12) 

(F-7) 

(F-31a) 


Now  we  note  that  is  required  for  the  velocity  calculation  in 

step  5,  because  At"  is  an  average  of  At"'*’^^^  and  At""^/^,  However,  the 
un+1/2  required  for  the  time  step  is  available  only  after  step  5. 

To  resolve  the  conflict  in  priority  over  calculations  of  Ac"'*'!/^ 
and  un+1/2,  we  separated  step  5  into  parts.  First,  the  right-hand  side 
of  (F-12)  is  evaluated  and  an  astinate  of  U  is  made  based  on 


.n-1/2 


n-1/2  5^ 

“ij  M,  ’^ikl  ^1 

j  XfK 


(F-33) 


This  velocity  estlsate  is  used  in  Eq.  (F-31)  to  obtain  the  tine  step. 
Finally  ia  computed  from  Eq.  (F-12).  In  the  code  these  calculations 

of  and  At^^^^  are  made  in  an  additional  loon  over  all  cells: 

this  loop  occurs  between  the  main  cell  and  node  loops. 
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Thi*  tsclmaclon  procedure  for  obceinlng  is  satisfactory 

because  we  need  only  correct  limits  on  the  velocity  to  obtain  safety 
from  the  time  step  calculation.  If  U.?**  represent  an  expansion,  the 


13 

velocity  terms  in  the  calculation  of  Ac 


n+1/2 


are  neglected 
n+1/2 


If  U 


est 


ij 


will  be  less  chan 
Hence  In  the  worst  case. 


represents  a  change  coward  compression,  Ac 
AC°~^^^,  and  the  change  is  overestimeted 
the  calculation  of  At””*"^^^  will  give  a  slight  underestimate. 

To  test  the  advantages  of  the  placement  of  the  time  step  celculatlon, 
a  projectile  Impact  simulation  was  made  chat  had  also  been  made  with 
TROTT.  In  the  TROTT  run,  the  time  step  calculation  with  the  convective 
correction  was  made  following  step  4  and  was  therefore  besed  on 
The  TROTT  calculations  halted  at  about  16  us  because  of  cell  tangling. 

The  corresponding  Composite  HEMP  calculation  was  halted  at  37  us  with 
the  cell  pattern  aho%m  in  Figure  S3  of  Section  Vtl.  The  cells  are  very 
distorted,  but  tangling  or  formation  of  negative  areas  has  not  occurred. 

An  examination  of  the  time  steps  for  Che  two  calculations  showed  no 
significant  differences.  So  the  new  procedure  does  not  require  excessively 
small  time  steps  to  produce  the  added  acabilicy.  Hence  Che  placement 
of  Che  time  step  calculation  In  Composite  HEMP  appears  to  be  a  clear 
improvement. 


Rezonlng 

A  subroutine  REZONH  has  been  developed  for  rezonlng.  The  rezoner 
operates  on  a  rectangular  grid  of  cells  such  as  chat  used  in  finite- 
difference  wave  propagation  computer  programs.  As  the  cells  are  distorted 
by  material  motion,  the  rezoner  modifies  the  cell  coordinate  positions 
to  produce  an  approximately  rectangular  grid.  REZONH  is  a  continuous 
rezoner;  it  makes  small  changes  in  some  of  the  coordinate  locations  at 
regular  intervals,  rather  chan  moving  all  coordinate  locations  at  a  few 
discrete  times  during  a  wave  propagation  problem.  From  our  review  of 
rezoners  used  with  STEALTH,  TOODY,  and  HEMP,  the  continuous  rezoner  was 
selected  as  being  more  reliable,  requiring  less  programmer  intervention 
in  running  a  problem,  and  being  easier  to  write.  (The  disadvantage  of 
Che  continuous  rezoner  is  that  it  tends  to  smear  properties  more  chan  a 
discrete  rezoner  does.) 
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The  rezoner  operates  on  Che  entire  grid,  selecting  alternate  nodes 
for  possible  motion.  For  each  selected  node,  the  distances  to  adjacent 
nodes  are  tested.  If  the  distances  are  not  approximately  equal,  the 
node  is  moved  to  make  the  distances  more  nearly  equal.  In  the  current 
rezoning,  the  moved  node  may  be  in  the  center  of  four  cells  or  on  a 
boundary  and  adjacent  to  one  or  two  cells.  Rezoning  cannot  yet  occur  on 
boundaries  between  materials  unless  that  boundary  is  also  a  slide  line. 

A  rezoning  pattern  requires  four  standard  computational  cycles. 
Following  the  first  cycle,  all  odd-numbered  nodes  are  rezoned.  After 
the  second  cycle,  even-numbered  nodes  are  rezoned.  After  the  third 
and  fourth  cycles,  odd-and  even-numbered  nodes,  respectively,  are 
rezoned  starting  with  the  last  node  and  proceeding  towards  the  first. 

Described  below  are  the  input  controls,  the  coordinate  moving 
algorithm,  the  procedure  for  allocating  stresses  and  other  properties 
to  "he  newly  configured  cells,  and  the  system  for  computing  velocities 
so  that  momentum  la  conserved.  Because  each  of  these  procedures  differs 
from  those  used  in  other  codes,  the  subroutine  represents  a  new  approach. 

The  REZONK  subroutine  requires  input  of  controls  for  the  frequency 
and  degree  of  rezoning  and  to  indicate  any  nodes  for  which  rezoning  is 
not  permitted.  The  input  parameters  are  NFREQ,  REZMIN,  and  REZWT.  NFREQ 
is  the  rezoning  frequency.  For  NFREQ  equal  to  10,  the  standard  rezoning 
pattern  occurs  following  the  first  four  computational  cycles,  there  is 
no  rezoning  for  six  cycles,  and  then  the  sequence  is  repeated.  REZMIN 
is  the  minimum  amount  of  relative  motion  that  the  rezoner  will 
consider;  for  lesser  rezoning  motions,  the  node  is  skipped  and  not  rezoned. 
REZWT  is  the  weighting  factor  determining  how  rapidly  the  nodes  are  moved 
toward  the  equal-distance  position  chosen  by  the  rezoner. 

For  each  node,  two  rezoning  steps  occur.  First,  the  nodal  motion 
Is  computed  and  then  the  cell  masses  and  other  properties  are 
calculated  for  the  new  cell  configuration.  The  equal  nodal  distance 
criteria  result  in  indicated  motion  AX  and  AY.  Next  AX  and  AY  are 
compared  with  REZMIN  to  determine  whether  the  rezoning  should  occur.  If 
the  criterion  is  exceeded,  the  node  is  moved  REZWT*AX  and  REZWT*AY  in 
the  X  and  Y  directions. 
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Moving  the  coordinate  points  necessitates  the  redefinition  of  the 

size  and  other  properties  of  all  the  contiguous  cells.  First,  the 

masses  Z  of  the  new  cells  are  computed  by  an  exact  sectioning  of  the 
n 

volumes  of  the  old  cells  based  on  the  new  and  old  coordinate  lines. 


Z  (i) 
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Z -(i) 
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'ij 


(F-34) 


and 


Cij  -  (F.35) 

where  Z  (i)  is  the  mass  of  the  illi  old  cell 
o 

is  the  density  of  the  old  cell 
is  the  overlap  of  the  i^  new  cell  and  the  old  cell. 

The  cell  index  k  is  determined  by  the  nature  of  the  overlap.  If  the 
new  i  cell  overlaps  the  old  j  cell,  then  is  positive  and  Pj^  •  . 

If  the  overlap  is  reversed,  so  that  the  old  i  cell  overlaps  the  new  j 
cell,  then  is  negative  and  pj^  - 

Mew  cell  densities  are  computed  as  the  ratio  of  ma's  to  cell 
volume,  whereas  cell  area  is  computed  directly.  By  a  mass-weighting 
calculation,  all  other  cell-centered  quantities  (e.g.,  stresses, 
internal  energy,  yield  strength,  fracture  state  quantities)  are  computed 
for  the  new  cells. 


For  example,  the  pressure  is  computed  from 

P„(i)  -  I  Cj  P^{j)  /  Z  Cj  (F-36) 

are  new  and  old  cell  pressures,  and 

as  above  if  mass  is  added  to  the  i—  cell  from  the 
1 cell 

0  if  mass  is  removed  from  the  cell 
“  ^ik  ^  ■  J* 


where  P  and 
n 


'ij 
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With  the  special  definition  of  the  weighting  factor  ,  the  cell 
quantities  are  unaltered  if  portions  of  the  cell  are  lost  in  rezoning, 
but  an  averaging  occurs  if  mass  is  added  to  the  cell.  Note  that 

Z„(l)  .  E  Cj  (F-37, 


The  velocity  of  the  new  coordinate  point  is  computed  by  requiring 
an  exact  conservation  of  linear  momentum.  For  this  calculation,  each 
of  the  four  (or  two  for  a  coordinate  point  along  an  edge)  adjacent 
cells  are  partitioned  into  four  submasses  z^,  as  they  would  be  for  a 
normal  momentum  calculation.  Each  of  the  submasses  is  assigned  to  its 
nearest  coordinate  point.  Then  the  momentum  balance  equation  is 

I  Vok  "  L  “k  "nk  ‘‘c^  "nk 


where 


'He' 


z 


ok’ 


u  are  the  velocities  at  all  the  old  coordinate  points,  and 
c 

the  velocity  (to  be  computed)  at  this  new  coordinate  point 

z  ,  are  old  and  new  submasses, 
nk 


The  first  sunmatlon  runs  over  all  16  submasses,  but  the  second  skips 
Che  four  submasses  around  the  central  (new)  coordinate.  The  last  summation 
is  for  the  four  central  masses.  With  this  equation,  the  velocity  u^ 
can  be  computed. 

The  rezoner  has  been  written  and  combined  with  Composite  HEMP  and 
a  calculation  of  a  projectile  impact  has  been  attempted.  The  impact  is 
one  in  which  large  shear  distortions  occur  so  that  rezoning  is  needed. 
REZONH  will  require  some  debugging  before  this  calculation  can  be 
completed. 


Slide  Lines 

The  slide  line  is  constructed  to  permit  relative  motion  between 
adjacent  materials.  A  slide  line  is  usually  required  to  represent  the 
interface  between  a  penetrator  and  target:  the  slide  allows  gradual 
contact  of  the  penetrator  and  target  surfaces  and  permits  the  penetrator 
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to  slip  along  the  target.  The  slide  line  also  represents  the  foimation 
and  relative  motion  of  a  distinct  slip  surface  within  the  target. 

A  single  slide  line  treatment  has  been  derived  and  written  into  the 
CYCLE  subroutine  of  Composite  HEMP.  The  coding  permits  a  variety  of 
slip  conditions:  free  sliding  with  no  shear  stress,  stick  conditions, 
a  shear  stress  with  a  Coulomb  friction  limit,  and  a  viscous  shear  stress. 
The  latter  two  may  be  combined  so  that  Coulomb  friction  determines  the 
shear  stress  above  which  sliding  can  occur.  For  sliding,  the  shear 
stress  is  given  by  a  viscous  expression.  In  any  case,  separation  and 
impact  can  occur  on  the  slide  line. 

The  slide  line  strategy  generally  follows  that  in  TOODY.  Of  the 
two  surfaces  making  up  the  slide  line,  one  is  designated  the  master  and 
the  other  the  slave.  As  in  TOODY,  the  treatment  of  the  two  surfaces 
is  nearly  identical  so  it  is  not  important  which  is  designated  the  master. 
The  TOODY  technique  also  provides  a  satisfactory  transfer  of  shear 
stress  across  the  slide. 

The  slide  calculations  occur  at  two  locations  In  CYCLE:  during 
the  momentum  computation  and  during  the  boundary  condition  treatment. 

For  the  momentum  calculation,  the  forces  applied  by  the  slave  nodes  on 
the  master  nodes  and  the  forces  of  the  masters  on  the  slaves  are 
computed.  Then  during  the  boundary  calculations,  the  slave  nodes  are 
forced  to  lie  on  the  master  surface.  Also  in  this  last  step,  if  any 
momentum  is  lost  or  gained,  it  is  stored  for  use  during  the  next  momentum 
calculation.  To  avoid  the  "chatter"  of  impact  and  separation  that  often 
occurs  on  slide  lines.  Impact  is  permitted  only  on  the  boundary  portion 
of  the  cycle  and  separation  only  on  the  momentum  portion.  For  separation, 
a  tensile  stress  criterion  must  be  met. 

The  slide  line  logic  has  been  constructed  with  continuous  rezoning 
in  mind.  With  rezonlng,  nodes  may  be  shifted  by  arbitrary  amounts,  may 


* 

Single  means  nonintersecting.  Several  such  single  slide  lines  can  occur 
in  a  problem. 
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penetrate  the  slide  line,  or  separate  from  it.  Therefore,  the  first 
step  in  each  slide  calculation  is  to  locate  the  point  with  respect  to 
the  nodes  on  the  mating  surface.  This  location  process  is  the  most 
difficult  of  the  several  slide  operations. 

Approaches  have  been  formulated  for  treating  intersecting  slide 
lines  and  wandering  slide  lines,  but  these  ideas  have  not  been  completely 
developed  or  coded.  The  intersecting  slide  lines  require  an  overseer 
section  in  CYCLE  to  assess  the  intersection  conditions  and  keep  track 
of  which  portions  of  the  bounding  surfaces  are  master,  which  are  slave, 
and  which  have  ocher  conditions.  This  overseer  can  be  added  later 
without  major  changes  to  the  existing  coding. 

The  wandering  or  unzipping  slide  lines  will  permit  a  plug  to  be 
formed  in  a  target.  In  our  plan,  the  direction  and  speed  with  which  the 
slide  line  is  formed  will  depend  on  shear  band  damage  in  adjacent  cells. 
The  slide  may  proceed  either  along  cell  boundaries  or  diagonally  across 
cells  as  Illustrated  in  Figure  FIO.  The  computations  associated  with 
construction  of  the  slide  line  and  generation  of  new  nodes  and  cells  is 
greatly  facilitated  by  the  finite-element  numbering  technique.  Following 
Che  renumbering  associated  with  unzipping  Che  slide,  the  new  cells  and 
nodes  can  be  treated  by  the  standard  algorithms  used  for  all  other  cells 
and  nodes. 

Currently,  only  the  single  slide  treatment  has  been  written,  and 
it  has  been  tested  only  for  Impact  of  a  projectile  onto  a  cylindrical 
wall.  (The  case  Is  illustrated  in  Figure  53  of  Section  VII.)  Both 
the  wall  and  the  projectile  surface  are  represented  by  a  series  of  node 
points.  Thus,  the  test  verifies  the  methods  used  to  locate  adjacent 
nodes  on  the  mating  surface,  to  detect  penetration  of  the  master  line 
by  slave  nodes,  and  to  slide  nodes  along  the  master  line. 

Further  tests  are  planned  for  the  single  slide  treatment,  and  then 
the  intersecting  and  wandering  slide  development  will  begin. 
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